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ABSTRACT 

We present the implications for cosmic inflation of the Planck measurements of the cosmic microwave background (CMB) anisotropies in both 
temperature and polarization based on the full Planck survey, which includes more than twice the integration time of the nominal survey used 
for the 2013 Release papers. The Planck full mission temperature data and a first release of polarization data on large angular scales measure the 
spectral index of curvature perturbations to be n^ = 0.968 ± 0.006 and tightly constrain its scale dependence to dn^^ldlnk = -0.003 ± 0.007 when 
combined with the Planck lensing likelihood. When the Planck high-^ polarization data is included, the results are consistent and uncertainties are 
further reduced. The upper bound on the tensor-to-scalar ratio is rg oo 2 < 0.11 (95 % CL). This upper limit is consistent with the 6-mode polarization 
constraint r < 0.12 (95 % CL) obtained from a joint analysis of the BICEP2/Keck Array and Planck data. These results imply that V((^) oc 0’ and 
natural inflation are now disfavoured compared to models predicting a smaller tensor-to-scalar ratio, such as 6’ inflation. We search for several 
physically motivated deviations from a simple power-law spectrum of curvature perturbations, including those motivated by a reconstruction of 
the inflaton potential not relying on the slow-roll approximation. We find that such models are not preferred, either according to a Bayesian model 
comparison or according to a frequentist simulation-based analysis. Three independent methods reconstructing the primordial power spectrum 
consistently recover a featureless and smooth 'P«(A:) over the range of scales 0.008 Mpc“* < k < 0.1 Mpc“*. At large scales, each method finds 
deviations from a power law, connected to a deficit at multipoles { « 20^0 in the temperature power spectrum, but at an uncompelling statistical 
significance owing to the large cosmic variance present at these multipoles. By combining power spectrum and non-Gaussianity bounds, we 
constrain models with generalized Lagrangians, including Galileon models and axion monodromy models. The Planck data are consistent with 
adiabatic primordial perturbations, and the estimated values for the parameters of the base ACDM model are not significantly altered when more 
general initial conditions are admitted. In correlated mixed adiabatic and isocurvature models, the 95 % CL upper bound for the non-adiabatic 
contribution to the observed CMB temperature variance is ItTnon-adil < 1.9 %, 4.0 %, and 2.9 % for cold dark matter (CDM), neutrino density, and 
neutrino velocity isocurvature modes, respectively. We have tested inflationary models producing an anisotropic modulation of the primordial 
curvature power spectrum finding that the dipolar modulation in the CMB temperature field induced by a CDM isocurvature perturbation is not 
preferred at a statistically significant level. We also establish tight constraints on a possible quadrupolar modulation of the curvature perturbation. 
These results are consistent with the Planck 2013 analysis based on the nominal mission data and further constrain slow-roll single-field inflationary 
models, as expected from the increased precision of Planck data using the full set of observations. 

Key words. Cosmology: theory - early Universe - inflation 




1. Introduction 

The precise measurements by Planck^ of the cosmic microwave 
background (CMB) anisotropies covering the entire sky and over 
a broad range of scales, from the largest visible down to a reso¬ 
lution of approximately 5', provide a powerful probe of cosmic 
inflation, as detailed in the Planck 2013 inflation paper (Planck 
Collaboration XXII, 2014, hereafter PCI13). In the 2013 re¬ 
sults, the robust detection of the departure of the scalar spec¬ 
tral index from exact scale invariance, i.e., < 1, at more 

than 5 cr confidence, as well as the lack of the observation of 
any statistically significant running of the spectral index, were 
found to be consistent with simple slow-roll models of inflation. 
Single-field inflationary models with a standard kinetic term 
were also found compatible with the new tight upper bounds 
on the primordial non-Gaussianity parameters /nl reported in 
Planck Collaboration XXVI (2014). No evidence of isocurva¬ 
ture perturbations as generated in multi-field inflationary models 
(PCI13) or by cosmic strings or topological defects was found 
(Planck Collaboration XXV, 2014). The Planck 2013 results 
overall favoured the simplest inflationary models. However, we 
noted an amplitude deficit for ^ < 40 whose statistical sig¬ 
nificance relative to the six-parameter base A cold dark matter 
(ACDM) model is only about 2 cr, as well as other anomalies on 
large angular scales but also without compelling statistical sig¬ 
nificance. The constraint on the tensor-to-scalar ratio, r < 0.12 at 
95 % CL, inferred from the temperature power spectrum alone, 
combined with the determination of n^, suggested models with 
concave potentials. 

This paper updates the implications for inflation in the light 
of the Planck full mission temperature and polarization data. The 
Planck 2013 cosmology results included only the nominal mis¬ 
sion, comprising the first 14 months of the data taken, and used 
only the temperature data. However, the full mission includes 
the full 29 months of scientific data taken by the cryogenically 
cooled high frequency instrument (HFI) (which ended when the 
^He/'^He supply for the final stage of the cooling chain ran out) 
and the approximately four years of data taken by the low fre¬ 
quency instrument (LFI), which covered a longer period than the 
HFI because the LFI did not rely on cooling down to 100 mK for 
its operation. For a detailed discussion of the new likelihood and 
a comparison with the 2013 likelihood, we refer the reader to 
Planck Collaboration XI (2015) and Planck Collaboration XIII 
(2015), but we mention here some highlights of the differences 
between the 2013 and 2015 data processing and likelihoods: (1) 
Improvements in the data processing such as beam characteriza¬ 
tion and absolute calibration at each frequency result in a better 
removal of systematic effects. (2) The 2015 temperature high-^ 
likelihood uses half-mission cross-power spectra over more of 
the sky, owing to less aggressive Galactic cuts. The use of po¬ 
larization information in the 2015 likelihood release contributes 
to the constraining power of Planck in two principal ways: (1) 
The measurement of the F-mode polarization at large angular 
scales (presently based on the 70 GHz channel) constrains the 
reionization optical depth, t, independently of other estimates 
using ancillary data; and (2) the measurement of the TE and 

* Corresponding authors: Martin Bucher, bucher@apc.univ-paris7.fr; 
Fabio Finelli, finelli@iasfbo.inaf.it 
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EE spectra at f > 30 at the same frequencies used for the TT 
spectra (100, 143, and 217 GHz) helps break parameter degen¬ 
eracies, particularly for extended cosmological models (beyond 
the baseline six-parameter model). A full analysis of the Planck 
low-/” polarization is still in progress and will be the subject of 
another forthcoming set of Planck publications. 

The Planck 2013 results have sparked a revival of interest 
in several aspects of inflationary models. We mention here a 
few examples without the ambition to be exhaustive. A lively 
debate arose on the conceptual problems of some of the infla¬ 
tionary models favoured by the Planck 2013 data (Ijjas et al., 
2013; Guth et al., 2014; Linde, 2014; Ijjas et al., 2014). The 
interest in the inflationary model originally proposed by 
Starobinsky (1980) increased, since its predictions for cosmo¬ 
logical fluctuations (Starobinsky, 1980; Mukhanov & Chibisov, 
1981) are compatible with the Planck 2013 results (PCI13). It 
has been shown that supergravity motivates a potential similar to 
the Einstein gravity conformal representation of the inflation¬ 
ary model in different contexts (Ellis et al., 2013a,b; Buchmiiller 
et al., 2013; Earakos et al., 2013; Eerrara et al., 2013b). A sim¬ 
ilar potential can also be generated by spontaneous breaking of 
conformal symmetry (Kallosh & Linde, 2013). 

The constraining power of Planck also motivated a com¬ 
parison between large numbers of inflationary models (Martin 
et al., 2014) and stimulated different perspectives on how best 
to compare theoretical inflationary predictions with observations 
based on the parameterized dependence of the Hubble parameter 
on the scale factor during inflation (Mukhanov, 2013; Binetruy 
et al., 2014; Garcia-Bellido & Roest, 2014). The interpretation of 
the asymmetries on large angular scales (Planck Collaboration 
XXIII, 2014) also prompted a re-analysis of the primordial 
dipole modulation (Lyth, 2013; Liddle & Cortes, 2013; Kanno 
et al., 2013) of curvature perturbations during inflation. 

Another recent development has been the renewed interest in 
possible tensor modes generated during inflation, sparked by the 
BICEP2 results (BICEP2 Collaboration, 2014a,b). The BICEP2 
team suggested that the B-mode polarization signal detected at 
50 <{ < 150 at a single frequency (150 GHz) might be of pri¬ 
mordial origin. However, a crucial step in this possible interpre¬ 
tation was excluding an explanation based on polarized thermal 
dust emission from our Galaxy. The BICEP2 team put forward a 
number of models to estimate the likely contribution from dust, 
but at the time relevant observational data were lacking, and this 
modelling involved a high degree of extrapolation. If dust polar¬ 
ization were negligible in the observed patch of 380 deg^, this in¬ 
terpretation would lead to a tensor-to-scalar ratio r - 0.2^q for 
a scale-invariant spectrum. A value of r ^ 0.2, as suggested by 
BICEP2 Collaboration (2014b), would have obviously changed 
the Planck 2013 perspective from which slow-roll inflationary 
models are favoured, and such a high value of r would also 
have required a strong running of the scalar spectral index, or 
some other modification from a simple power-law spectrum, to 
reconcile the contribution of gravitational waves to temperature 
anisotropies at low multipoles with the observed TT spectrum. 

The interpretation of the B-mode signal in terms of gravita¬ 
tional waves alone presented in BICEP2 Collaboration (2014b) 
was later cast in doubt by Planck measurements of dust polariza¬ 
tion at 353 GHz (Planck Collaboration Int. XIX, 2014; Planck 
Collaboration Int. XX, 2014; Planck Collaboration Int. XXI, 
2014; Planck Collaboration Int. XXII, 2014). The Planck mea¬ 
surements characterized the frequency dependence of intensity 
and polarization of the Galactic dust emission, and moreover 
showed that the polarization fraction is higher than expected in 
regions of low dust emission. With the help of the Planck mea- 
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surements of Galactic dust properties (Planck Collaboration Int. 
XIX, 2014), it was shown that the interpretation of the B-mode 
polarization signal in terms of a primordial tensor signal plus a 
lensing contribution was not statistically preferred to an expla¬ 
nation based on the expected dust signal at 150 GHz plus a lens¬ 
ing contribution (see also Flauger et al., 2014a; Mortonson & 
Seljak, 2014). Subsequently, Planck Collaboration XXX (2014) 
extrapolated the Planck B-mode power spectrum of dust polar¬ 
ization at 353 GHz over the multipole range 40 < ^ < 120 to 
150 GHz, showing that the B-mode polarization signal detected 
by BICEP2 could be entirely due to dust. 

More recently, a B1CEP2/Keck Anay-Planck (BKP) joint 
analysis (BICEP2/Keck Array and Planck Collaborations, 2015, 
herafter BKP) combined the high sensitivity B-mode maps from 
BICEP2 and Keck Array with the Planck maps at higher fre¬ 
quencies where dust emission dominates. A study of the cross¬ 
correlations of all these maps in the B1CEP2 field found the 
absence of any statistically significant evidence for primordial 
gravitational waves, setting an upper limit of r < 0.12 at 
95 % CL (BKP). Although this upper limit is numerically al¬ 
most identical to the Planck 2013 result obtained combining 
the nominal mission temperature data with WMAP polarization 
to remove parameter degeneracies (Planck Collaboration XVI, 
2014; Planck Collaboration XXII, 2014), the BKP upper bound 
is much more robust against modifications of the inflationary 
model, since B-modes are insensitive to the shape of the pre¬ 
dicted scalar anisotropy pattern. In Sect. 13 we explore how the 
recent BKP analysis constrains inflationary models. 

This paper is organized as follows. Section 2 briefly re¬ 
views the additional information on the primordial cosmologi¬ 
cal fluctuations encoded in the polarization angular power spec¬ 
trum. Section 3 describes the statistical methodology as well 
as the Planck and other likelihoods used throughout the pa¬ 
per. Sections 4 and 5 discuss the Planck 2015 constraints on 
scalar and tensor fluctuations, respectively. Section 6 is dedi¬ 
cated to constraints on the slow-roll parameters and provides a 
Bayesian comparison of selected slow-roll inflationary models. 
In Sect. 7 we reconstruct the inflaton potential and the Hubble 
parameter as a Taylor expansion of the inflaton in the observ¬ 
able range without relying on the slow-roll approximation. The 
reconstruction of the curvature perturbation power spectrum is 
presented in Sect. 8. The search for parameterized features is 
presented in Sect. 9, and combined constraints from the Planck 
2015 power spectrum and primordial non-Gaussianity derived 
in Planck Collaboration XVII (2015) are presented in Sect. 10. 
The analysis of isocurvature perturbations combined and corre¬ 
lated with curvature perturbations is presented in Sect. 11. In 
Sect. 12 we study the implications of relaxing the assumption 
of the statistical isotropy of the power spectrum of primordial 
fluctuations. We discuss two examples of anisotropic inflation in 
light of the tests of isotropy performed in Planck Collaboration 
XVI (2015). Section 14 presents some concluding remarks. 

2. What new information does poiarization 
provide? 

This section provides a short theoretical overview of the extra in¬ 
formation provided by polarization data over that of temperature 
alone. (More details can be found in White et al. (1994); Ma & 
Bertschinger (1995); Bucher (2014), and references therein.) In 
Sect. 2 of the Planck 2013 inflation paper (PCI13), we gave an 
overview of the relation between the inflationary potential and 
the three-dimensional primordial scalar and tensor power spec¬ 
tra, denoted as PKik) and Ptik), respectively. (The scalar variable 



Fig. 1. Comparison of transfer functions for the scalar and ten¬ 
sor modes. The CMB transfer functions A^^(k) and A^^(k), 
where ^ - T,E, B, define the linear transformations mapping 
the primordial scalar and tensor cosmological perturbations to 
the CMB anisotropies as seen by us on the sky today. These 
functions are plotted for two representative values of the mul¬ 
tipole number: ^ = 2 (in black) and € - bS (in red). 


7? is defined precisely in Sect. 3). We shall not repeat the discus¬ 
sion there, instead referring the reader to PCI 13 and references 
therein. 

Under the assumption of statistical isotropy, which is pre¬ 
dicted in all simple models of inflation, the two-point corre¬ 
lations of the CMB anisotropies are described by the angular 
power spectra C^^, and where £ is the multipole 
number. (See Kamionkowski et al. (1997); Zaldarriaga & Seljak 
(1997); Seljak & Zaldarriaga (1997); Hu & White (1997); Hu 
et al. (1998) and references therein for early discussions elu¬ 
cidating the role of polarization.) In principle, one could also 
envisage measuring and but in theories where parity 
symmetry is not explicitly or spontaneously broken, the expecta¬ 
tion values for these cross spectra (i.e., theoretical cross spectra) 
vanish, although the observed realizations of the cross spectra 
are not exactly zero because of cosmic variance. 

The CMB angular power spectra are related to the three- 
dimensional scalar and tensor power spectra via the transfer 
functions A^^(k) and A*^^(k), so that the contributions from 
scalar and tensor perturbations are 

f Af.s(^) Pnik) (1) 
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and 

cf £ f Kjk) A' nik), (2) 


respectively, where = T, E, B. The scalar and tensor pri¬ 
mordial perturbations are uncorrelated in the simplest models, 
so the scalar and tensor power spectra add in quadrature, mean¬ 
ing that 






(3) 


Roughly speaking, the form of the linear transformations encap¬ 
sulated in the transfer functions A^^(A:) and A^^(A:) probe the 
late-time physics, whereas the primordial power spectra V-Rik) 
and Vtik) are solely determined by the primordial universe, per¬ 
haps not so far below the Planck scale if large-held inhation 
turns out to be correct. 

To better understand this connection, it is useful to plot and 
compare the shapes of the transfer functions for representative 
values of ( and characterize their qualitative behavior. Referring 
to Fig. 1, we emphasize the following qualitative features: 


1. For the scalar mode transfer functions, of which only A^ j{k) 
and i;ik) are non-vanishing (because to linear order, a 
three-dimensional scalar mode cannot contribute to the B 
mode of the polarization), both transfer functions start to rise 
at more or less the same small values of k (due to the centrifu¬ 
gal barrier in the Bessel differential equation), but A^^(k) 
falls off much faster at large k, and thus smooths sharp fea¬ 
tures in Vuik) to a smaller extent, than A^j,(k). This means 
that polarization is more powerful for reconstructing possi¬ 
ble sharp features in the scalar primordial power spectrum 
than temperature, provided that the required signal-to-noise 
is available. 

2. For the tensor modes, A^^(k) starts rising at about the same 
small k as A* ^(k) and A® ^{k) but falls off faster with increas¬ 
ing k than A^j,(k). On the other hand, the polarization com¬ 
ponents, A^g(k) and Aj,^(k)„ have a shape completely dif¬ 
ferent from any of the other transfer functions. The shape of 
A^g(k) and A^^(k) is much wider in ln(k) than the scalar po¬ 
larization transfer function, with a variance ranging from 0.5 
to 1.0 decades. These functions exhibit several acoustic os¬ 
cillations which appear as an envelope to a rapidly varying 
carrier frequency. 


Regarding the scalar primordial cosmological perturbations, 
the power spectrum of the fs-mode polarization provides an im¬ 
portant consistency check. As we explore in Sects. 8 and 9, to 
some extent the fit of the temperature power spectrum can be 
improved by allowing a complicated form for the primordial 
power spectrum (relative to a simple power law), but the and 
Cf ^ power spectra provide an important cross-check. Moreover, 
in multi-field inflationary models, in which isocurvature modes 
may have been excited (possibily correlated amongst themselves 
as well as with the adiabatic mode), polarization information 
provides a powerful way for breaking degeneracies (see, e.g., 
Bucher etal. (2001)). 

The inability of scalar modes to generate B-mode polariza¬ 
tion (apart from the effects of lensing) has an important conse¬ 
quence. For the primordial tensor modes, polarization informa¬ 
tion, especially information concerning the B-mode polarization, 
offers powerful potential for discovery or for establishing upper 
bounds. Planck 2013 and WMAP established upper bounds on 
a possible tensor mode contribution using alone, but these 
bounds crucially relied on an assumed simple form for the scalar 


primordial power spectrum. For example, as reported in PC113, 
when a simple power law was generalized to allow for running, 
the bound on the tensor contribution degraded by approximately 
a factor of two. The new joint BICEP2/Keck Anay-Planck upper 
bound (see Sect. 13), however, is much more robust and cannot 
be avoided by postulating baroque models that alter the scale 
dependence of the scalar power spectrum. 


3. Methodology 

This section describes updates to the formalism used to describe 
cosmological models and the likelihoods used with respect to the 
Planck 2013 inflation paper (PCI13). 

3.1. Cosmological model 

The cosmological models that predict observables, such as the 
CMB anisotropies, rely on inputs specifying the conditions and 
physics at play during different epochs of the history of the 
Universe. The primordial inputs describe the power spectrum of 
the cosmological perturbations at a time when all the observable 
modes were situated outside the Hubble radius. The inputs from 
this epoch consist of the primordial power spectra, which may 
include scalar curvature perturbations, tensor perturbations, and 
possibly also isocurvature modes and their correlations. The late 
time (i.e., z < 10"^) cosmological inputs include parameters such 
as Wc, T, which determine the conditions when the 

primordial perturbations become imprinted on the CMB and also 
the evolution of the Universe between last scattering and today, 
affecting primarily the angular diameter distance. Finally, there 
is a so-called “nuisance” component, consisting of parameters 
that determine how the measured CMB spectra are contaminated 
by unsubtracted Galactic and extragalactic foreground contam¬ 
ination. The focus of this paper is on the primordial inputs and 
how they are constrained by the observed CMB anisotropy, but 
we cannot completely ignore the other non-primordial parame¬ 
ters because their presence and uncertainties must be dealt with 
in order to correctly extract the primordial information of inter¬ 
est here. 

As in PCI13, we adopt the minimal six-parameter spatially 
flat base ACDM cosmological model as our baseline for the late¬ 
time cosmology, mainly altering the primordial inputs, i.e., the 
simple power-law spectrum parameterized by the scalar ampli¬ 
tude and spectral index for the adiabatic growing mode, which 
in this minimal model is the only late-time mode excited. This 
model has four free non-primordial cosmological parameters 
(wb, Wc, 6 mc> t) (for a more detailed account of this model, we 
refer the reader to Planck Collaboration XIII (2015)). On occa¬ 
sion, this assumption will be relaxed in order to consider the im¬ 
pact of more complex alternative late-time cosmologies on our 
conclusions about inflation. Some of the commonly used cosmo¬ 
logical parameters are defined in Table 1. 

3.2. Primordial spectra of cosmological fluctuations 

In inflationary models, comoving curvature {'R) and tensor (h) 
fluctuations are amplified by the nearly exponential expansion 
from quantum vacuum fluctuations to become highly squeezed 
states resembling classical states. Formally, this quantum me¬ 
chanical phenomenon is most simply described by the evolu¬ 
tion in conformal time, rj, of the mode functions for the gauge- 
invariant inflaton fluctuation, 6<p, and for the tensor fluctuation. 
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Table 1. Primordial, baseline, and optional late-time cosmological parameters. 


Parameter Definition 


As. Scalar power spectrum amplitude (at k, = 0.05 Mpc^*) 

Us. Scalar spectral index (at k, - 0.05 Mpc^', unless otherwise stated) 

dus/dlnfc. Running of scalar spectral index (at k, = 0.05 Mpc^', unless otherwise stated) 

d^njdlnk^ . Running of running of scalar spectral index (at k, = 0.05 Mpc^*) 

r . Tensor-to-scalar power ratio (at k, = 0.05 Mpc^', unless otherwise stated) 

U)..Tensor spectrum spectral index (at k, = 0.05 Mpc~‘)_ 

tUb. B ary on density today 

tUc. Cold dark matter density today 

6mc . Approximation to the angular size of sound horizon at last scattering 

T..Thomson scattering optical depth of reionized intergalactic medium_ 

Aeff. Effective number of massive and massless neutrinos 

Sm,, . Sum of neutrino masses 

Tp. Fraction of baryonic mass in primordial helium 

SIk . Spatial curvature parameter 

Wde. Dark energy equation of state parameter (i.e., pde/Pde) (assumed constant) 


h: 


(ajk)" + 



ayk = 0 , 


(4) 


with (x,y) = (a^lH,6(p) for scalars, and (x,y) = (a ,h) for ten¬ 
sors. Here a is the scale factor, primes indicate derivatives with 
respect to 77 , and (p and H - aja am the proper time derivative 
of the inflaton and the Hubble parameter, respectively. The cur¬ 
vature fluctuation, 'R, and the inflaton fluctuation, dcp, are related 
via R - H5(pl(l>. Analytic and numerical calculations of the pre¬ 
dictions for the primordial spectra of cosmological fluctuations 
generated during inflation have reached high standards of pre¬ 
cision, which are more than adequate for our purposes, and the 
largest uncertainty in testing specific inflationary models arises 
from our lack of knowledge of the history of the Universe be¬ 
tween the end of inflation and the present time, during the so- 
called “epoch of entropy generation.” 

This paper uses three different methods to compare inflation¬ 
ary predictions with Planck data. The first method consists of a 
phenomenological parameterization of the primordial spectra of 
scalar and tensor perturbations as: 


Rnik) = 


= AJ- 


Ws-1+4 d«s/dln^ln(^//:*)+7 d^Ws/dln(fc^(ln(^/^,))^+... 




dni/dln^ln(A:/A:*)+.. 


, (5) 


, (6) 


where As (At) is the scalar (tensor) amplitude and (wt), 
dns/dlnk(d«t/dlnk), and are the scalar (tensor) 

spectral index, the running of the scalar (tensor) spectral index, 
and the running of the running of the scalar spectral index, re¬ 
spectively. R denotes the comoving curvature perturbation for 
adiabatic initial conditions, h the amplitude of the two polariza¬ 
tion states (-H, x) of gravitational waves, and k, the pivot scale. 
Unless otherwise stated, the tensor-to-scalar ratio. 


" rnikS 


(7) 


is fixed to which is the relation that holds when inflation is 

driven by a single slow-rolling scalar field with a standard kinetic 


^ When running is considered, we fix «, = -r(2 - r/8 - Ms)/ 8 and 
drilld\nk = r(r/8 + ns - l)/8. 


term. We will use a parameterization analogous to Eq. (5) with 
no running for the power spectra of isocurvature modes and their 
correlations in Sect. 11. 

The second method exploits the analytic dependence of the 
slow-roll power spectra of primordial perturbations in Eqs. (5) 
and ( 6 ) on the values of the Hubble parameter and the hierar¬ 
chy of its time derivatives, known as the Hubble flow-functions 
(HEE): e\ - e;+i s eiUHet), with i > 1. We will use 

the analytic power spectra calculated up to second order us¬ 
ing the Green’s function method (Gong & Stewart, 2001; Leach 
et ah, 2002) (see Habib et al. 2002, Martin & Schwarz 2003, and 
Casadio et al. 2006 for alternative derivations). The spectral in¬ 
dices and the relative scale dependence in Eqs. (5) and ( 6 ) are 
given in terms of the HEEs by: 

ns - 1 = - 2ei - 62 - 2ef - (2C -H 3) ei 62 - Ceie^, (8) 
dn Jd In k = - 2ei 62 - 62 63 , (9) 

n, =-261-26^-2 (C-b 1) 61 62, (10) 

dwt/dlnk = - 26162 , (11) 

where C = ln2 + je - ^ ~ -0.7296 (je is the Euler-Mascheroni 
constant). See Appendix of PCI13 for more details. Primordial 
spectra as functions of the 6 , will be employed in Sect. 6 , 
and the expressions generalizing Eqs. ( 8 ) to (11) for a general 
Lagrangian p((p,X), where X = -g^’'5^05v0/2, will be used 
in Sect. 10. The good agreement between the first and second 
method as well as with alternative approximations of slow-roll 
spectra is illustrated in the Appendix of PCI13. 

The third method is fully numerical, suitable for models 
where the slow-roll conditions are not well satisfied and an¬ 
alytical approximations for the primordial fluctuations are not 
available. Two different numerical codes, the inflation module 
of Lesgourgues & Valkenburg (2007) as implemented in CLASS 
(Lesgourgues, 2011; Bias et al., 2011) and ModeCode (Adams 
et al., 2001; Peiris et al., 2003; Mortonson et al., 2009; Easther 
& Peiris, 2012), are used in Sects. 7 and 10, respectively.^ 

Conventions for the functions and symbols used to describe 
inflationary physics are defined in Table 2. 

3.3. Planck data 

The Planck data processing proceeding from time-ordered 
data (TOD) to maps has been improved for this 2015 re- 

^ http://class-code.net, http://modecode.org 
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Table 2. Conventions and definitions for inflation physics 


Parameter Definition 


(p . Inflaton 

. Inflaton potential 

a . Scale factor 

t . Cosmic (proper) time 

6X . Fluctuation of X 

X = dXjdt . Derivative with respect to proper time 

X’ = dXjdri . Derivative with respect to conformal time 

X^ = dXjdp . Partial derivative with respect to p 

Mpi. Reduced Planck mass {= 2.435 x 10** GeV) 

Q . Scalar perturbation variable 

. Gravitational wave amplitude of (+, x)-polarization component 

X, . X evaluated at Hubble exit during inflation of mode with wavenumber k, 

Xe. X evaluated at end of inflation 

ey = MpjV'^/2V'^. First slow-roll parameter for P(0) 

rjy = . Second slow-roll parameter for P(0) 

^y = .. Third slow-roll parameter for V(0) 

my = . Fourth slow-roll parameter for V((f)) 

ei = -HI . First Hubble hierarchy parameter 

£n+i = £nlHe„ . (n -t- l)st Hubble hierarchy parameter (where u > 1) 

N(t) = H . Number of e-folds to end of inflation 


lease in various aspects (Planck Collaboration II, 2015; Planck 
Collaboration VII, 2015). We refer the interested reader to 
Planck Collaboration II (2015) and Planck Collaboration VII 
(2015) for details, and we describe here two of these improve¬ 
ments. The absolute calibration has been improved using the 
orbital dipole and more accurate characterization of the Planck 
beams. The calibration discrepancy between Planck and WMAP 
described in Planck Collaboration XXXI (2014) for the 2013 
release has now been greatly reduced. At the time of that re¬ 
lease, a blind analysis for primordial power spectrum recon- 
straction described a broad feature at f =» 1800 in the tem¬ 
perature power spectrum, which was most prominent in the 
217x217 GHz auto-spectra (PCI13). In work done after the 
Planck 2013 data release, this feature was shown to be asso¬ 
ciated with imperfectly subtracted systematic effects associated 
with the 4-K cooler lines, which were particularly strong in the 
first survey. This systematic effect was shown to potentially lead 
to 0.5 cr shifts in the cosmological parameters, slightly increas¬ 
ing Ws and H(), similarly to the case in which the 217x217 chan¬ 
nel was excised from the likelihood (Planck Collaboration XV, 
2014; Planck Collaboration XVI, 2014). The Planck likelihood 
(Planck Collaboration XI, 2015) is based on the full mission data 
and comprises temperature and polarization data (see Fig. 2). 


Planck \ov\j-£ likelihood 

The Planck low-^ temperature-polarization likelihood uses 
foreground-cleaned LFI70 GHz polarization maps together with 
the temperature map obtained from the Planck 30 to 353 GHz 
channels by the Commander component separation algorithm 
over 94 % of the sky (see Planck Collaboration IX (2015) for fur¬ 
ther details). The Planck polarization map uses the LFI 70 GHz 
(excluding Surveys 2 and 4) low-resolution maps of Q and U 
polarization from which polarized synchrotron and thermal dust 
emission components have been removed using the LFI 30 GHz 
and HFI 353 GHz maps as templates, respectively. (See Planck 
Collaboration XI (2015) for more details.) The polarization map 
covers the 46 % of the sky outside the lowP polarization mask. 


The low-f” likelihood is pixel-based and treats the tempera¬ 
ture and polarization at the same resolution of 3?6, or HEALpix 
(Gorski et al., 2005) Wide = 16. Its multipole range extends from 
^ = 2 to f = 29 in TT, TE, EE, and BB. In the 2015 Planck pa¬ 
pers the polarization part of this likelihood is denoted as “lowP.”"* 
This Planck low-^ likelihood replaces the Planck temperature 
low-f Gibbs module combined with the WMAP 9-year low-^ 
polarization module used in the Planck 2013 cosmology papers 
(denoted by WP), which used lower resolution polarization maps 
at Wide = 8 (about 7?3). With this Planck-on\y low-^ likelihood 
module, the basic Planck results presented in this release are 
completely independent of external information. 

The Planck low multipole likelihood alone implies r = 
0.067 ± 0.022 (Planck Collaboration XI, 2015), a value smaller 
than the value inferred using the WP polarization likelihood, 
T = 0.089 + 0.013, used in the Planck 2013 papers (Planck 
Collaboration XV, 2014). See Planck Collaboration XIII (2015) 
for the important implications of this decrease in t for reioniza¬ 
tion. However, the LFI 70 GHz and WMAP polarization maps 
are in very good agreement when both are foreground-cleaned 
using the HFI 353 GHz map as a polarized dust template (see 
Planck Collaboration XI (2015) for further details). Therefore, 
it is useful to construct a noise-weighted combination to obtain 
a joint P/anckAVMAP low-resolution polarization data set, also 
described in Planck Collaboration XI (2015), using as a polariza- 


'* In this paper we use the conventions introduced in Planck 
Collaboration XIII (2015). We adopt the following labels for likeli¬ 
hoods: (i) Planck TT denotes the combination of the TT likelihood at 
multipoles £> 30 and a low-7 temperature-only likelihood based on the 
CMB map recovered with Commander; (ii) Planck TT-lowT denotes 
the TT likelihood at multipoles 7 > 30; (hi) Planck TT-l-lowP further in¬ 
cludes the Planck polarization data in the low-7 likelihood, as described 
in the main text; (iv) Planck TE denotes the likelihood at 7 > 30 using 
the TE spectrum; and (v) Planck TT,TE,EE-h1owP denotes the combi¬ 
nation of the likelihood at 7 > 30 using TT, TE, and EE spectra and 
the low-7 multipole likelihood. The label “t prior” denotes the use of a 
Gaussian prior t = 0.07 ± 0.02. The labels “lowT,P” and ”lowEB” de¬ 
note the low-7 multipole likelihood and the Q, U pixel likelihood only, 
respectively. 
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Fig. 2. Planck TT (top), high-^ TE (centre), and high-^ EE (bot¬ 
tom) angular power spectra. Here De = {({ + l)C[l{2n). 


tion mask the union of the WMAP P06 and Planck lowP polar¬ 
ization masks and keeping 74 % of the sky. The polarization part 
of the combined low-multipole likelihood is called lowP-rWP 
This combined low-multipole likelihood gives t = 0.071);g g[j 
(Planck Collaboration XI, 2015). 

Planck high-^ likelihood 

Following Planck Collaboration XV (2014), and Planck 
Collaboration XI (2015) for polarization, we use a Gaussian 
approximation for the high-f” part of the likelihood (30 < { < 
2500), so that 

- iog4c|C(e)) = i(c - cio))^ Mr\c - ao)), (i2) 

where a constant offset has been discarded. Here C is the data 
vector, C{0) is the model prediction for the parameter value vec¬ 
tor 6, and M is the covariance matrix. For the data vector, we 
use 100 GHz, 143 GHz, and 217 GHz half-mission cross-power 
spectra, avoiding the Galactic plane as well as the brightest point 


sources and the regions where the CO emission is the strongest. 
We retain 66 % of the sky for 100 GHz, 57 % for 143 GHz, and 
47 % for 217 GHz for the T masks, and respectively 70 %, 50 %, 
and 41 % for the Q, U masks. Following Planck Collaboration 
XXX (2014), we do not mask for any other Galactic polarized 
emission. All the spectra are corrected for the beam and pixel 
window functions using the same beam for temperature and po¬ 
larization. (For details see Planck Collaboration XI (2015).) 

The model for the cross-spectra can be written as 




c^^e) -H 


(13) 


where C‘^™*’(fl) is the CMB power spectrum, which is indepen¬ 
dent of the frequency, C^^y(0) is the foreground model contribu¬ 
tion for the cross-frequency spectrum fi x v, and c^ is the cal¬ 
ibration factor for the fj. x /u spectrum. The model for the fore¬ 
ground residuals includes the following components: Galactic 
dust, clustered CIB, tSZ, kSZ, tSZ correlations with CIB, and 
point sources, for the TT foreground modeling; and for polar¬ 
ization, only dust is included. All the components are modelled 
by smooth Ce templates with free amplitudes, which are deter¬ 
mined along with the cosmological parameters as the likelihood 
is explored. The tSZ and kSZ models are the same as in 2013 
(see Planck Collaboration XV, 2014), while the CIB and tSZ- 
CIB correlation models use the updated CIB models described 
in Planck Collaboration XXX (2014). The point source contam¬ 
ination is modelled as Poisson noise with an independent am¬ 
plitude for each frequency pair. Finally, the dust contribution 
uses an effective smooth model measured from high frequency 
maps. Details of our dust and noise modelling can be found in 
Planck Collaboration XI (2015). The dust is the dominant fore¬ 
ground component for TT at { < 500, while the point source 
component, and for 217x217 also the CIB component, dom¬ 
inate at high {. The other foreground components are poorly 
determined by Planck. Finally, our treatment of the calibration 
factors and beam uncertainties and mismatch are described in 
Planck Collaboration XI (2015). 

The covariance matrix accounts for the correlation due to 
the mask and is computed following the equations in Planck 
Collaboration XV (2014), extended to polarization in Planck 
Collaboration XI (2015) and references therein. The hducial 
model used to compute the covariance is based on a joint ht of 
base ACDM and nuisance parameters obtained with a previous 
version of the matrix. We iterate the process until the parame¬ 
ters stop changing. For more details, see Planck Collaboration 
XI (2015). 

The joint unbinned covariance matrix is approximately of 
size 23 000 x 23 000. The memory and speed requirements for 
dealing with such a huge matrix are signihcant, so to reduce its 
size, we bin the data and the covariance matrix to compress the 
data vector size by a factor of 10. The binning uses varying bin 
width with = 5 for 29 < ^ < 100, = 9 for 99 < ^ < 2014, 

and A{ - 33 for 2013 < { < 2509, and a weighting in {{{ + 1) 
to flatten the spectrum. Where a higher resolution is desirable, 
we also use a more finely binned version (“bin3”, unbinned up 
to = 80 and At - 2 beyond that) as well as a completely 
unbinned version (“binl”). We use odd bin sizes, since for an 
azimuthally symmetric mask, the correlation between a multi¬ 
pole and its neighbours is symmetric, oscillating between posi¬ 
tive and negative values. Using the base ACDM model and sin¬ 
gle parameter classical extensions, we confirmed that the cos¬ 
mological and nuisance parameter fits with or without binning 
are indistinguishable. 
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Planck CMB bispectrum 

We use measurements of the non-Gaussianity amplitude /nl 
from the CMB bispectrum presented in Planck Collaboration 
XVII (2015). Non-Gaussianity constraints have been obtained 
using three optimal bispectrum estimators (separable template 
btting (also known as “KSW”), binned, and modal). The maps 
analysed are the Planck 2015 full mission sky maps, both in 
temperature and in E polarization, as cleaned with the four 
component separation methods SMICA, SEVEM, NILC, and 
Commander. The map is masked to remove the brightest parts 
of the Galaxy as well as the brightest point sources and cov¬ 
ers approximately 70 % of the sky. In this paper we mainly ex¬ 
ploit the joint constraints on equilateral and orthogonal non- 
Gaussianity (after removing the ISW-lensing bias), = 

-16 + 70, /°f ° = -34 + 33 from T only, and = -3.7 + 43, 

^ ^ reference, the 

constraints on local non-Gaussianity are = 2.5 + 5.7 from 
T only, and = 0.8 + 5.0 from T and E (68 % CL). Starting 
from a Gaussian /NL-likelihood, which is an accurate assump¬ 
tion in the regime of small primordial non-Gaussianity, we use 
these constraints to derive limits on the sound speed of the infla- 
ton fluctuations (or other microscopic parameters of inflationary 
models) (Planck Collaboration XXIV, 2014). The bounds on the 
sound speed for various models are then used in combination 
with Planck power spectrum data. 

Planck CMB lensing data 

Some of our analysis includes the Planck 2015 lensing likeli¬ 
hood, presented in Planck Collaboration XV (2015), which uti¬ 
lizes the non-Gaussian trispectrum induced by lensing to esti¬ 
mate the power spectrum of the lensing potential, This sig¬ 
nal is extracted using a full set of temperature- and polarization- 
based quadratic lensing estimators (Okamoto & Hu, 2003) ap¬ 
plied to the SMICA CMB map over approximately 70% of the 
sky, as described in Planck Collaboration IX (2015). We have 
used the conservative bandpower likelihood, covering multi¬ 
poles 40 < 7” < 400. This provides a measurement of the lensing 
potential power at the 40 cr level, giving a 2.5 %-accurate con¬ 
straint on the overall lensing power in this multipole range. The 
measurement of the lensing power spectrum used here is approx¬ 
imately twice as powerful as the measurement used in our pre¬ 
vious 2013 analysis (Planck Collaboration XXII, 2014; Planck 
Collaboration XVII, 2014), which used temperature-only data 
from the Planck nominal mission data set. 

3.4. Non-Planck data 
BAO data 

B ary on acoustic oscillations (BAO) are the counterpart in the 
late-time matter power spectrum of the acoustic oscillations 
seen in the CMB multipole spectrum (Eisenstein et al., 2005). 
Both originate from coherent oscillations of the photon-baryon 
plasma before these two components become decoupled at re¬ 
combination. Measuring the position of these oscillations in the 
matter power spectra at different redshifts constrains the expan¬ 
sion history of the universe after decoupling, thus removing de¬ 
generacies in the interpretation of the CMB anisotropies. 

In this paper, we combine constraints on Dy{z)lr^ (the ratio 
between the spherically averaged distance scale Dy to the effec¬ 
tive survey redshift, z, and the sound horizon, r^) inferred from 


6 dFGRS data (Beutler et al., 2011) at z = 0.106, the SDSS- 
MGS data (Ross et al., 2014) at z = 0.15, and the SDSS-DRll 
CM ASS and LOWZ data (Anderson et al., 2014) at redshifts 
z = 0.57 and 0.32. (For details see Planck Collaboration XIII 
(2015).) 

Joint BICEP2/Keck Array and Planck constraint on r 

Since the Planck temperature constraints on the tensor-to-scalar 
ratio are close to the cosmic variance limit, the inclusion of data 
sets sensitive to the expected B-mode signal of primordial gravi¬ 
tational waves is particularly useful. In this paper, we provide re¬ 
sults including the joint analysis cross-correlating BICEP2/Keck 
Array observations and Planck (BICEP2/Keck Array and Planck 
Collaborations, 2015, hereafter BKP). Combining the more sen¬ 
sitive BICEP2/Keck Array B-mode polarization maps in the ap¬ 
proximately 400 deg^ BICEP2 held with the Planck maps at 
higher frequencies where dust dominates allows a statistical 
analysis taking into account foreground contamination. Using 
BB auto- and cross-frequency spectra between BICEP2/Keck 
An-ay (150GHz) and Planck (217 and 353GHz), BKP And a 
95 % upper limit of rg.os < 0.12. 

3.5. Parameter estimation and model comparison 

Much of this paper uses a Bayesian approach to parameter es¬ 
timation, and unless otherwise specified, we assign broad top- 
hat prior probability distributions to the cosmological parame¬ 
ters listed in Table 1. We generate posterior probability distri¬ 
butions of the parameters using either the Metropolis-Hastings 
algorithm implemented in CosmoMC (Lewis & Bridle, 2002) or 
MontePython (Audren et al., 2013), the nested sampling al¬ 
gorithm MultiNest (Feroz & Hobson, 2008; Feroz et al., 
2009, 2013), or PolyChord, which combines nested sam¬ 
pling with slice sampling (Handley et al., 2015). The latter 
two also compute Bayesian evidences needed for model com¬ 
parison. Nevertheless, values are often provided as well 
(using CosmoMC’s implementation of the BOBYQA algo¬ 
rithm (Powell, 2009) for maximizing the likelihood), and other 
parts of the paper employ frequentist methods when appropriate. 

4. Constraints on the primordial spectrum of 
curvature perturbations 

One of the most important results of the Planck nominal mission 
was the determination of the departure from scale invariance for 
the spectrum of scalar perturbations at high statistical signifi¬ 
cance (Planck Collaboration XVI, 2014; Planck Collaboration 
XXII, 2014). We now update these measurements with the 
Planck full mission data in temperature and polarization. 

4.1. Tilt of the curvature power spectrum 

For the base ACDM model with a power-law power spectrum 
of curvature perturbations, the constraint on the scalar spectral 
index, «s, with the Planck full mission temperature data is 

= 0.9655 + 0.0062 (68 % CL, Planck TT-tlowP). (14) 

This result is compatible with the Planck 2013 results («s = 
0.9603 + 0.0073 (Planck Collaboration XV, 2014; Planck 
Collaboration XVI, 2014); see Fig. 3 for the accompanying 
changes in r, Tlbh^, and 0 mc)- The shift towards higher values 
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Parameter 

TT+lowP 

TT+lowP+lensing 

TT+lowP+BAO 

TT,TE,EE+lowP 

siy 

0.02222 ± 0.00023 

0.02226 ± 0.00023 

0.02226 + 0.00020 

0.02225 + 0.00016 


0.1197 + 0.0022 

0.1186 + 0.0020 

0.1190 + 0.0013 

0.1198 + 0.0015 

lOO^MC 

1.04085 + 0.00047 

1.04103 + 0.00046 

1.04095 + 0.00041 

1.04077 + 0.00032 

r 

0.078 + 0.019 

0.066 + 0.016 

0.080 + 0.017 

0.079 + 0.017 

ln(10'%) 

3.089 + 0.036 

3.062 ± 0.029 

3.093 + 0.034 

3.094 + 0.034 

Hs 

0.9655 ± 0.0062 

0.9677 ± 0.0060 

0.9673 + 0.0045 

0.9645 + 0.0049 

Ho 

67.31 + 0.96 

67.81 ± 0.92 

67.63 + 0.57 

67.27 + 0.66 


0.315 + 0.013 

0.308 + 0.012 

0.3104 + 0.0076 

0.3156 + 0.0091 


Table 3. Confidence limits on the parameters of the base ACDM model, for various combinations of Planck 2015 data, at the 68% 
confidence level. 


for Hs with respect to the nominal mission results is due to sev¬ 
eral improvements in the data processing and likelihood which 
are discussed in Sect. 3, including the removal of the 4-K cooler 
systematics. For the values of other cosmological parameters in 
the base ACDM model, see Table 3. We also provide the results 
for the base ACDM model and extended models online.^ 

When the Planck high-^ polarization is combined with tem¬ 
perature, we obtain 

Us = 0.9645 + 0.0049 (68 % CL, Planck TT,TE,EEh-1owP), 

(15) 

together with t = 0.079 + 0.017 (68% CL), which is consis¬ 
tent with the TTh-IowP results. The Planck high-f polarization 
pulls up T to a slightly higher value. When the Planck lensing 
measurement is added to the temperature data, we obtain 

Ms = 0.9677 + 0.0060 (68 % CL, Planck TT-blowP-blensing), 

(16) 

with T = 0.066 + 0.016 (68% CL). The shift towards slightly 
smaller values of the optical depth is driven by a marginal pref¬ 
erence for a smaller primordial amplitude. As, in the Planck 
lensing data (Planck Collaboration XV, 2015). Given that the 
temperature data provide a sharp constraint on the combination 
e^^Ms, a slightly lower As requires a smaller optical depth to 
reionization. 

4.2. Viability of the Harrison-Zeidovich spectrum 

Even though the estimated scalar spectral index has risen slightly 
with respect to the Planck 2013 release, the assumption of 
a Harrison-Zeidovich (HZ) scale-invariant spectrum (Harrison, 
1970; Peebles & Yu, 1970; Zeldovich, 1972) continues to be 
disfavoured (with a modest increase in significance, from 5.1 cr 
in 2013 to 5.6cr today), because the error bar on n^ has de¬ 
creased. The value of n^ inferred from the Planck 2015 tem¬ 
perature plus large-scale polarization data lies 5.6 standard de¬ 
viations away from unity (with a corresponding A;^f^ = 29.9), 
if one assumes the base ACDM late-time cosmological model. 
If we consider more general reionization models, parameterized 
by a principal component analysis (Mortonson & Hu, 2008) in¬ 
stead of T (where reionization is assumed to have occured in¬ 
stantaneously), we find tS.x^ - 14.9 for n^ - 1. Previously, 
simple one-parameter extensions of the base model, such as 
ACDM-bAeff (where Aeff is the effective number of neutrino 
flavours) or ACDM-nTp (where Yp is the primordial value of the 
helium mass fraction), could nearly reconcile the Planck tem¬ 
perature data with n^ - They now lead to A;^f^ = 7.6 and 9.3, 
respectively. Eor any of the cosmological models that we have 

^ http://www.cosmos.esa.int/web/planck/pla 


considered, the by which the HZ model is penalized with 
respect to the tilted model has increased since the 2013 analy¬ 
sis (PCI 13) thanks to the constraining power of the full mission 
temperature data. Adding Planck high-/’ polarization data further 
disfavours the HZ model; in ACDM, the;^f^ increases by 57.8, for 
general reionization we obtain A;^f^ = 41.3, and for ACDM-nAeff 
and ACDMh-Tp we find A;^f^ = 22.5 and 24.0, respectively. 

4.3. Running of the spectral index 

The running of the scalar spectral index is constrained by the 
Planck 2015 full mission temperature data to 

= -0.0084 + 0.0082 (68 % CL, Planck TT+lowP). (17) 

dmk 

The combined constraint including high-^ polarization is 

= -0.0057 + 0.0071 (68 % CL, Planck TT,TE,EE+lowP). 

dmk 

(18) 

Adding the Planck CMB lensing data to the tempera¬ 
ture data further reduces the central value for the running, 
i.e., dn^ldlnk = -0.0033 + 0.0074 (68% CL, Planck 

TT+lowP+lensing). 

The central value for the running has decreased in magn- 
tude with respect to the Planck 2013 nominal mission (Planck 
Collaboration XVI (2014) found dnjdlnk = -0.013 + 0.009; 
see Fig. 4), and the improvement of the maximum likelihood 
with respect to a power-law spectrum is smaller, x. -0.8. 
Among the different effects contributing to the decrease in the 
central value of the running with respect to the Planck 2013 re¬ 
sult, we mention a change in HFI beams at /” < 200 (Planck 
Collaboration XIII, 2015). Nevertheless, the deficit of power at 
low multipoles in the Planck 2015 temperature power spectrum 
contributes to a preference for slightly negative values of the run¬ 
ning, but with low statistical significance. 

The Planck constraints on n^ and dnsfdlnk are remarkably 
stable against the addition of the BAG likelihood. The combina¬ 
tion with BAG shifts n^ to slighly higher values and shrinks its 
uncertainty by about 30 % when only high-/’ temperature is con¬ 
sidered, and by only about 15% when high-/’ temperature and 
polarization are combined. In slow-roll inflation, the running of 
the scalar spectral index is connected to the third derivative of 
the potential (Kosowsky & Turner, 1995). As was the case for 
the nominal mission results, values of the running compatible 
with the Planck 2015 constraints can be obtained in viable infla¬ 
tionary models (Kobayashi & Takahashi, 201 1). 
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Fig. 3. Comparison of the marginalized joint 68 % and 95 % CL constraints on («s, t) (left panel), (n ^, Q.bh^) (middle panel), and 
(Hs ,0mc) (right panel), for Planck 2013 (grey contours), Planck TT+lowP (red contours), Planck TT+lowP+lensing (green con¬ 
tours), and Planck TT,TE,EEh-1owP (blue contours). 



0.945 0.960 0.975 



0.945 0.960 0.975 



0.945 0.960 0.975 


-1-1-r 

I-1 Planck 20r3 

^ * P/anc/( TT+lowP 

t I P/anc/( TT+lowP+lensing 

I-1 Planck TT,TE,EE+lowP 



0.92 0.94 0.96 0.98 1.00 


Fig. 4. Marginalized joint 68 % and 95 % CL for (n^ , dus/d In k) 
using Planck TT-rlowP and Planck TT,TE,EEh-1owP. For com¬ 
parison, the thin black stripe shows the prediction for single-held 
monomial chaotic inhationary models with 50 < Nt < 60. 


ning, which would require new dedicated A^-body simulations. 
Hence we do not include Lycr constraints here. 

The fact that the data cannot accommodate a signihcant run¬ 
ning, but are compatible with a larger running of the running, 
will also be discussed in the inhaton potential reconstructions of 
Sect. 7. 


4.4. Suppression of power on the largest scales 

Although not statistically signihcant, the trend for a negative 
running or positive running of the running observed in the last 
subsection was driven by the lack of power in the Planck tem¬ 
perature power spectrum at low multipoles, already mentioned 
in the Planck 2013 release. This dehcit could potentially be ex¬ 
plained by a primordial spectrum featuring a depletion of power 
only at large wavelengths. Here we investigate two examples of 
such models. 

We hrst update the analysis (already presented in PCI 13) of 
a power-law spectrum multiplied by an exponential cut-off: 


Vnik) = mk) 


1 1 - exp 



( 20 ) 


When the running of the running is allowed to boat, the 
Planck TTh-IowP {Planck TT,TE,EEh-1owP) data give: 

Ms = 0.9569 + 0.0077 (0.9586 + 0.0056), 

dnjdlnk = 0.01l!fJ;[||^ (0.009 ± 0.010), (68 % CL) (19) 

d\jdlnP = 0.029!JJ;J^|^ (0.025 + 0.013), 

at the pivot scale = 0.05 Mpc '. Allowing for running of 
the running provides a better ht to the temperature spectrum 
at low multipoles, such that Ax^ ~ -4.8 (-4.9) for TT-tlowP 
(TT,TE,EEh-1owP). 

Note that the inclusion of small-scale data such as Lya might 
further constrain the running of the spectral index and its deriva¬ 
tive. The recent analysis of the BOSS one-dimensional Lya hux 
power spectrum presented in Palanque-Delabrouille et al. (2014) 
and Rossi et al. (2014) was optimised for measuring the neutrino 
mass. It does not include constraints on the spectral index run¬ 


This simple parameterization is motivated by models with a 
short inhationary stage in which the onset of the slow-roll phase 
coincides with the time when the largest observable scales ex¬ 
ited the Hubble radius during inhation. The curvature spectrum 
is then strongly suppressed on those scales. We apply top-hat 
priors on the parameter Tc, controlling the steepness of the cut¬ 
off, and on the logarithm of the cut-off scale, k^. We choose 
prior ranges Tc e [0,10] and ln(kc/Mpc“') e [-12,-3]. For 
Planck TTh-IowP (Planck TT,TE,EEh-1owP), the best-ht model 
has Tc = 0.50 (0.53), \n{kJMpc-^) = -7.98 (-7.98), n, = 
0.9647 (0.9649), and improves the effective x^ by a modest 
amount, Ax^ ~ -3.4 (-3.4). 

As a second model, we consider a broken-power-law spec¬ 
trum for curvature perturbations: 




( , \Ws-l+5 

if k < kb, 
ifk>kb. 


(21) 
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with Alow = A^ikblKY^ to ensure continuity at A: = kb- Hence 
this model, like the previous one, has two parameters, and also 
suppresses power at large wavelengths when 6 > 0. We assume 
top-hat priors on d 6 [0,2] and ln(A:b/Mpc“*) e [-12,-3], and 
standard uniform priors for ln(10'°As) and n^. The best ht to 
Planck TTh-IowP {Planck TT,TE,EEh-1owP) is found for «s = 
0.9658 (0.9647), 6 = 1.14 (1.14), and ln(kb/Mpc^') = -7.55 
(-7.57), with a very small improvement of ~ -1.9 
(- 1 . 6 ). 

We conclude that neither of these two models with two extra 
parameters is preferred over the base ACDM model. (See also 
the discussion of a step inflationary potential in Sect. 9.1.1.) 


5. Constraints on tensor modes 

In this section, we focus on the Planck 2015 constraints on tensor 
perturbations. Unless otherwise stated, we consider that the ten¬ 
sor spectral index satisfies the standard inflationary consistency 
condition to lowest order in slow-roll, nt = -r/8. We recall that r 
is defined at the pivot scale k* = 0.05 Mpc *. However, for com¬ 
parison with other studies, we also report our bounds in terms of 
the tensor-to-scalar ratio ro.oo 2 at K - 0.002 Mpc“'. 


5 .1. Planck 2015 upper bound on r 

The constraints on the tensor-to-scalar ratio infered from the 
Planck full mission data for the ACDM-rr model are: 


^0.002 < 0.10 
^0.002 < 0.11 
^0.002 < 0.11 
ro.oo 2 < 0.10 


(95 % CL, Planck TT-rlowP), 

(95 % CL, Planck TT-rlowP-rlensing), 
(95 % CL, Planck TTh-IowPh-BAO) , 
(95 % CL, Planck TT,TE,EEh-1owP) . 


( 22 ) 

(23) 

(24) 

(25) 


Table 4 also shows the bounds on n^ in each of these cases. 

These results slightly improve over the constraint ro.oo 2 < 
0.12 (95% CL) derived from the Planck 2013 temperature 
data in combination with WMAP large-scale polarization data 
(Planck Collaboration XVI, 2014; Planck Collaboration XXII, 
2014). The constraint obtained by Planck temperature and po¬ 
larization on large scales is tighter than Planck Z?-mode 95 % CL 
upper limit from the 100 and 143 GHz HEI channels, r < 0.27 
(Planck Collaboration XI, 2015). The constraints on r reported in 
Table 4 can be translated into upper bounds on the energy scale 
of inflation at the time when the pivot scale exits the Hubble 
radius, using 




37r^As 

2 


rMl, = (1.88 X 10'® GeV)"'— . 
P' 0.10 


(26) 


This gives an upper bound on the Hubble parameter during in¬ 
flation of HtlMpi < 3.6 X 10 ® (95 % CL) for Planck TT-rlowP 
These bounds are relaxed when allowing for a scale depen¬ 
dence of the scalar and tensor spectral indices. In that case, we 
assume that the tensor spectral index and its running are fixed by 
the standard inflationary consistency condition at second order 
in slow-roll. We obtain 


ro.oo 2 <0.18 (95 % CL, Planck TT-rlowP), (27) 

dfi 

-0.013+°;®^° (68 % CL, P/anck TTh-IowP), (28) 

with = 0.9667 + 0.0066 (68 % CL). At the standard pivot 
scale, k* = 0.05 Mpc^', the bound is stronger (r < 0.17 at 95 % 
CL), because k* is closer to the scale at which Ms and r decorre¬ 
late. The constraint on ro.oo 2 in Eq. (27) is 21 % tighter compared 



0.92 0.94 0.96 0.98 1.00 


Fig. 5. Marginalized joint confidence contours for 
(Hs ,dns/dlnk), at the 68% and 95% CL, in the presence 
of a non-zero tensor contribution, and using Planck TT-rlowP or 
Planck TT,TE,EEh-1owP. Constraints from the Planck 2013 data 
release are also shown for comparison. The thin black stripe 
shows the prediction of single-field monomial inflation models 
with 50 < A« < 60. 

with the Planck 2013 results. The mean value of the running in 
Eq. (28) is higher (lower in absolute value) than with Planck 
2013 by 45 %. Figures 5 and 6 clearly illustrate this significant 
improvement with respect to the previous Planck data release. 
Table 4 shows how bounds on (r, n^, dn^ldlnk) are affected by 
the lensing reconshTiction, BAG, or high-f polarization data. The 
tightest bounds are obtained in combination with polarization: 

^0.002 < 0.15 

(95 % CL, Planck TT,TE,EEh-1owP), (29) 

= -0.009 + 0.008 
Ink 

(68 % CL, Planck TT,TE,EE+lowP), (30) 

with ws = 0.9644 + 0.0049 (68 % CL). 

Neither the Planck full mission constraints in Eqs. (22)-(25), 
nor those including a running in Eqs. (27) and (29), are com¬ 
patible with the interpretation of the BICEP2 B-mode polariza¬ 
tion data in terms of primordial gravitational waves (BICEP2 
Collaboration, 2014b). Instead, they are in excellent agree¬ 
ment with the results of the BICEP2/Keck Anay-Planck cross¬ 
correlation analysis, as discussed in Sect. 13. 

5.2. Dependence of the r constraints on the low-( 
likelihood 

The constraints on r discussed above are further tightened by 
adding WMAP polarization information on large angular scales. 
The Planck measurement of CMB polarization on large angular 
scales at 70 GHz is consistent with the WMAP 9-year one, based 
on the K, Q, and V bands (at 30, 40, and 60 GHz, respectively), 
once the Planck 353 GHz channel is used to remove the dust 
contamination, instead of the theoretical dust model used by the 
WMAP team (Page et ah, 2007). (For a detailed dicussion, see 
Planck Collaboration XI (2015).) By combining Planck TT data 
with LEI 70 GHz and WMAP polarization data on large angular 
scales, we obtain a 35 % reduction of uncertainty, giving r = 
0.074 ±0.012 (68 % CL) and n^ = 0.9660 ± 0.060 (68 % CL) 
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Model 

Parameter 

Planck TT+lowP 

Planck TT+lowP+lensing 

Planck TT+lowP+BAO 

Planck TT,TE,EE+lowP 

ACDM+r 

^0.002 

-2Aln£,^,x 

0.9666 ± 0.0062 

< 0.103 

0 

0.9688 + 0.0061 

< 0.114 

0 

0.9680 ± 0.0045 

< 0.113 

0 

0.9652 + 0.0047 

< 0.099 

0 

ACDM+r 

+d«s/dlnk 

ris 

^0.002 

r 

dfis/dlnk 

-2Alni;max 

0.9667 ± 0.0066 

< 0.180 
< 0.168 

-0.81 

0.9690 ± 0.0063 

<0.186 

<0.176 

-0 0076+° °°®^ 

-0.08 

0.9673 ± 0.0043 

< 0.176 

< 0.166 

-0.0125 + 0.0091 

-0.87 

0.9644 + 0.0049 

<0.152 
< 0.149 

-0.0085 + 0.0076 

-0.38 


Table 4. Constraints on the primordial perturbation parameters for ACDM+r and ACDM+r+dns/dlnA: models from Planck. 
Constraints on the spectral index and its dependence on the wavelength are given at the pivot scale of - 0.05 Mpc '. 



0.94 0.96 0.98 1.00 

Primordial tilt (rig) 


Fig. 6. Marginalized joint confidence contours for {n ^, r), at the 
68 % and 95 % CL, in the presence of running of the spectral 
indices, and for the same data combinations as in the previous 
figure. 

for the base ACDM model. When tensors are added, the bounds 
become 

ro.oo 2 < 0.09 (95 % CL, Planck TT+lowP+WP),(31) 

ns = 0.9655 + 0.058 (68 % CL, Planck TT+lowP+WP),(32) 
T = 0.073:|;° °J^ (68 % CL, Planck TT+lowP+WP).(33) 

When tensors and running are both varied, we obtain ro.oo 2 < 
0.14 (95% CL) and dnjdlnk = -0.010 + 0.008 (68% CL) 
for Planck TT+lowP+WP. These constraints are all tighter than 
those based on Planck TT+lowP only. 

5.3. The tensor-to-scalar ratio and the tow-€ deficit in 
temperature 

As noted previously (Planck Collaboration XV, 2014; Planck 
Collaboration XVI, 2014; Planck Collaboration XXII, 2014), the 
low-f’ temperature data display a slight lack of power compared 
to the expectation of the best-fit tensor-free base ACDM model. 
Since tensor fluctuations add power on small scales, the effect 
will be exacerbated in models allowing r > 0. 

In order to quantify this tension, we compare the observed 
constraint on r to that inferred from simulated Planck data. In 
the simulations, we assume the underlying fiducial model to be 


tensor-free, with parameters close to the base ACDM best-fit val¬ 
ues. We limit the simulations to mock temperature power spec¬ 
tra only and fit these spectra with an exact low-f likelihood for 
2 <(< 29 (see Perotto et al. (2006)), and a high-^ Gaussian 
likelihood for 30 < f < 2508 based on the frequency-combined, 
foreground-marginalized, unbinned Planck temperature power 
spectrum covariance matrix. Additionally, we impose a Gaussian 
prior of T = 0.07 + 0.02. 

Based on 100 simulated data sets, we find an expectation 
value for the 95 % CL upper limit on the tensor-to-scalar ratio of 
fia- ~ 0.260. The corresponding constraint from real data (using 
low-^ Commander temperature data, the frequency-combined, 
foreground-marginalized, unbinned Planck high-f TT power 
spectrum, and the same prior on r as above) reads r < 0.123, 
confirming that the actual constraint is tighter than what one 
would have expected. However, the actual constraint is not ex¬ 
cessively unusual: out of the 100 simulations, 4 lead to an even 
tighter bound, corresponding to a significance of about 2cr. 
Thus, under the hypothesis of the base ACDM cosmology, the 
upper limit on r that we get from the data is not implausible as a 
chance fluctuation of the low multipole power. 

To illustrate the contribution of the low-^ temperature power 
deficit to the estimates of cosmological parameters, we show as 
an example in Fig. 7 how n^ shifts towards lower values when 
the { < 30 temperature information is discarded (we will refer to 
this case as “Planck TT-lowT”). The shift in n^ is approximately 
-0.005 (or -0.003 when the lowP likelihood is replaced by a 
Gaussian prior t = 0.07 + 0.02). Comparing with Sect. 4.4, we 
notice that these shifts are larger than when fitting a primordial 
spectrum suppressed on large scales to the data. 

Figure 8 displays the posterior probability for r for various 
combinations of data sets, some of which exclude the { <30 TT 
data. This leads to the very conservative bounds r < 0.24 and 
r < 0.23 at 95 % CL when combined with the lowP likelihood 
or with the Gaussian prior t = 0.07 + 0.020, respectively. 

5.4. Relaxing assumptions on the late-time cosmological 
evolution 

As in the Planck 2013 release (PCI 13), we now ask how robust 
the Planck results on the tensor-to-scalar ratio are against as¬ 
sumptions on the late-time cosmological evolution. The results 
are summarized in Table 5, and some particular cases are illus¬ 
trated in Fig. 9. Constraints on r turn out to be remarkably stable 
for one-parameter extensions of the ACDM-i-r model, with the 
only exception the ACDM-i-r-i-Q/j' case in the absence of the late¬ 
time information from Planck lensing or BAG data. The weak 
trend towards Q/f < 0 , i.e. towards a positively curved (closed) 
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Extended model. 

Parameter 

Planck TT-l-lowP 

Planck TT±lowP 

Planck TT,TE,EE 

ACDM-l-r-l- 


-l-lensing 

±BAO 

±lowP 

-l-general reionization 

r 

ns 

r< 0.11 

0.975 ± 0.006 

r < 0.10 

0.971 ± 0.005 

< 0.10 

0.968 ± 0.005 


r 

< 0.14 

< 0.12 

< 0.11 

+Aeff 

ns 

n Q77+0.016 
' -0.017 

0.972 ± 0.009 

0.964 ±0.010 


Aeff 

o 74 - 1 - 0.30 
-0 35 

3.19 ±0.24 



r 

< 0.14 

< 0.12 

< 0.12 

+ THe 

lls 

0.975 ± 0.007 

0.973 ± 0.009 

0.969 ± 0.008 


Tue 

0.258 ± 0.022 

0.257 ± 0.022 

0.252 ±0.014 


r 

< 0.11 

< 0.11 

< 0.11 

+my 

ns 

0.963 ± 0.007 

0.967 ± 0.005 

0.962 ± 0.005 


niy 

< 0.67 

< 0.21 

<0.58 


r 

< 0.15 

r < 0.11 

< 0.15 

+Q.g 

ns 

0.971 ±0.007 

0.971 ± 0.007 

0.969 ± 0.005 


Qk 

-o.oo 8 !;;!!>« 

-0.001 ± 0.003 

-0 

'-'•'-'^^-0 n?n 


r 

< 0.14 

< 0.11 

< 0.12 

+w 

ns 

0.969 ± 0.006 

0.967 ± 0.006 

0.966 ± 0.005 


w 

_1 4 ^+ 0.20 

A •‘+^-0.40 


-1 


r 

< 0.20 

< 0.18 

<0.19 

+Q.K+dns/dlnk 

ns 

0.971 ±0.007 

0.969 ± 0.007 

0.969 ± 0.005 

dn^jdlnk 

-0.006 ± 0.009 

-0.013 ±0.009 

-0.004 ± 0.008 


fl/C 

u.uuu^ noQ 

-0.001 ± 0.003 

-0 


r 

r < 0.14 

r<0.13 

< 0.12 

+Aeff+meff 

ih 

Weff 

0 980+° °'“ 

< 0.27 

00,_1 

00 0 

V 

d 

0 968+° °°° 

U. 5 'UO_gQ(,g 

<0.83 


Aeff 

< 3.45 

< 3.73 

< 3.47 


Table 5. Constraints on extensions of the ACDM+r cosmological model for Planck TT+lowP+lensing, Planck TT+lowP+BAO, and 
Planck TT,TE,EE+lowP Eor each model we quote 68 % CL, unless 95 % CL upper bounds are reported. 




Fig. 7. One dimensional posterior probabilities for Wj for the base 
ACDM model obtained by excluding temperature multipoles for 
f < 30 (“TT-lowT”), while either keeping low-^ polarization 
data, or in addition replacing them with a Gaussian prior on r. 


universe from the temperature and polarization data alone, and 
the well-known degeneracy between D.k and HolPl^a lead to a 
slight suppression of the Sachs-Wolfe plateau in the scalar tem¬ 
perature spectrum. This leaves more room for a tensor compo¬ 
nent. 


Fig. 8. One dimensional posterior probabilities for r for various 
data combinations, either including or not including temperature 
multipoles for { < 30, and compared with the baseline choice 
(Planck TTh-IowP, black curve). 


This further degeneracy when r is added builds on the neg¬ 
ative values for the curvature allowed by Planck TT-rlowP, 
PIk = -0.052+jJ “g at 95 % CL (Planck Collaboration XIII, 
2015). The exploitation of the information contained in 
the Planck lensing likelihood leads to a tighter constraint, 
fl/c - -0.005)^q q}® at 95 % CL, which improves on the Planck 
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Fig. 9. Marginalized joint 68 % and 95 % CL for (n ^, ro. 002 ) 
using Planck TT+lowP+BAO (upper panel) and Planck 
TT,TE,EE+lowP (lower panel). 


6. Implications for single-field slow-roll inflation 

In this section we study the implications of Planck 2015 con¬ 
straints on standard slow-roll single-field inflationary models. 

6.1. Constraints on slow-roll parameters 

We first present the Planck 2015 constraints on slow-roll pa¬ 
rameters obtained through the analytic perturbative expansion 
in terms of the Hubble flow functions (HEEs) e, for the primor¬ 
dial spectra of cosmological fluctuations during slow-roll infla¬ 
tion (Stewart & Lyth, 1993; Gong & Stewart, 2001; Leach et al., 
2002). When restricting to first order in e,-, we obtain 

ei < 0.0068 (95 % CL, Planck TT-tlowP), (34) 

62 = 0.029!°;“ (68 % CL, Planck TT-tlowP). (35) 

When high-^ polarization is included we obtain e\ < 0.0066 at 
95 % CL and 62 = 0.030!° °°’ at 68 % CL. When second-order 
contributions in the HEEs are included, we obtain 

61 < 0.012 (95 % CL, Planck TTh-IowP) , (36) 

62 = 0.031!° °|° (68 % CL, Planc/fc TTh-IowP) , (37) 

-0.41 < 63 < 1.38 (95 % CL, Planck TTh-IowP) . (38) 

When high-f polarization is included we obtain 61 < 0.011 at 
95 % CL, 62 = 0.032!° °^^ at 68 % CL, and -0.32 < 63 < 0.89 at 
95 % CL. 


I I Planck 2013 CS3 Planck TT-i-lowP i I Planck TT,TE,EE+lowP 



o -0.04 -0.02 0.00 0.02 


Vv 


2013 results (Q.k = -0.007!° °[* at 95% CL). However, due to 
the remaining degeneracies left by the uncertainties in polariza¬ 
tion on large angular scales, a full appreciation of the improve¬ 
ment due to the full mission temperature and lensing data can be 
obtained by using 1 owPh-WP, which leads to = -0.003!g ^[4 
at 95 % CL. Let us also Anally note that negative values allowed 
for the curvature are decreased when the running is combined, 
showing that the low-^ temperature deficit could also contribute 
to the estimate of these parameters. 

The trend found for ACDM-nrH-G/j- is even clearer when spa¬ 
tial curvature and the running of the spectral index are varied at 
the same time. Then, the Planck temperature plus polarization 
data are compatible with r values as large as 0.19 (95 % CL), at 
the cost of an almost 4 cr deviation from spatial flatness (which, 
however, disappears as soon as lensing or BAG data are consid¬ 
ered). 


Fig. 10. Marginalized joint 68 % and 95 % CL regions for ( 61 , 62 ) 
(top panel) and (ey , rjv) (bottom panel) for Planck TTh-IowP (red 
contours), Planck TT,TE,EEh-1owP (blue contours), and com¬ 
pared with the Planck 2013 results (grey contours). 


The potential slow-roll parameters are obtained as derived 
parameters by using their exact expressions as function of 6 , 
(Leach et al., 2002; Finelli et al., 2010): 


ev = 


Vv 


, (1 - T + I) 

2y2 


(1 -1) 


_ *^2 _ ^=1 , 3 E 1 E 2 _ '2 _ CM 

^G-y-i- + “6“ T2 “6“ 


26 ? 


5£iE2 


1 - ^ 
^ 3 


(39) 
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(40) 




















Planck Collaboration: Constraints on inflation 


15 








y2 


4ef - 3eie2 + 


3efe2 


1 _ IL j_ ^ Z' ^ ^ 

^ 3 6 I < 9 ^ ^9 ^3 

7 ^2^3 ^2^2 £2 ^3 £4 

._,,e,e3 + —+ —+ ^ 


eiei 


(41) 


where V((p) is the inflaton potential, the subscript (p denotes the 
derivative with respect to p, and Mpi = ( 87 rG)“'^^ is the reduced 
Planck mass (see also Table 2). 

By using Eqs. (39) and (40) with 63 = 64 = 0 and the pri¬ 
mordial power spectra to lowest order in the HFFs, the derived 
constraints for the first two slow-roll potential parameters are: 


ey < 0.0068 (95 % CF, Planck TT-blowP), (42) 

77y= -0.010+°;°“ (68 % CF, Ptonck TTh-IowP) . (43) 

When high-^ polarization is included we obtain ey < 0.0067 at 
95 % CF and rjv = -0.010+° °“ at 68 % CF. By using Fqs. (39), 
(40), and (41) with 64 = 0 and the primordial power spectra to 
second order in the FIFFs, the derived constraints for the slow- 
roll potential parameters are: 

ey< 0.012 (95 % CF, P/anck TTh-IowP) , (44) 

77 y= -0.0080!°;°“^ (68 % CF, P/anck TTh-IowP) , (45) 

^y= 0 . 0070 +° °“5 (68 % CF, P/anck TTh-IowP) . (46) 

When high-Z” polarization is included we obtain ey < 0.011 at 
95 % CF, and t]v = -0.0092+° °°“ ^2 ^ o.0044+° °“ 7 ^ 

at 68 % CF. 

In Figs. 10 and 11 we show the 68 % CF and 95 % CF of the 
HFFs and the derived potential slow-roll parameters with and 
without the inclusion of high-f polarization, comparing with the 
Planck 2013 results. 


6.2. Implications for selected inflationary models 

The predictions to lowest order in the slow-roll approximation 
for (Hs, r) of a few inflationary models with a representative un¬ 
certainty for the entropy generation stage (50 < Nt < 60) are 
shown in Fig. 12. 

In the following we discuss the implications of Planck 
TTh-IowPh-BAO data for selected slow-roll inflationary models 
by taking into account the uncertainties in the entropy genera¬ 
tion stage. We model these uncertainties by two parameters, as 
in PCI13: the energy scale pth by which the Universe has ther- 
malized, and the parameter wint which characterizes the effec¬ 
tive equation of state between the end of inflation and the energy 
scale specified by pth. For each inflationary model we provide 
in Table 6 and in the main text the value with respect to 
the base ACDM model and the Bayesian evidence with respect 
to the inflationary model (Starobinsky, 1980), computed by 
CosmoMC connected to CAMB using MultiNest as the sam¬ 
pler. We use the primordial power spectra of cosmological fluc¬ 
tuations generated during slow-roll inflation parameterized by 
the HFFs, e,, to second order, which can be expressed in terms 
of the number of e-folds to the end of inflation, A*, and the pa¬ 
rameters of the considered inflationary model, using modified 
routines of the public code ASPIC® (Martin et al., 2014). For 
the number of e-folds to the end of inflation (Fiddle & Feach, 


® http : //cp3 . irmp . ucl . ac . be/~ringeval/aspic . 
html 


2003; Martin & Ringeval, 2010) we use the expression (PCI13) 
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where pend is the energy density at the end of inflation, aoTfo is 
the present Hubble scale, V* is the potential energy when k, left 
the Hubble radius during inflation, and g* is the number of ef¬ 
fective bosonic degrees of freedom at the energy scale pth. We 
consider the pivot scale k» = 0 . 002 Mpc^', g* = 10 ^, eend = 1 , 
and a logarithmic prior on pth (in the interval [( 10 ° GeV)"^,pend])- 
We have validated the slow-roll approach by cross-checking 
the Bayes factor computations against the fully numerical in¬ 
flationary mode equation solver ModeCode coupled to the 
PolyChord sampler. 


Power-law potentials 

We first investigate the class of inflationary models with a single 
monomial potential (Finde, 1983): 

in which inflation occurs for large values of the inflaton (p > Mpp 
The predictions for the scalar spectral index and the tensor- 
to-scalar ratio at first order in the slow-roll approximation are 
«s - 1 ~ -2(n + 2)/(4At -I- n) and r x 16n/(4At + n), respec¬ 
tively. By assuming a dust equation of state (i.e., Wint = 0) prior 
to thermalization, the cubic and quartic potentials are strongly 
disfavoured by InB = -11.6 andlnB = -23.3, respectively. The 
quadratic potential is moderately disfavoured by In B = -4.7. 
Other values, such as n = 4/3, 1, and 2/3, motivated by ax- 
ion monodromy (Silverstein & Westphal, 2008; McAllister et al., 
2010), are compatible with Planck data with Wint = 0. 

Small modifications occur when considering the effective 
equation of state parameter, Wint = (n - 2) f(n + 2 ), defined by 
averaging over the coherent oscillation regime which follows in¬ 
flation (Turner, 1983). The Bayes factors are slightly modified 
when Wint is allowed to float, as can be seen from Table 6 . 


Hilltop models 

In hilltop models (Boubekeur & Fyth, 2005), with potential 

V(,^)^A4|l-^ + ...j, (49) 

the inflaton rolls away from an unstable equilibrium. The 
predictions to first order in the slow-roll approximation are 
r ^ 8 p°(Mpi///) 2 x 2 '’^ 2 /(l - xPf and n, - 1 ^ -2p{p - 

l)(Mpi/ju)2x^-V(l - xP) - 3r/8, where x = (p,lp. As in PCI13, 
the ellipsis in Fq. (49) and in what follows indicates higher-order 
terms that are negligible during inflation but ensure positiveness 
of the potential. 

By sampling logjo(/ 7 /Mpi) within the prior [0.30,4.85] for 
p - 2, we obtain log[Q(ju/Mpi) > 1.02 (1.05) at 95% CF and 
InB = -2.6 (-2.4) for Wint = 0 (allowing Wint to float). 

An exact potential which could also apply after inflation, 
instead of the approximated one in Eq. (49), might be needed 
for a better comparison among different models. Hilltop mod¬ 
els in Eq. (49) approximate a linear potential V{^) oc (p for 
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Fig. 11. Marginalized joint 68% and 95% CL regions for (ei , 62 ,^ 3 ) (top panels) and {ey ,t]v (bottom panels) for Planck 
TT+lowP (red contours), Planck TT,TE,EE+lowP (blue contours), and compared with the Planck 2013 results (grey contours). 
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filMpi » 1. By considering a double-well potential, V((p) — 
A^[l -0^/(2/i^)]^, instead, we obtain a slightly worse Bayes fac¬ 
tor than the hilltop p - 2 model, InZ? = -3.1 (-2.3) for Wint = 0 
(wint allowed to vary). This different result can be easily under¬ 
stood. Although the double-well potential is equal to the hilltop 
model for (p <s: p, it approximates V{<p) oc for pIMpi » 1. 
Since a linear potential is a better ht to Planck than 0^, the ht 
of the double-well potential is therefore worse than the hilltop 
p — 2 case for pjMpx » 1, and this partially explains the slightly 
different Bayes factors obtained. 

In the p - A case, we obtain log[g(ju/Mpi) > 1.05 (1.02) at 
95% CL and InB = -2.8 (-2.6) for Wint = 0 (allowing Wint to 
float), assuming a prior range [-2,2] for logjQ(/r/Mpi). 

Natural inflation 

In natural inflation (Freese et ah, 1990; Adams et ah, 1993) a 
non-perturbative shift symmetry is invoked to suppress radiative 
corrections, leading to the periodic potential 


As for the hilltop models described earlier, the ellipsis indicates 
possible higher-order terms that are negligible during inflation 
but ensure positiveness of the potential. These models predict 
r ss and n^ - \ ~ -q^e^‘i'>’/^<‘'(2 + 

e^q‘f’lMpiy(^l - jjj ^ slow-roll trajectory characterized 

by ~ fi(f>IMpi) - /(0end/Afpi), with f{x) = (e«-‘‘ - qx)lq^. 
By sampling log[o(q'/Mpi) within the prior [-3,3], we obtain a 
Bayes factor of -0.6 for Wint = 0 (-0.9 when Wint is allowed to 
vary). 

Spontaneously broken SUSY 

Hybrid models (Copeland et ah, 1994; Linde, 1994) predict¬ 
ing Ws > 1 are strongly disfavoured by the Planck data, as for 
the first cosmological release (PCI13). An example of a hybrid 
model predicting ns < 1 is the case in which slow-roll infla¬ 
tion is driven by loop corrections in spontaneously broken su¬ 
persymmetric (SUSY) grand unified theories (Dvali et ah, 1994) 
described by the potential 


y(0) = A^ 


1 + cos 



(50) 


y(0) = A4[l+a„log(0/Mpi)], (53) 


where / is the scale which determines the curvature of the po¬ 
tential. We sample log[g(//Mpi) within the prior [0.3,2.5] as in 
PCI13. We obtain logio(//Mpi) > 0.84 (> 0.83) at 95 % CL and 
InB = -2.4 (-2.3) for Wint = 0 (allowing Wint to vary). 

Note that the super-Planckian value for / required by ob¬ 
servations is not necessarily a problem for this class of models. 
When several fields with a cosine potential as in Eq. (50) and 
scales fi appear in the Lagrangian, an effective single field infla¬ 
tionary trajectory can be found for a suitable choice of param¬ 
eters (Kim et ah, 2005). In such a setting, the super-Planckian 
value of the effective scale / required by observations can be 
obtained even if the original scales statisfy /]■ « Mp\ (Kim et ah, 
2005). 

D-brane inflation 

Inflation can be caused by physics in extra dimensions. If 
the standard model of particle physics is confined to our 3- 
dimensional brane, the distance between our brane and anti- 
brane can drive inflation. We consider the following parameteri¬ 
zation for the effective potential driving inflation: 

V(^) = a4(i_L+ ...). (51) 

We sample logjglju/Mpi) within the prior [-6,0.3]. We consider 
p - A (Kachru et ah, 2003; Dvali et ah, 2001) and p - 2 (Garcia- 
Bellido et ah, 2002). The predictions for r and n^ can be obtained 
from the hilltop case with the substitution p —> -p. These mod¬ 
els agree with the Planck data with a Bayes factor of In B = -0.4 
(-0.6) and In B = -0.7 (-0.9) for p - A and p — 2, respectively, 
for Wint = 0 (allowing wint to vary). 

Exponential potentials 

Exponential potentials are ubiquitous in inflationary models mo¬ 
tivated by supergravity and string theory (Goncharov & Linde, 
1984; Stewart, 1995; Dvali & Tye, 1999; Burgess et ah, 2002; 
Cicoli et ah, 2009). We restrict ourselves to the analysis of the 
following class of potentials: 

y(0) = A"^ (l - H-...) . (52) 


where ah > 0 is a dimensionless parameter. Note that for 
a/, « 1, this model leads to the same predictions as the power- 
law potential for p « 1 to lowest order in the slow-roll approxi¬ 
mation. By sampling log[g(a/,) on a flat prior [-2.5,1], we obtain 
a Bayes factor of -2.2 for Wint = 0 (-1.7 when Wint is allowed to 
vary). 

inflation 

The first inflationary model proposed (Starobinsky, 1980), with 
action 

(54) 

is still within the Planck 68 % CL constraints, as it was for the 
Planck 2013 release (PCI13). This model corresponds to the po¬ 
tential 

V{f) = A'^ (l - g-V2/30/Mpi^^ (55) 

in the Einstein frame, which leads to the slow-roll predictions 
Ms - 1 ~ -2/A and r « 12/A^ (Starobinsky, 1980; Mukhanov & 
Chibisov, 1981). 

After the Planck 2013 release, several theoretical develop¬ 
ments supported the model in Eq. (54) beyond the original mo¬ 
tivation of including quantum effects at one-loop (Starobinsky, 
1980). No-scale supergravity (Ellis et ah, 2013a), the large field 
regime of superconformal D-term inflation (Buchmiiller et ah, 
2013), or recent developments in minimal supergravity (Earakos 
et ah, 2013; Eerrara et ah, 2013b) can lead to a generalization 
of the potential in Eq. (55) (see Ketov & Starobinsky (201 1) for 
a previous embedding of inflation in FiR.) supergravity). The 
potential in Eq. (55) can also be generated by spontaneous break¬ 
ing of conformal symmetry (Kallosh & Linde, 2013). This infla¬ 
tionary model has A;^f^ 0.8 (0.3) larger than the base ACDM 

model with no tensors for Wint = 0 (for Wint allowed to vary). We 
obtain 54 < A, < 62 (53 < A* < 64) at 95 % CL for Wint = 0 
(for Wint allowed to vary), compatible with the theoretical pre¬ 
diction, At = 54 (Starobinsky, 1980; Vilenkin, 1985; Gorbunov 
& Panin, 2011). 
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a attractors 

We now study two classes of inflationary models motivated by 
recent developments in conformal symmetry and supergravity 
(Kallosh et al., 2013). The first class has been motivated by con¬ 
sidering a vector rather than a chiral multiplet for the inflaton in 
supergravity (Ferrara et ah, 2013a) and corresponds to the po¬ 
tential (Kallosh et ah, 2013) 

V((P) = (l - g-V2^/(V3SMp,)j^ _ 

To lowest order in the slow-roll approximation, these models 
predict r 64/[3q'( 1 - g V 20 /(V 3 ^Mp,)) 2 ] n, - 1 - 8(1 -H 

e'^'^/('^“pi))/[ 3 a(l - g 2 / 20 /(V 3 aMpi)^ 2 j Qjj inflationary trajec¬ 
tory characterized hy N ^ g(<^/Mpi) - g(0end/4^pi) with g(x) = 

(3a^e - Vto.r)/4. The relation between N and cp can be 

inverted through the use of the Lambert functions, as done for 
other potentials (Martin et al., 2014). By sampling log[g(Q:^) on 
a flat prior [0,4], we obtain log[o(a^) <1.7 (2.0) at 95% CL 
and a Bayes factor of -1.8 (-2) for Wint = 0 (for Wint allowed to 
vary). 

The second class of models has been called super-conformal 
a attractors (Kallosh et al., 2013) and can be seen as originat¬ 
ing from a different generating function with respect to the first 
class. This second class is described by the following potential 
(Kallosh etal., 2013): 


Vicp) = tanh^'" 



(57) 


Table 6. Results of the inflationary model comparison. We pro¬ 
vide with respect to base ACDM and Bayes factors with 
respect to inflation. 


Inflationary model 

o 

% 

a 

II 

In Sox 

Wint = 0 Wint t 0 

R + R-I( 6 M^) 

+0.8 

+0.3 


+0.7 

n = 2/3 

+6.5 

+3.5 

-2.4 

-2.3 

n = 1 

+6.2 

+5.5 

-2.1 

-1.9 

n = 4/3 

+6.4 

+5.5 

-2.6 

-2.4 

n = 2 

+8.6 

+8.1 

-A.l 

-4.6 

« = 3 

+22.8 

+21.7 

-11.6 

-11.4 

n = A 

+43.3 

+41.7 

-23.3 

-22.7 

Natural 

+7.2 

+6.5 

-2.4 

-2.3 

Hilltop (p = 2) 

+4.4 

+3.9 

-2.6 

-2.4 

Hilltop (p = 4) 

+3.7 

+3.3 

-2.8 

-2.6 

Double well 

+5.5 

+5.3 

-3.1 

-2.3 

Brane inflation (p = 2) 

+3.0 

+2.3 

-0.7 

-0.9 

Brane inflation (p = 4) 

+2.8 

+2.3 

-0.4 

-0.6 

Exponential inflation 

+0.8 

+0.3 

-0.7 

-0.9 

SB SUSY 

+0.7 

+0.4 

-2.2 

-1.7 

Supersymmetric a-model 

+0.7 

+0.1 

-1.8 

-2.0 

Superconformal (m = 1) 

+0.9 

+0.8 

-2.3 

-2.2 

Superconformal {m 1) 

+0.7 

+0.5 

-2.4 

-2.6 


Barvinsky et al. (2008) and Bezrukov & Shaposhnikov (2009) 
for the inclusion of loop corrections). It is therefore in agree¬ 
ment with the Planck constraints, as for the first cosmological 
data release (PCI 13). 


This is the simplest class of models with spontaneous breaking 
of conformal symmetry, and for a - m - \ reduces to the origi¬ 
nal model introduced by Kallosh & Linde (2013). The potential 
in Eq. (57) leads to the following slow-roll predictions (Kallosh 
et al.,2013): 


48am 

r ~ -, 

4mN^ + 2Ng(a, m) -i- 3am 

—SmN — 6am + 2Ng(a, m) 

_ 1 % _ 

AmN'^ + 2Ng(a, m) + 3am 

where g(a, m) - -\/3a(4m^ + 3a). The predictions of this second 
class of models interpolate between those of a large-field chaotic 
model, y(0) cx (fp-”', for a » 1 and the Rp model for a c 1. 

For a we adopt the same priors as the previous class in 
Eq. (56). By fixing m - 1, we obtain logjQ(a^) < 2.3 (2.5) at 
95 % CL and a Bayes factor of -2.3 (-2.2) for Wint = 0 (when 
Wint is allowed to vary). When m is allowed to vary as well with 
a flat prior in the range [0,2], we obtain 0.02 < m < 1 (m < 1) 
at 95 % CL for Wint = 0 (when Wint is allowed to vary). 

Non-minimally coupled inflaton 

Inflationary predictions are quite sensitive to a non-minimal cou¬ 
pling, ^R4P, of the inflaton to the Ricci scalar. One of the most 
interesting effects due to ^ 4= 0 is to reconcile the quartic poten¬ 
tial V{(p) - A(p‘^ jA with Planck observations, even for ^ c 1. 

The Higgs inflation model (Bezrukov & Shaposhnikov, 
2008), in which inflation occurs with V{(p) - AicfP - (p^PjA and 
^ » 1 for (p » (pQ, leads to the same predictions as the R^ model 
to lowest order in the slow-roll approximation at tree level (see 


(58) 

(59) 


Let us summarize our findings for Planck IowPh-BAO: 

- Monomial potentials with integer n > 2 are strongly dis¬ 
favoured with respect to R^. 

- The Bayes factor prefers R^ over chaotic inflation with 
monomial quadratic potential by odds of 110:1 under the 
assumption of a dust equation of state during the entropy 
generation stage. 

- R^ inflation has the strongest evidence among the models 
considered here. However, care must be taken not to overin¬ 
terpret small differences in likelihood lacking statistical sig¬ 
nificance. 

- The models closest to R^ in terms of evidence are brane infla¬ 
tion and exponential inflation, which have one more param¬ 
eter than R^. Both brane inflation considered in Eq. (51) and 
exponential inflation in Eq. (52) approximate the linear po¬ 
tential for a large portion of parameter space (for ju/Mpi » 1 
and q 'S> 1, respectively). For this reason these models have 
a higher evidence (although not at a statistically significant 
level) compared to those approximate a quadratic potential, 
as do a-attractors, for instance. 

- In the models considered here, the Ax^ obtained by allowing 
w to vary is modest, i.e., less than approximately 1.6 (with re¬ 
spect to Wint = 0). The gain in the logarithm of the Bayesian 
evidence is even smaller, since an extra parameter is added. 


7. Reconstruction of the potential and analysis 
beyond slow-roll approximation 

In the previous section, we derived constraints on several types 
of inflationary potentials, assumed to account for the inflaton dy¬ 
namics between the time at which the largest observable scales 
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crossed the Hubble radius during inflation and the end of infla¬ 
tion. The full shape of the potential was used in order to identify 
when inflation ends, and thus the held value <pt when the pivot 
scale crosses the Hubble radius. 

In section 6 of PCI13, we explored another approach, con¬ 
sisting of reconstructing the inflationary potential within its ob¬ 
servable range without making any assumptions concerning the 
inflationary dynamics outside that range. Indeed, given that the 
number of e-folds between the observable range and the end of 
inflation can always be adjusted to take a realistic value, any po¬ 
tential shape giving a primordial spectrum of scalar and tensor 
perturbations in agreement with observations is a valid candi¬ 
date. Inflation can end abruptly by a phase transition, or can last 
a long time if the potential becomes very flat after the observ¬ 
able region has been crossed. Moreover, there could be a short 
inflationary stage responsible for the origin of observable cos¬ 
mological perturbations, and another inflationary stage later on 
(but before nucleosynthesis), thus contributing to the total N^- 

In section 6 of PCI13, we performed this analysis with a full 
integration of the inflaton and metric perturbation modes, so that 
no slow-roll approximation was made. The only assumption was 
that primordial scalar perturbations are generated by the fluctu¬ 
ations of a single inflaton field with a canonical kinetic term. 
Since, in this approach, one is only interested in the potential 
over a narrow range of observable scales (centered around the 
held value (p, when the pivot scale crosses the Hubble radius), it 
is reasonable to test relatively simple potential shapes, described 
by a small number of free parameters. 

This approach gave very similar results to calculations based 
on the standard slow-roll analysis. This agreement can be ex¬ 
plained by the fact that the Planck 2013 data already preferred a 
primordial spectrum very close to a power law, at least over most 
of the observable range. Hence the 2013 data excluded strong 
deviations from slow-roll inflation, which would either produce 
a large running of the spectral index or imprint more compli¬ 
cated features on the primordial spectrum. However, this con¬ 
clusion did not apply to the largest scales observable by Planck, 
for which cosmic variance and slightly anomalous data points 
remained compatible with significant deviations from a simple 
power-law spectrum. The most striking result in section 6 of 
PCI 13 was the fact that, when giving enough freedom to the 
functional form of the inflation potential, the results were com¬ 
patible with a rather steep potential at the beginning of the ob¬ 
servable window, leading to “not-so-slow” roll during the first 
few observable e-folds. This explains the shape of the potential 
in figure 14 of PCI 13 for a Taylor expansion at order n = 4 and 
in the region where (p- p, < -0.2. However, such features were 
only partially explored because the method used for potential re¬ 
construction did not allow for an arbitrary value of the inflation 
velocity (p at the beginning of the observable window. Instead, 
our code imposed that the inflaton already tracked the inflation¬ 
ary attractor solution when the largest observable modes crossed 
the Hubble scale. 

Given that the Planck 2015 data establish even stronger con¬ 
straints on the primordial power spectrum than the 2013 re¬ 
sults, it is interesting to revisit the reconstruction of the po¬ 
tential V{(p). Section 7.1 presents some new results follow¬ 
ing the same approach as in PCI 13 (explained previously in 
Lesgourgues & Valkenburg (2007) and Mortonson et al. (201 1)). 
But in the present work, we also present some more general 
results, independent of any assumption concerning the initial 
held velocity p when the inflaton enters the observable win¬ 
dow. Following previous studies (Kinney, 2002; Kinney et al., 
2006; Peiris & Easther, 2006a; Easther & Peiris, 2006; Peiris 



Fig. 13. Posterior distributions for the first four potential slow- 
roll parameters when the potential is Taylor-expanded to 
nth order, using Planck TTh-IowPh-BAO (filled contours) or 
TT,TE,EEh-1owP (dashed contours). The primordial spectra are 
computed beyond any slow-roll approximation. 


& Easther, 2006b, 2008; Lesgourgues et al., 2008; Powell &, 
Kinney, 2007; Hamann et al., 2008; Norena et al., 2012), we re¬ 
construct the Hubble function H{(p), which determines both the 
potential, V{(p), through 

V{CP) = ml, hH4>) - 2Ml, [H'ipt . (60) 

and the solution p(t), through 

P = -2Ml,H'(p) , (61) 

with H'{(p) - dHld(p. Note that these two relations are exact. In 
Sect. 7.2, we fit H{(p) directly to the data, implicitly including 
all canonical single-field models in which the inflaton is rolling 
not very slowly (e not much smaller than one) just before en¬ 
tering the observable window, and the issue of having to start 
sufficiently early in order to allow the initial transient to decay 
is avoided. The only drawback in reconstructing H{(p) is that 
one cannot systematically test the most simple analytic forms 
for V{(p) in the observable range (for instance, polynomials of 
order n - 1,3,5,... in (0 - p^)). But our goal in this section is 
to check how much one can deviate from slow-roll inflation in 
general, independently of the shape of the underlying inflation 
potential. 

7. 1. Reconstruction of a smooth inflation potentiai 

Eollowing exactly the approach of PCI 13, we Taylor-expand the 
inflaton potential around (p - p^ to order n - 2, 3,4. To obtain 
faster converging Markov chains, instead of imposing flat pri¬ 
ors on the Taylor coefficients {V, V^, ..., we sample the 

potential slow-roll (PSR) parameters {ey,? 7 y,^y, zrr^), related to 
the former as indicated in Table 2. We stress that this is just a 
choice of prior, and does not imply that we are using any kind 
of slow-roll approximation in the calculation of the primordial 
spectra. 

The results are given in Table 7 (for Planck TTh-IowPh-BAO) 
and Pig. 13 (for the same data set, and also for Planck 
TT,TE,EEh-1owP). The second part of Table 7 shows the corre¬ 
sponding values of the spectral parameters ns, dnsfdlnk, and 
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Table 7. Numerical reconstruction of the potential slow- 
roll parameters beyond any slow-roll approximation, when 
the potential is Taylor-expanded to nth order, using Planck 
TTh-IowPh-BAO. We also show the corresponding bounds on 
some related parameters (here Us, dn^ldlnk, and ro .002 are de¬ 
rived from the numerically computed primordial spectra). All er¬ 
ror bars are at the 95 % CL. The effective value and Bayesian 
evidence logarithm (In B) of model n are given relative to model 
n - 1 (assuming flat priors for and in the range [-1,1]). 


^ 0 . 002 , as measured for each numerical primordial spectrum (at 
the pivot scales k - 0.05 Mpc ' for the scalar and 0.002 Mpc“^ 
for the tensor spectra), as well as the reionization optical depth. 
We also show in Fig. 14 the derived distribution of each coeffi¬ 
cient Vi (with a non-flat prior), and in Fig. 15 the reconstructed 
shape of the best-fit inflation potentials in the observable win¬ 
dow. Finally, the posterior distribution of the derived parameters 
ro.oo 2 and n^ is displayed in Fig. 16. 

Figure 13 shows that bounds are very similar when temper¬ 
ature data are combined with either high-^ polarization data or 
BAO data. This gives a hint of the robustness of these results. For 
both data sets, the error bars on the PSR parameters are typically 
smaller by a factor of 1.5 than in PCI13. 

Since potentials with n - 2 cannot generate a significant rvin- 
ning, the bounds on (Uj, ro. 002 ) and the best-fit models are very 
similar to those obtained with the ACDM-rr model in Sect. 5 and 
Table 4. On the other hand, in the n - 3 model, results follow 
the trend of the previous ACDM+r+dn^/dlnk analysis: the data 
prefer potentials with and of the same sign, generating 
a significant negative running (as can be seen in Fig. 16). This 
trend for is driven by the fact that a scalar spectrum with 
negative running reduces the power on large scales, and provides 
a better fit to low-^ temperature multipoles. However such a run¬ 
ning also suppresses power on small scales, so cannot be too 
large. 

The n = 4 case possesses a new feature: the potential 
has more freedom to generate complicated shapes which would 
roughly correspond to a running of the running of the tilt (as 
studied in Sect. 4). The best-fit models now have and 
of opposite sign, and a large positive The preferred com¬ 

bination of these parameters allows for even more suppression 
of power on large scales, while leaving small scales nearly un¬ 
changed. This can be seen very well in the shape of the scalar pri¬ 
mordial spectrum corresponding to the best-fit models, for both 
data sets Planck TTh-IowPh-BAO and Planck TT,TE,EEh-1owP. 



10“ Vo 10“ Vi 10“ V2 10“ V3 


Fig. 14. Posterior distribution for the coefficients of the infla¬ 
tion potential Taylor-expanded to nth order, in natural units 
(where VSwMpi = 1), reconstructed beyond any slow-roll ap¬ 
proximation, using Planck TTh-IowPh-BAO (filled contours) or 
TT,TE,EEh-1owP (dashed contours). The plot shows only half 
of the results; the other half is symmetric, with opposite signs 
for V^ and V^^^. Note that, unlike Pig. 13, the parameters shown 
here do not have flat priors, since they are mapped from the slow- 
roll parameters. 


These two best-fit models are very similar, but in Pig. 17 we 
show the one for Planck TT,TE,EEh-1owP, for which the trend is 
even more pronounced. Interestingly, the preferred models are 
such that power on large scales is suppressed in the scalar spec¬ 
trum and balanced by a small tensor contribution, of roughly 
^ 0.002 ~ 0.05. This particular combination gives the best fit to the 
low-^ data, shown in Pig. 18, while leaving the high-f tempera¬ 
ture spectrum identical to the best fit of the base ACDM model. 
Inflation produces such primordial perturbations with the fam¬ 
ily of green potentials displayed in Pig. 15. At the beginning of 
the observable range, the potential is very steep (ev(0) decreases 
from 0(1) to 0(10“^)), and produces a low amplitude of curva¬ 
ture perturbations (allowing a rather large tensor contribution, 
up to ro.oo 2 ~ 0.3). Then there is a transition towards a second 
region with a much smaller slope, leading to a nearly power-law 
curvature spectrum with the usual tilt value n^ =» 0.96. In Pig. 15, 
one can check that the height of the n - A potentials varies in a 
definite range, while the n - 2 and 3 potentials can have arbitrar¬ 
ily small amplitude at the pivot scale, reflecting the posteriors on 
ey or r. 

However, the improvement in x\g between the base ACDM 
and n = 4 models is only 2.2 (for Planck TTh-IowPh-BAO) or 4.3 
(for Planck TT,TE,EEh-1owP). This is very marginal and brings 
no significant evidence for these complicated models. This con¬ 
clusion is also supported by the calculation of the Bayesian evi¬ 
dence ratios, shown in the last line of Table 7 (under the assump¬ 
tion of flat priors in the range [-1,1 ] for f ^ and uZy): the evidence 
decreases each time that a new free parameter is added to the po¬ 
tential. At the 95 % CL, ro .002 is still compatible with zero, and 
so are the higher order PSR parameters and m\,. More free¬ 
dom in the inflaton potential allows fitting the data better, but 
under the assumption of a smooth potential in the observable 
range, a simple quadratic form provides the best explanation of 
the Planck observations. 
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Fig. 15. Observable range of the best-fit inflaton potentials, when 
V{4>) is Taylor expanded to the nth order around the pivot value 
0», in natural units (where VSTrMpi = 1), assuming a flat prior on 
ey, Tjv, and and using Planck TTh-IowPh-BAO. Potentials 
obtained under the transformation ((f) - (pt) (<^» - <p) leave the 
same observable signature and are also allowed. The sparsity of 
potentials with a small Vq = V((f>t) comes from the flat prior on 
ey rather than on In(yo); in fact, Vq is unbounded from below in 
the n = 2 and 3 results. The axis ranges are identical to those in 
Fig. 20, to make the comparison easier. 


With the Planck TTh-IowPh-BAO and TT,TE,EEh-1owP 
datasets, models with a large running or running of the running 
can be compatible with an unusually large value of the optical 
depth, as can be seen in Table 7. Including lensing informa¬ 
tion allows breaking the degeneracy between the optical depth 
and the primordial amplitude of scalar perturbations. Hence the 
Planck lensing data could be used to strengthen the conclusions 
of this section. 

Since in the n - A model, slow roll is marginally satisfied at 
the beginning of observable inflation, the reconstruction is very 
sensitive to the condition that there is an attractor solution at that 
time. Hence this case can in principle be investigated in a more 
conservative way using the P[((p) reconstruction method of the 
next section. 


7.2. Reconstruction of a smooth Hubble function 


In this section, we assume that the shape of the function H(cf>) 
is well captured within the observable window by a polynomial 
of order n (corresponding to a polynomial inflaton potential of 
order 2n): 


(j)' 


(=0 


(62) 


We vary n between 2 and 4. To avoid parameter degeneracies, as 
in the previous section we assume flat priors not on the Taylor 
coefficient //,, but on the Hubble slow-roll (HSR) parameters, 
which are related according to 


en = 2M; 
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Fig. 16. Posterior distribution for the tensor-to-scalar ratio (at 
k - 0.002Mpc“^) and for the running parameter dn^ldXnk 
(at k - 0.05 Mpc“^), for the potential reconstructions in 

Sects. 7.1 and 7.2. The Virp) reconstruction gives the solid curves 
for Planck TTh-IowPh-BAO, or dashed for TT,TE,EEh-1owP. 
The H((p) reconstruction gives the dotted curves for Planck 
TTh-IowPh-BAO, or dashed-dotted for TT,TE,EEh-1owP. The 
tensor-to-scalar ratio appears as a derived parameter, but by tak¬ 
ing a flat prior on either ey or en, we implicitly also take a nearly 
flat prior on r. The same applies to dn^ldink. 
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3 

4 
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< 0.020 
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0-082-ao36 

u.u^u_oo43 

0 096^^ 

W.S 

'-'•^^^■^-0.0093 


'-'•7U/_qqjq 

\ri7dn^ld\nk 

-0 25 

_13-fl8 

A-2’_i9 

0+21 

°-21 

ro.oo 2 

< 0.11 

< 0.16 

< 0.32 

Akeff 


= -0-6 

A43 = -2-3 


Table 8. Numerical reconstruction of the Hubble slow-roll pa¬ 
rameters beyond any slow-roll approximation, using Planck 
TTh-IowPh-BAO. We also show the corresponding bounds on 
some related parameters (here n^, dn^ldtnk, and ro.oo 2 are de¬ 
rived from the numerically computed primordial spectra). All 
error bars are at the 95 % CL. The effective value of model n 
is given relative to model n - 1. 


This is just a choice of prior. This analysis does not rely on the 
slow-roll approximation. 

Table 8 and Eig. 19 show our results for the reconstructed 
HSR parameters. Eigure 20 shows a representative sample of 
potential shapes V((p - (p,) derived using Eq. (60), for a sample 
of models drawn randomly from the chains, for the three cases 
n = 1,2,3. 
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Fig. 17. Primordial spectra (scalar and tensor) of the best-fit V{(p) 
model with n - A, for the Planck TT,TE,EEh-1owP data set, com¬ 
pared to the primordial spectrum (scalar only) of the best-fit base 
ACDM model. The best-ht potential is initially very steep, as can 
be seen in Pig. 15 (note the typical shape of the green curves). 
The transition from “marginal slow roll” (ev(0) between 0.01 
and 1) to “full slow roll” {ey{^) of order 0.01 or smaller) is re¬ 
sponsible for the suppression of the large-scale scalar spectrum. 



Fig. 19. Posterior distributions for the hrst four Hubble slow- 
roll parameters, when H{(p) is Taylor-expanded to nth order, us¬ 
ing Planck TTh-IowPh-BAO (filled contours) or TT,TE,EEh-1owP 
(dashed contours). The primordial spectra are computed beyond 
any slow-roll approximation. 



Fig. 20. Same as Fig. 15, when the Taylor expansion to nth or¬ 
der is performed on H{(p) instead of F(0), and the potential is 
inferred from Eq. (60). 


Fig. 18. Temperature spectrum (total, scalar contribution, ten¬ 
sor contribution) of the best-ht V{(p) model with n - 4, for the 
Planck TT,TE,EEh-1owP data set, compared to the temperature 
spectrum (scalar contribution only) of the best-ht base model. 
We also show the Planck low-^ temperature data, which is driv¬ 
ing the small differences between the two best-ht models. 


Most of the discussion of Sect. 7.1 also applies to this 
section, and so will not be repeated. Results for Planck 
TTh-IowPh-BAO and TT,TE,EEh-1owP are still very similar. The 
n - 2 case still gives results close to ACDM-tr, and the n - 3 
case to ACDM+r+dn^i/dlnk. The type of potential preferred in 
the n = 4 case is very similar to the n = 4 analysis of the previous 
section, for the reasons explained in Sect. 7.1. There are, how¬ 
ever, small differences, because the range of parametric forms 
for the potential explored by the two analyses differ. In the H{(p) 


reconstruction, the underlying potentials V((f)) are not polynomi¬ 
als. In the hrst approximation, they are close to polynomials of 
order 2n, but with constraints between the various coefficients. 
The main two differences with respect to the results of Sect. 7.1 
are as follows: 

- The reconstructed potential shape for n = 4 at the beginning 
of the observable window. Figure 20 shows that even steeper 
potentials are allowed than for the V((f>) method, with an 
even greater excursion of the inhaton held between Hubble 
crossing for the largest observable wavelengths and the pivot 
scale. This is because the H(<p) reconstruction does not rely 
on attractor solutions and automatically explores all valid po¬ 
tentials regardless of their initial held velocity. 

- The best-ht models are different, since they do not explore 
the same parametric families of potentials. In particular, for 
n = 4, the best-ht models have a negligible tensor contri- 














































Planck Collaboration: Constraints on inflation 


23 


bution, but still have a thick distribution tail towards large 
tensor-to-scalar ratios, so that the upper bound on ro .002 is as 
high as in the previous n - 4 models, ro .002 < 0.32. 


Note that tir^ can be significantly larger than unity for n - 4. 
This does not imply violation of slow roll within the observable 
range. By assumption, for all accepted models, en must remain 
smaller than unity over that range. In fact, for most of the green 
potentials visible in Fig. 20, we checked that en either has a max¬ 
imum very close to unity near the beginning of the observable 
range or starts from unity. So the best-fit models (maximizing 
the power suppression at low multipoles) correspond either to 
inflation of short duration, or to models nearly violating slow 
roll just before the observable window. However, such peculiar 
models are not necessary for a good fit. Table 8 shows that the 
improvement in^^ as n increases is negligible. 

In summary, this section further establishes the robust¬ 
ness of our potential reconstruction and two main conclusions. 
Firstly, under the assumption that the inflaton potential is smooth 
over the observable range, we showed that the simplest para¬ 
metric forms (involving only three free parameters including 
the amplitude no deviation from slow roll, and nearly 

power-law primordial spectra) are sufficient to explain the data. 
No high-order derivatives or deviations from slow roll are re¬ 
quired. Secondly, if one allows more freedom in the potential— 
typically, two more parameters—it is easy to decrease the large- 
scale primordial spectrum amplitude with an initial stage of 
“marginal slow roll” along a steep branch of the potential fol¬ 
lowed by a transition to a less steep branch. This type of 
model can accommodate a large tensor-to-scalar ratio, as high 
as ro.oo2 ~ 0.3. 


8. Pik) reconstruction 

In PCI 13 (section 7) we presented the results of a penalized like¬ 
lihood reconstruction, seeking to detect any possible deviations 
from a homogeneous power-law form (i.e., P(k) oc for the 
primordial power spectrum (PPS) for various values of a smooth¬ 
ing parameter, A. In the initial March 2013 preprint version of 
that paper, we reported evidence for a feature at moderate sta¬ 
tistical significance around k 0.15 Mpc However, in the 
November 2013 revision we retracted this finding, because sub¬ 
sequent tests indicated that the feature was no longer statistically 
significant when more aggressive cuts were made to exclude sky 
survey rings where contamination from electromagnetic inter¬ 
ference from the 4-K cooler were largest, as indicated in the 
November 2013 “Note Added.” 

In this section we report on results using the 2015 Cjt like¬ 
lihood (Sect. 8.1) using essentially the same methodology as de¬ 
scribed in PCI13. (See Gauthier & Bucher (2012) and references 
therein for more technical details.) This method is also extended 
to include the EE and TE likelihoods in Sect. 8.1.2. As part of 
this 2015 release, we include the results of two other methods 
(see Sects. 8.2 and 8.3) to search for features. We find that all 
three methods yield broadly consistent reconstructions and reach 
the following main conclusion: there is no statistically signifi¬ 
cant evidence for any features departing from a simple power- 
law (i.e., Pik) cx PPS. Given the substantial differences 

between these methods, it is satisfying to observe this conver¬ 
gence. 


8.1. Method I: penalized likelihood 

8.1.1. Update with 2015 temperature likelihood 

We repeated the same maximum likelihood analysis used to re¬ 
construct the PPS in PCI 13 using the updated Planck TT like¬ 
lihood. Since we are interested in deviations from the nearly 
scale invariant model currently favoured by the parametric ap¬ 
proach, we replaced the true PPS Puik) by a fiducial power-law 
spectrum P^^{k) = Aj(k/k,)”* *“', modulated by a small deviation 
function f(k): 

pK(k)^P^°\k)exp[f(k)]. (65) 


The deviation function f(ky was represented by B-spline ba¬ 
sis functions parameterized by control points f = 
which are the values of f(k) along a grid of knot points a:, = In kj. 

Naively, maximizing the Planck likelihood with respect to 
f results in over-fitting to cosmic variance and noise in the 
data. Furthermore, due to the limited range of scales over which 
Planck measures the anisotropy power spectrum, the likelihood 
is very weakly dependent on f{k) at extremely small and large 
scales. To address these issues, the following two penalty func¬ 
tions were added to the Planck likelihood: 


(^) 

X k'min pA-oo 

Ak /^(k) + a I d/c f^(k). 


( 66 ) 


The first term on the right-hand side of Eq. (66) is a rough¬ 
ness penalty, which disfavours /(k) that “wiggle” too much. The 
last two terms drive /(/c) to zero for scales below and above 
^max- The values of A and a represent the strengths of the re¬ 
spective penalties. The exact value of a is unimportant as long 
as it is large enough to drive /(k) close to zero on scales outside 
['('min, R'max]- However, the magnitude of the roughness penalty. A, 
controls the smoothness of the reconstruction. 

Since the anisotropy spectrum depends linearly on the PPS, 
the Newton-Raphson method is well suited to optimizing with 
respect to f. However, a maximum likelihood analysis also 
has to take into account the cosmological parameters, 0 s 
{//o, Qb/j^, These additional parameters are not easy to 

include in the Newton-Raphson method since it is difficult to 
evaluate the derivatives dCcld&, d^C(ld&^, etc., to the accuracy 
required by the method. Therefore a non-derivative method, such 
as the downhill simplex algorithm, is best suited to optimiza¬ 
tion over these parameters. Unfortunately the downhill simplex 
method is inefficient given the large number of control points in 
our parameter space. Since each method has its drawbacks, we 
combined the two methods to draw on their respective strengths. 
We define the function Ai as 

M(0)= min {-21ni:(0,f)-i-f^R(/l,a)f}. (67) 


Given a set of non-PPS cosmological parameters 0, M. is the 
value of the penalized log likelihood, minimized with respect 
to f using the Newton-Raphson method. The function M. is in 

’ The definition of f(k) used here differs from that of PCI 13 in that 
exp(/) is used in place of 1 - 1 - /, to ensure that the reconstructed primor¬ 
dial power spectrum is always non-negative. 

* Due to the high correlation between t and Aj, t is not included as a 
free parameter. Any change in t can be almost exactly compensated for 
by a change in Aj. We fix t to its best-fit fiducial model value. 
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Fig. 21. Planck TT likelihood primordial power spectrum (PPS) 
reconstruction results. Top four panels: Reconstruction of the de¬ 
viation f{k) using four different roughness penalties. The red 
curves represent the best-fit deviation f{k) using the Planck TT 
likelihood, fik) - 0 would represent a perfectly featureless spec¬ 
trum with respect to the hducial PPS model, which is obtained 
from the best-ht base ACDM model with a power-law PPS. The 
vertical extent of the dark and light green error bars indicates 
the +1 cr and +2 cr errors, respectively. The width of the error 
bars represents the minimum reconstructible width (the mini¬ 
mum width for a Gaussian feature so that the mean square devi¬ 
ation of the reconstruction is less than 10%). The grey regions 
indicate where the minimum reconstructible width is undefined, 
indicating that the reconstruction in these regions is untrustwor¬ 
thy. The hatched region in the T = 10^ plot shows where the 
fixing penalty has been applied. These hashed regions are not 
visible in the other three reconstructions, for which /<■„;„ lies out¬ 
side the range shown in the plots. For all values of the roughness 
penalty, all data points are within 2 cr of the hducial PPS ex¬ 
cept for the deviations around k 0.002 Mpc * in the A - 10^ 
and A - 10^ reconstructions. Lower three panels: ±lcr error 
bars of the three non-PPS cosmological parameters included in 
the maximum likelihood reconstruction. All values are consis¬ 
tent with their respective best-ht hducial model values indicated 
by the dashed lines. 


turn minimized with respect to 0 using the downhill simplex 
method. In contrast to the analysis done in PCI13, the Planck 
\qw-{ likelihood has been modihed so that it can be included 
in the Newton-Raphson minimization. Thus the reconstructions 
presented here extend to larger scales than were considered in 
2013. 


Fig. 22. Planck TT,TE,EEh-1owTEB likelihood primordial 
power spectrum reconstruction results. Top four panels: 
Reconstruction of the deviation f{k) using four different rough¬ 
ness penalties. As in Eig. 21, the red curves represent the best- 
ht deviation f{k) and the height and width of the green error 
bars represent the error and minimum reconstructible width, re¬ 
spectively. Eor all values of the roughness penalty, the deviations 
are consistent with a featureless spectrum. Lower three panels: 
+1 cr error bars of the three non-PPS cosmological parameters 
included in the maximum likelihood reconstruction. All values 
are consistent with their respective best-ht hducial model values 
(indicated by the dashed lines). 


Eigure 21 shows the best-ht PPS reconstruction using the 
Planck TTh-IowTEB likelihood. The penalties in Eq. (66) in¬ 
troduce a bias in the reconstruction by smoothing and other¬ 
wise deforming potential features in the power spectrum. To as¬ 
sess this bias, we dehne the “minimum reconstructible width” 
(MRW) to be the minimum width of a Gaussian feature needed 
so that the integrated squared difference between the feature and 
its reconstruction is less than 1 % of the integrated square of 
the input Gaussian, which is equivalent to 10 % rms. Due to the 
combination of the roughness and hxing penalties, it is impos¬ 
sible to satisfy the MRW criterion too close to /Cmin and 
Wherever the MRW is undefined, the reconstruction is substan¬ 
tially biased and therefore suspect. An MRW cannot be defined 
too close to the endpoints /Cmin and for two reasons: (1) 
lack of data; and (2) if a feature is too close to where the hx¬ 
ing penalty has been applied, the hxing penalty distorts the re¬ 
construction. Consequently a larger roughness penalty decreases 
the range over which an MRW is well dehned. The grey shaded 
areas in Eig. 21 show where the MRW is undehned and thus the 
reconstruction cannot be trusted. The cutoffs A-min and /Cmax have 
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been chosen to maximize the range over which an MRW is de- 
hned for a given value of A. The 1 cr and 2 cr error bars in Fig. 21 
are estimated using the Hessian of the log-likelihood evaluated 
at the best-ht PPS reconstruction. More detailed MC investiga¬ 
tions suggest that the nonlinear corrections to these error bars 
are small. 

For the A - 10^ and 10^ cases of the TT reconstruction, no 
deviation exceeds 2 cr, so we do not comment on the probability 
of obtaining a worse ht. For the other cases, we use the maxi¬ 
mum of the deviation, expressed in cr, of the plotted points as a 
metric of the quality of ht; then using Monte Carlo simulations 
we compute the p-value, or the probability to obtain a worse 
ht, according to this metric. For A - 10^ and 10"^, we obtain p- 
values of 0.304 and 0.142, respectively, corresponding to 1.03 
and 1.47 cr. We thus conclude that the observed deviations are 
not statistically signihcant. 

8.1.2. Penalized likelihood results with polarization 

In Fig. 22 the best-ht reconstruction of the PPS from 
the Planck TT,TE,EE likelihood is shown. We observe that 
the reconstruction including polarization broadly agrees with 
the reconstruction obtained using temperature only. Eor the 
Planck TT,TE,EE likelihood, we obtain for T = 10^, 10"^, and 
10^ the p-values 0.166, 0.107, and 0.045, respectively, corre¬ 
sponding to 1.38, 1.61, and 2.00 cr, and likewise conclude the 
absence of any statistically signihcant evidence for deviations 
from a simple power-law scalar primordial power spectrum. 

8.2. Method ll: Bayesian model comparison 

In this section we model the PPS 'P-R{k) using a nested family of 
models where 'P-R{k) is piecewise linear in the ln(!P)-ln(k) plane 
between a number of knots, A^knots, that is allowed to vary. The 
question arises as to how many knots one should use, and we ad¬ 
dress this question using Bayesian model comparison. A family 
of priors is chosen where both the horizontal and vertical posi¬ 
tions of the knots are allowed to vary. We examine the “Bayes 
factor” or “Bayesian evidence” as a function of Aknots to decide 


log'pR(k) 



Fig. 23. Linear spline reconstruction. The primordial power 
spectrum is reconstructed using Aknots interpolation points 
{{ki, Vj) : i - 1,2,... A^knotsl- The end knots are hxed in k but al¬ 
lowed to vary in V, whereas the internal knots can vary subject to 
the constraint that k\ < k 2 < ■ ■ ■ < The function PKik) is 

constructed within the range [ki, by interpolating logarith¬ 
mically between adjacent knots (i.e., linearly in log-log space). 
Outside this range the function is extrapolated logarithmically. 
The function P-Rik-, {ki,Pi]) thus has 2Aknots - 2 parameters. 



Parameter range 

Prior type 

10-^ Mpc^* 

V 

V 

II 

••• =0.3Mpc- 

' log uniform (sorted) 

2 < ln| 

10‘“'Pi),... 


log uniform 


2 < Aknots < 10 

integer uniform 


0.019 < < 0.025 

uniform 


0.095 < < 0.145 

uniform 


1.03 < lOO^Mc < 1.05 

uniform 


0.01 < T < 0.4 

uniform 


Table 9. Prior for moveable knot positions. The Pr posi¬ 
tions are distributed in a log-uniform manner across a wide 
range. The k positions are also log-uniformly distributed across 
the entire range needed by CosmoMC and are sorted so that 
ki < ... < When we marginalize over the number of 

knots, A^knots, we assume a uniform prior between 2 and 10. 


how many knots are statistically justihed. A similar analysis has 
been performed by Vazquez et al. (2012) and Aslanyan et al. 
(2014). In addition, we marginalize over all possible numbers of 
knots to obtain an averaged reconstruction weighted according 
to the Bayesian evidence. 

The generic prescription is illustrated in Pig. 23. Vknots knots 
{{ki,Pi): i - 1 ,..., Vknots) are placed in the (k,PR) plane and 
the function PR{k) is constructed by logarithmic interpolation (a 
linear interpolation in log-log space) between adjacent points. 
Outside the horizontal range {k \, the function is extrapolated 
using the outermost interval. 

Within this framework, base ACDM arises when 
A^knots = 2—in other words, when there are two boundary 
knots and no internal knots, and the parameters P\ and P 2 
(in place of As and n^) parameterize the simple power-law 
PPS. There are also, of course, the four standard cosmological 
parameters (Qb/r^, 1OO0 mc^ and r), as well the numerous 

foreground parameters associated with the Planck high-^ 
likelihood, all of which are unrelated to the PPS. This simplest 
model can be extended iteratively by successively inserting an 
additional internal knot, thus requiring with each iteration two 
more variables to parameterize the new knot position. 

We run models for a variety of numbers of internal knots, 
Mnt = A^knots - 2, evaluating the evidence for Vint. Under the as¬ 
sumption that the prior is justihed, the most likely, or preferred, 
model is the one with the highest evidence. Evidences are eval¬ 
uated using the PolyChord sampler (Handley et al., 2015) in 
CAMS and CosmoMC. The use of PolyChord is essential, as 
the posteriors in this parameterization are often multi-modal. 
Also, the ordered log-uniform priors on the k, are easy to im¬ 
plement within the PolyChord framework. All runs were per¬ 
formed with 1000 live points, oversampling the semi-slow and 
fast parameters by a factor of 5 and 100, respectively. 

Priors for the reconstruction and cosmological parameters 
are detailed in Table 8.2. We report evidence ratios with re¬ 
spect to the base ACDM case. The cosmological priors re¬ 
main the same for all models, and this part of the prior has al¬ 
most no impact on the evidence ratios. The choice of prior on 
the reconstruction parameters {Ph does affect the Bayes factor. 
CosmoMC, however, puts an implicit prior on all models by ex¬ 
cluding parameter choices that render the internal computational 
approximations in CAMB invalid. The baseline prior for the ver¬ 
tical position of the knots includes all of the range allowed by 
CosmoMC, so slighly increasing this prior range will not affect 
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Fig. 24. Bayesian movable knot reconstructions of the primordial power spectrum Pnik) using Planck XT data. The plots indicate 
our knowledge of the PPS P{V-R{k)\k,N) for a given number of knots. The number of internal knots Nint increases (left to right and 
top to bottom) from 0 to 8. For each k-slice, equal colours have equal probabilities. The colour scale is chosen so that darker regions 
correspond to lower-cr conhdence intervals. 1 <j and 2 <j conhdence intervals are also sketched (black curves). The upper horizontal 
axes give the approximate corresponding multipoles via £ ^ where Dkc is the comoving distance to recombination. 



Fig. 25. Bayes factor (relative to the base ACDM model) as a 
function of the number of knots for three separate runs. Solid 
line: Planck TT. Dashed line: Planck TT,TE,EE. Dotted line: 
Planck TT, with priors on the P parameters reduced in width 
by a factor of 2 (2.5 < ln(10“’!P) < 3.5). 
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Fig. 26. Bayesian reconstruction of the primordial power spec¬ 
trum averaged over different values of Mint (as shown in 
Fig. 24), weighted according to the Bayesian evidence. The re¬ 
gion 30 < £ < 2300 is highly constrained, but the resolution is 
lacking to say anything precise about higher £. At lower £, cos¬ 
mic variance reduces our knowledge of V-R{k). The weights as¬ 
signed to the lower Aint models outweigh those of the higher 
models, so no oscillatory features are visible here. 
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the evidence ratios. If one were to reduce the prior widths sig- 
nihcantly, the evidence ratios would be increased. The allowed 
horizontal range includes all k-scales accessible to Planck. Thus, 
altering this width would be unphysical. 

After completion of an evidence calculation, PolyChord 
generates a representative set of samples of the posterior for each 
model P{&} = /"(©Idata, Aint). We may use this to calculate a 
marginalized probability distribution for the PPS: 

P{\ogP<R\k,Niat)^ ^6(\ogPK-\ogPK(t,&))P{@)d@. ( 68 ) 

This expression encapsulates our knowledge of at each value 
of k for a given number of knots. Plots of this PPS posterior are 
shown in Fig. 24 using Planck TT data. 

If one considers the Bayesian evidences of each model. 
Fig. 25 shows that although no model is preferred over base 
ACDM, the case Aint = 1 is competitive. This model is anal¬ 
ogous to the broken-power-law spectrum of Sect. 4.4, although 
the models differ signihcantly in terms of the priors used. In this 
case, the additional freedom of one knot allows a reconstruc¬ 
tion of the suppression of power at low £. Adding polarization 
data does not alter the evidences signihcantly, although Aint = 1 
is strengthened. We also plot a Planck TT run, but with the re¬ 
duced vertical priors 2.5 < ln(lO'°!P) < 3.5. As expected, this 
increases the evidence ratios, but does not alter the above con¬ 
clusion. 

For increasing numbers of internal knots, the Bayesian evi¬ 
dence monotonically decreases. Occam’s razor dictates, there¬ 
fore, that these models should not be preferred, due to their 
higher complexity. However, there is an intriguing stable oscilla¬ 
tory feature, at 20 < f < 50, that appears once there are enough 
knots to reconstruct it. This is a qualitative feature predicted by 
several inhationary models (discussed in Sect. 9), and a possible 
hint of new physics, although its statistical signihcance is not 
compelling. 

A full Bayesian analysis marginalizes over all models 
weighted according to the normalized evidence so that 

P{\ogrn\k) - ^ P(\ogrn\k, Ai„t)Z^„„ (69) 

as indicated in Fig. 26. This reconstruction is sensitive to how 
model complexity is penalized in the prior distribution. 

8.3. Method III: cubic spline reconstruction 

In this section we investigate another reconstruction algorithm 
based on cubic splines in the ln(k)-ln!P?j plane, where (unlike 
for the approach of the previous subsection) the horizontal po¬ 
sitions of the knots are uniformly spaced in ln(A:) and fixed. A 
prior on the vertical positions (described in detail below) is cho¬ 
sen and the reconstructed power spectrum is calculated using 
CosmoMC for various numbers of knots. This method differs 
from the method in Sect. 8.1 in that the smoothness is controlled 
by the number of discrete knots rather than by a continuous 
parameter of a statistical model having a well-defined contin¬ 
uum limit. With respect to the Bayesian model comparison of 
Sect. 8.2, the assessment of model complexity differs because 
here the knots are not movable. 

Let the horizontal positions of the n knots be given by kh, 
where b - spaced so that kt+x/kh is independent of b. 

We single out a “pivot knot” b-p, so that kp - kt - 0.05 Mpc“*, 
which is the standard scalar power spectrum pivot scale. For 
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Fig. 27. Reconstructed power spectra applied to the Planck 2015 
data using 12 knots (with positions marked as A at the bottom 
of each panel) with cubic spline interpolation. Mean spectra as 
well as sample trajectories are shown for scalars and tensors, and 
+ 1 cr and +2cr limits are shown for the scalars. Top: uniform 
prior 0 < r < 1. Middle: hxed r = 0.1; Bottom: hxed r - 0.01. 
Data sets; Planck TTH-lowPH-BAOH-SNH-HSTH-Zre>6 prior. Diec is 
the comoving distance to recombination. 


a given number of knots n we choose k\ and k„ so that the 
interval of relevant cosmological scales, taken to extend from 
10 Mpc“^ to 0(1) Mpc ', is included. We now define the prior 
on the vertical knot coordinates. For the pivot point, we dehne 
InAs = where InAj has a uniformly distributed prior, 

and for the other points with b p,'Nt dehne the derived vari- 
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Fig. 28. Reconstructed 12-knot power spectra with polarization included. Data sets in common: lowPH-BAOH-SNH-HSTH-Zie>6 prior. 


able 


qb = In 


/ 'Pnih) \ 

yPuM^^b)] 


(70) 


where VnM = Here the spectral index risfid is 

fixed. A uniform prior is imposed on each variable qt, (b + p) 
and the constraint -1 < < 1 is also imposed to force the re¬ 

construction to behave reasonably near the endpoints, where it is 
hardly constrained by the data. The quantity InPuik) is interpo¬ 
lated between the knots using cubic splines with natural bound¬ 
ary conditions (i.e., the second derivatives vanish at the first 
and the last knots). Outside [ki,k„] we set 'Puik) = 

(for k < k\) and Puik) — e^"!P'K,fid(^) (for k > kn). For most 
knots, we found that the upper and lower bounds of the qt prior 
hardly affect the reconstruction, since the data sharpens the al¬ 
lowed range significantly. However, for superhorizon scales (i.e., 
k < 10“^ Mpc^') and very small scales (i.e., k > 0.2 Mpc^'), 
which are only weakly constrained by the cosmological data, the 
prior dominates the reconstruction. For the results here, a fidu¬ 
cial spectral index - 0.967 for 'PfiM was chosen, which is 
close to the estimate from Planck TTh-IowPh-BAO. A different 
choice of Ws.fid leads to a trivial linear shift in the qt. 

The possible presence of tensor modes (see Sect. 5) has the 
potential to bias and introduce additional uncertainty in the re¬ 
construction of the primordial scalar power spectrum as parame¬ 
terized above. Obviously, in the absence of a detection of tensors 
at high statistical significance, it is not sensible to model a pos¬ 
sible tensor contribution with more than a few degrees of free¬ 


dom. A complicated model would lead to prior-dominated re¬ 
sults. We therefore use the power-law parameterization, Pi{k) - 
rA^ikjkPf', where the consistency relation n, = -r/8 is enforced 
as a constraint. 

Primordial tensor fluctuations contribute to CMB temper¬ 
ature and polarization angular power spectra, in particular at 
spatial scales larger than the recombination Hubble length, k < 
(aH)i-sc ~ 0.005 Mpc *. If a large number of knots in lnP<ji(k) 
is included over that range, then a modified can mimic a ten¬ 
sor contribution, leading to a near-degeneracy. This can lead to 
large uncertainty in the tensor amplitude, r. Once r is measured 
or tightly constrained in B-mode experiments, this near degener¬ 
acy will be broken. As examples here, we do allow r to float, but 
also show what happens when r is constrained to take the values 
r = 0.1 and r — 0.01 in the reconstruction. 

Figure 27 shows the reconstruction obtained using the 2015 
Planck TTh-IowP likelihood, BAO, SNIa, HST, and a Zie > 6 
prior. Including these ancillary likelihoods improves the con¬ 
straint on the PPS by helping to fix the cosmological parame¬ 
ters (e.g., Hq, t, and the late-time expansion history), which in 
this context may be regarded as nuisance parameters. These re¬ 
sults were obtained by modifying CosmoMC to incorporate the 
n-knot parameterization of the PPS. Here 12 knots were used 
and the mean reconstruction as well as the 1 cr and 2 cr limits are 
shown. Some 1 cr sample trajectories (dashed curves) are also 
shown to illustrate the degree of correlation or smoothing of 
the reconstruction. The tensor trajectories are also shown, but, 
as explained above, have been constrained to be straight lines. 
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Fig. 29. Reconstructed power spectra with the base ACDM 
best-fit subtracted. The mean spectra shown are for the floating 
r and the two fixed r cases with 12 cubic spline knots. These 
should be contrasted with the running best-fit mean (green) and 
the similar looking uniform case in which t has been lowered 
from its best-fit base ACDM value to 0.04. Data points are the 
Planck 2015 Commander {{ < 30) and Plik {{ > 30) temper¬ 
ature power spectrum. 


In the top panel r is allowed to freely float, and a wide range 
of r is allowed because of the near-degeneracy with the low-k 
scalar power. Two illustrative values of fixed r (i.e., r - 0.1 and 
r = 0.01) are also shown, to give an idea of how much the re¬ 
construction is sensitive to variations in r within the range of 
presently plausible values. 

The reconstructions using the 2013 Planck likelihood in 
place of the 2015 likelihood are broadly consistent with the re¬ 
construction shown in Fig. 27. To demonstrate robustness with 
respect to the interpolation scheme we tried using linear interpo¬ 
lation instead of cubic splines and found that the reconstruction 
was consistent provided enough knots (i.e., Wknot ~ 14) were 
used. At intermediate k the reconstruction is consistent with a 
simple power law, corresponding to a straight line in Fig. 27. 
We observe that once k drops, so that the effective multipole 
being probed is below about 60, deviations from a power law 
appear, but the dispersion in allowed trajectories also rises as 
a consequence of cosmic variance. The power deficit dip at 
k ~ 0.002 Mpc“^ (i.e., 4 = ~ 30, where Djec is the co¬ 

moving distance to recombination) is largely driven by the power 
spectrum anomaly in the { 20-30 range that has been evident 

since the early spectra from WMAP (Bennett et ah, 2011), and 
verified by Planck. 

We also explore the impact of including the Planck polariza¬ 
tion likelihood in the reconstruction. Figure 28 shows the recon¬ 
structed power spectra using various combinations of the polar¬ 
ization and temperature data. The f < 30 treatments are the same 
in all cases, so this is mainly a test of the higher k region. What 
is seen is that, except at high k, the EE polarization data also 
enforce a nearly uniform n^, consistent with that from TT, over 
a broad k-range. When rF is used alone, or TE and EE are used 
in combination, the result is also very similar. The upper right 
panel shows the constraints from all three spectra together, and 
the errors on the reconstruction are now better than those from 
TT alone. 

It is interesting to examine how the TT power spectrum ob¬ 
tained using the above reconstructions compares to the CMB 
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Fig. 30. The degeneracy between t and the knot variables and 
q 4 in the 12-knot case shown in Fig. 27. 


data, in particular around the range £ =» 20-30, corresponding 
roughly to k^ ^ 1.5 x 10“^ Mpc^'. In Fig. 29 the differences in 
from the best-fit simple power-law model are plotted for 
various assumptions concerning r. We see that a better fit than 
the power-law model can apparently be obtained around £ xi 20- 
30. We quantify this improvement below. 

Due to the degeneracy of scalar and tensor contributions to 
the significance of the low-f” anomaly depends on the ten¬ 
sor prior and whether polarization data are used. For k < 10“^, 
once more degeneracy appears: the shape of also depends 
on the reionization optical depth, t. In Fig. 29 we also show the 
effect of replacing the best-fit t for tilted base ACDM with a 
low value, while keeping unchanged. A low t bends 

downward at 7” < 10. For the 12-knot (or similar) runs, if t is al¬ 
lowed to run into (nonphysically) small values t < 0.04, a slight 
rise in P^ik) at k ~ 3 x 10“"^Mpc“* is preferred to compensate 
the low-T effect. This degeneracy can be broken to a certain ex¬ 
tent using low-redshift data: Zw > 6 from quasar observations 
(Becker et ah, 2001), BAG (SDSS), Supernova (JLA), and HST 
It is evident that allowing n^ to run is not what the data 
prefer. The best-fit running is also shown in Fig. 29. The k-space 
!PK(k)-response shown in Fig. 27 shows that running does not 
capture the shape of the low-^ residuals. 

We have shown that the cubic spline reconstruction studied 
in this section consistently produces a dip in q 4 , corresponding to 
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Table 10. Reducedand /?-values for low-A: knots (5 knots) and high-k knots (6 knots, pivot knot excluded), with the null hypothesis 
being the best-fit power law spectrum. Low-z data refers to BAOH-SNH-HSTH-Zre>6 prior. In all cases lowP data are used. 


r prior 

low-z data 

Planck data 

loW-^T^uced 

low-k p-value 

high-k ^r-Luced 

high-k p-value 

q-i constraint 

^4 constraint 

0 < r < 1 

used 

TT 

0.95 

0.45 

0.17 

0.98 

-0.07 ± 0.28 

-0.39 ± 0.20 

r = 0.01 

used 

TT 

1.13 

0.34 

0.09 

0.997 

0.01 ± 0.24 

-0.23 ±0.12 

r = 0.01 

not used 

TT 

0.89 

0.49 

0.36 

0.90 

0.10 ±0.24 

-0.23 ±0.12 

r = 0.1 

used 

TT 

1.70 

0.13 

0.12 

0.994 

-0.04 ± 0.26 

-0.28 ±0.13 

r = 0.1 

not used 

TT 

1.46 

0.20 

0.38 

0.89 

0.05 ± 0.27 

-0.28 ±0.13 

r = 0.1 

used 

TT,TE,EE 

1.71 

0.13 

0.17 

0.985 

-0.02 ± 0.25 

-0.30 ±0.12 

r = 0.1 

used 

TE,EE 

1.72 

0.13 

0.38 

0.89 

0.06 ± 0.25 

-0.32 ±0.15 

r = 0.1 

used 

TE 

1.80 

0.11 

0.26 

0.95 

-0.02 ± 0.27 

-0.17 ±0.16 

r = 0.1 

used 

EE 

1.78 

0.11 

0.18 

0.98 

0.09 ± 0.25 

-0.39 ±0.16 

r = 0.1 

used 

TT-l-lensing 

1.54 

0.17 

0.05 

0.9995 

0.05 ± 0.25 

-0.27 ±0.13 





A:[Mpc 


Fig. 31. Slow-roll parameter e for reconstructed trajectories us¬ 
ing 12 knots (marked as A at the bottom of the figure) with cubic 
spline interpolation. The mean values are shown for floating r 
and r fixed to be 0.1 and 0.01. Sample 1 cr trajectories shown for 
the floating r case show wide variability, which is significantly 
diminished if r is fixed to r = 0.1, as shown. 



Fig. 32. Reconstructed single-field inflaton potentials from the 
cubic spline power spectra mode expansion using 12 knots. 


k ss 1.5 X 10“^ Mpc. We now turn to the question of whether this 
result is real or simply the result of cosmic variance. To assess 
the statistical significance of the departures of the mean recon¬ 
struction from a simple power law, we calculate the low-k and 
high-k reduced for the five qh values for scales below and six 
qi, values (b + p) for scales above SQ/Dj-g^, respectively, indicat¬ 
ing the corresponding p-values (i.e., probability to exceed), for 
various data combinations, in Table 10. The high-k fit is better 
than expected for reasons that we do not understand, but we at¬ 
tribute this situation to chance. The low-k region shows a poor fit, 
but in no case does the p-value fall below 10 %. Therefore, even 
though the low-k dip is robust against the various choices made 
for the reconstruction, we conclude that it is not statistically sig¬ 
nificant. The plot for the knot position of the dip (corresponding 
to q^) in Fig. 30 does not contradict this conclusion. 

Because of the r degeneracy associated with the scalar 
power, it is best when quoting statistics to use the fixed r 
cases, although for completeness we show the floating r case 
as well. There is also a smaller effect associated with the t 
degeneracy, and the values quoted have restricted the redshift 
of reionization to exceed 6. The value Zk = 6.5 was used in 
Planck Collaboration XIII (2015). The significance of the low-k 
anomaly is meaningful only if an explicit r prior and low-redshift 
constraint on t have been applied. 

Finally, we relate the reconstructed 'Puik) calculated above 
to the trajectories of the slow-roll parameter e - -Hj 
plotted as a function of k (see Fig. 31). We also plot in Fig. 32 the 
reconstructed inflationary potential in the region over which the 
inflationary potential is constrained by the data. Here canonical 
single-field inflation is assumed, and the value of r enters solely 
to fix the height of the potential at the pivot scale. This is not 
entirely self-consistent, but justified by the lack of constraining 
power on the tensors at present. 

8.4. Power spectrum reconstruction summary 

The three non-parametric methods for reconstructing the primor¬ 
dial power spectrum explored here support the following two 
conclusions: 

1. Except possibly at low k, over the range of k where the CMB 
data best constrains the form of the primordial power spec¬ 
trum, none of the three methods And any statistically signifi¬ 
cant evidence for deviations from a simple power law form. 
The fluctuations seen in this regime are entirely consistent 
with the expectations from cosmic variance and noise. 
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2. At low k, all three methods reconstruct a power deficit at 
k ai 1.5-2.0 X 10“^ Mpc which can be linked to the dip in 
the TT angular power spectrum at f « 20-30. This agree¬ 
ment suggests that the reconstruction of this “anomaly” is 
not an artefact of any of the methods, but rather inherent in 
the CMB data itself. However, the evidence for this feature 
is marginal since it is in a region of the spectrum where the 
fluctuations from cosmic variance are large. 


9. Search for parameterized features 

In this section, we explore the possibility of a radical depar¬ 
ture from the near-scale-invariant power-law spectrum, !P^(k) = 
As(k/k,)"*“\ of the standard slow-roll scenario for a selection of 
theoretically motivated parameterizations of the spectrum. 


Model 

Parameter 

Prior range 



[0,2] 

Step 

logic 

[-5,0] 


InXs 

[-1,5] 


^log 

[0,0.5] 

Log osc. 

logic ^log 

[0,2.1] 


<Plog 

[0,2;r] 


•^lin 

[0,0.5] 

Linear osc. 

log 10 tUlin 

[0,2] 

[-1,1] 


'Lin 


^lin 

[0,2;r] 

Cutoff 

logic (kc/Mpc^‘) 

[-5, -2] 


Table 11. Parameters and prior ranges. 


9.1. Models 


9.1.1. Step in the inflaton potential 


A sudden, step-like feature in the inflaton potential (Adams 
et ak, 2001) or the sound speed (Achucarro et al., 2011) will lead 
to a localized oscillatory burst in the scalar primordial power 
spectrum. A general parameterization describing both a tanh- 
step in the potential and in the warp term of a DBI model was 
proposed in Miranda & Hu (2014): 

InP^^ik) = exp [ln!P°(k) -r Io(k) + \n[l+ l]{k ))], (71) 

where the first- and second-order terms are given by 


Io^Jl,'Wo(klK)D 

1 


I X, 


Ji = 


V2 


-(l-ns) + [^sAVi(k/k,)+ 


^2 ^2(k/k,) + ^3 ^iWks)] D 


(“I 


with window functions 


'Wo(x) - |^(l8v - 6x^)cos2v-H (l5x^ - 9 ) sin2vj, (74) 


(72) 


(73) 


2x4 
1 
x4 


AVi(x) = {3(xcos X - sinx) [3xcos x -H (2x^ - 3) sinxjj, 

(75) 


^lix) - (sin X - X cos xy 
x4 


(76) 


- ^ [3 H- 2x^ “ (3 “ 4x^)cos(2x) - 6xsin(2x)j, 


(77) 


and damping function 

Dix) = (78) 

sinhx 

Due to the high complexity of this model, we focus on the limit¬ 
ing case of a step in the potential (^2 - -^3 = 0). 


2002; Bozza et ak, 2003), or, approximately, in the axion mon- 
odromy model, explored in more detail in Sect. 10. We assume 
a constant modulation amplitude and use 


=!P“(k)n +^,ogCos 


tuiog In 



(79) 


9.1.3. Linear oscillations 

A modulation linear in k can be obtained, e.g., in boundary effec¬ 
tive field theory models (Jackson & Shiu, 2013), and is typically 
accompanied by a scale-dependent modulation amplitude. We 
adopt the parameterization used in Meerburg & Spergel (2014), 
which allows for a strong scale dependence of the modulation 
amplitude: 


= !P«(k) 


1 -I- ,;?iii, 



cos 



(80) 


9.1.4. Cutoff model 


If today’s largest observable scales exited the Hubble radius be¬ 
fore the inflaton field reached the slow-roll attractor, the ampli¬ 
tude of the primordial power spectrum is typically strongly sup¬ 
pressed at low k. As an example of such a model, we consider 
a scenario in which slow roll is preceded by a stage of kinetic 
energy domination. The resulting power spectrum was derived 
by Contaldi et ak (2003) and can be expressed as 

InP^^ik) = lnf’°(k) + In ^ |Ce - Dd j , (81) 


with 



9.1.2. Logarithmic oscillations 

Logarithmic modulations of the primordial power spectrum 
generically appear, e.g., in models with non-Bunch-Davies ini¬ 
tial conditions (Martin & Brandenberger, 2001; Danielsson, 


(2^ 

where H;, denotes the Hankel function of the second kind. The 
power spectrum in this model is exponentially suppressed for 
wavenumbers smaller than the cutoff scale k^. and converges to 
a standard power-law spectrum for k :s> k^, with an oscillatory 
transition region for k> k^. 
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Model 

Planck TT-l-lowP 

Planck TT,TE,EE-l-lowP 

PTE 


InS 

ALeff 

InS 

Step 

-8.6 

-0.3 

-7.3 

-0.6 

0.09 

Log osc. 

-10.6 

-1.9 

-10.1 

-1.5 

0.24 

Linear osc. 

-8.9 

-1.9 

-10.9 

-1.3 

0.50 

Cutoff 

-2.0 

-0.4 

-2.2 

-0.6 

0.12 


Table 12. Improvement in fit and Bayes factors with respect 
to power-law base ACDM for Planck TT-rlowP and Planck 
TT,TE,EEh-1owP data, as well as approximate probability to ex¬ 
ceed the observed (p-value), constructed from simulated 
Planck TTh-IowP data. Negative Bayes factors indicate a prefer¬ 
ence for the power-law model. 
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Fig. 33. Distribution of from 400 simulated Planck 

TTh-IowP data sets. 


9.2. Analysis and results 

We use MultiNest to evaluate the Bayesian evidence for the 
models and establish parameter constraints, and to roughly iden¬ 
tify the global maximum likelihood region of parameter space. 
The features model best-ht parameters and InX are then ob¬ 
tained with the help of the CosmoMC minimization algorithm, 
taking narrow priors around the MultiNest best ht. Eor the 
parameters of the features models, we assign flat prior probabili¬ 
ties. The prior ranges are listed in Table 11 . Note that throughout 
this section, for the sake of maximizing sensitivity to very sharp 
features, the unbinned (“binl”) versions of the high-f TT and 
TT,TE,EE likelihoods are used instead of the standard, binned 
versions. 

Since the features considered here can lead to broad distor¬ 
tions of the CMB angular power spectrum degenerate with the 
late-time cosmological parameters (Miranda & Hu, 2014), in all 
cases we simultaneously vary primordial parameters and all the 
ACDM parameters, but we keep the foreground parameters hxed 
to their best-ht values for the power-law base ACDM model. 

We present the Bayes factors with respect to the power-law 
base ACDM model and the improvement over the power- 
law model in Table 12. Eor our choice of priors, none of the 
features models are preferred over a power-law spectrum. The 
best-ht power spectra are plotted in Eig. 34. While the cutoff and 
step model best hts reproduce the large-scale suppression 
20-30 also obtained by direct power spectrum reconstruction in 


Fig. 34. Best-ht power spectra for the power-law (black curve), 
step (green), logarithmic oscillation (blue), linear oscillation (or¬ 
ange), and cutoff (red) models using Planck TT-tlowP data. The 
brown curve is the best ht for a model with a step in the warp 
and potential (Eqs. (71)-(78)). 


Sect. 8, the oscillation models prefer relatively high frequencies 
that are beyond the resolution of the reconstruction methods. 

In addition to the four features models we also show in 
Fig. 34 the best ht of a model allowing for steps in both in- 
Haton potential and warp (brown line); note particularly the 
strong resemblance to the reconstructed features of the previ¬ 
ous section. The effective A;^f^ for this model is -12.1 (-11.5) 
for Planck TT-tlowP {Planck TT,TE,EEh-1owP) data, at the cost 
of adding hve new parameters, resulting in a In-Bayes factor of 
-0.8 (-0.4). A similar phenomenology can be also be found for 
the case of a sudden change in the slope of the inhaton potential 
(Starobinsky, 1992; Choe et ah, 2004), which yields a best-ht 
Ay^ - -4.5(-4.9) for two extra parameters. 

As shown in Table 13, constraints on the remaining cosmo¬ 
logical parameters are not signihcantly affected when allowing 
for the presence of features. 

For the cutoff and step models, the inclusion of Planck 
small-scale polarization data does not add much in terms of di¬ 
rect sensitivity; the best hts lie in the same parameter region 
as for Planck TT-tlowP data and the and Bayes factors 
are not subject to major changes. The two oscillation models’ 
Planck TTh-IowP best hts, on the other hand, also predict a 
non-negligible signature in the polarization spectra at high €. 
Therefore, if the features were real, one would expect an ad¬ 
ditional improvement in Ay^ for Planck TT,TE,EEh-1owP. This 
is not the case here. Though the linear oscillation model’s max¬ 
imum Ay^ does increase, the local Ay^ in the Planck TT-tlowP 
best-ht regions is in fact reduced for both models, and the global 
likelihood maxima occur at different frequencies (logjg wiog = 
1.25 and logiQtuiin = 1.02) compared to their Planck TT-tlowP 
counterparts. 

In addition to the Bayesian evidence analysis, we also ap¬ 
proach the matter of the statistical relevance of the features 
models from a frequentist statistics perspective, in order to give 
the Ay^^g numbers a quantitative interpretation. Assuming that 
the underlying V-Rik) was actually a featureless power law, we 
can ask how large an improvement to InX the different fea¬ 
tures models would yield on average, just by overhtting scat¬ 
ter from cosmic variance and noise. For this purpose, we sim¬ 
ulate Planck power spectrum data sets consisting of tempera- 
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Parameter 

Step 

Log osc. 

Linear osc. 

Cutoff 

Power law 

100 

2.23 ± 0.02 

2.22 ± 0.02 

2.23 ± 0.02 

2.23 ± 0.02 

2.23 ± 0.02 

lOtjj. 

1.20 ±0.02 

1.20 ±0.02 

1.20 ±0.02 

1.19 ±0.02 

1.19 ±0.02 

100 Smc 

1.0409 ± 0.0004 

1.0409 ± 0.0004 

1.0409 ± 0.0004 

1.0410 ±0.0005 

1.0409 ±0.0005 

T 

0.083 ±0.015 

0.082 ±0.015 

0.084 ±0.014 

0.086 ±0.017 

0.085 ±0.016 

In (10*%) 

3.10 ±0.03 

3.10 ±0.03 

3.10 ±0.03 

3.11 ±0.03 

3.10 ±0.03 


0.966 ± 0.005 

0.970 ± 0.007 

0.967 ± 0.004 

0.968 ± 0.005 

0.968 ± 0.005 


0.374 





logio (fcs/Mpc^‘) 

-3.10 





InXs 

0.342 





^log 


0.0278 




logiQ tt>iog 


1.51 




(^iog/2;r 


0.634 




^lin 



0.0292 



logio Wlin 



1.73 



^lin 



0.662 



Ain/2^ 



0.554 



logio (^c/Mpc^‘) 




-3.44 



Table 13. Best-fit features parameters and parameter constraints on the remaining cosmological parameters for the four features 
models for Planck TT-rlowP data. Note that the foreground parameters have been fixed to their power-law base ACDM best-fit 
values. 


ture and polarization up to ^ = 29 and unbinned temperature for 
30 < f < 2508, taking as input fiducial spectra the power-law 
base ACDM model’s best-fit spectra. 

For each of these simulated Planck data sets, we perform 
the following procedure: (i) find the power-law ACDM model’s 
best-fit parameters with CosmoMC’s minimization algorithm, 
(ii) hx the non-primordial parameters (wb, Wc, t) to their re¬ 
spective best-fit values, (iii) using MultiNest, hnd the best fit 
of the features models,'^ and (iv) extract the effective be¬ 
tween power-law and features models. 

The resulting distributions are shown in Fig. 33. Compared 
to the real data A;^^^ values from Table 12, they are biased to¬ 
wards lower values, since we do not vary the late-time cos¬ 
mological parameters in the analysis of the simulated data. 
Nonetheless, the observed improvements in the fit do not appear 
to be extraordinarily large, with the respective (conservative) p- 
values ranging between 0.09 and 0.50. 

These observations lead to the conclusion that even though 
some of the peculiarities seen in the residuals of the Planck data 
with respect to a power-law primordial spectrum may be ex¬ 
plained in terms of primordial features, none of the simple model 
templates considered here are required by Planck data. The sim¬ 
plicity of the power-law spectrum continues to give it an edge 
over more complicated initial spectra and the most plausible ex¬ 
planation for the apparent features in the data remains that we 
are just observing fluctuations due to cosmic variance at large 
scales and noise at small scales. 

10. Implications of P/anc^ bispectrai constraints 
on inflationary modeis 

The combination of power spectrum constraints and primor¬ 
dial non-Gaussianity (NG) constraints, such as the Planck upper 
bound on the NG amplitude /nl (Planck Collaboration XVII, 
2015), can be exploited to limit extensions to the simplest stan- 

* Due to the multimodal nature of the posterior, usual minimization 
routines perform poorly here. 


dat'd single-field models of slow-roll inflation. The next sub¬ 
section considers inflationary models with a non-standard ki¬ 
netic term (Garriga & Mukhanov, 1999), where the inflaton 
Lagrangian is a general function of the scalar inflaton field and 
its hrst derivative, £ = P((f>,X), where X = -g^’'5^0ciy(^/2 
(Garriga & Mukhanov, 1999; Chen et ah, 2007). Sect. 10.2 fo¬ 
cuses on a specihc example of a single-held model of inha- 
tion with more general higher-derivative operators, the so-called 
“Galileon inhation.” Sect. 10.3 presents constraints on axion 
monodromy inhation. See Planck Collaboration XVII (2015) for 
the analysis of other interesting non-standard inhationary mod¬ 
els, including warm inhation (Berera, 1995), whose /nl predic¬ 
tions can be constrained by Planck . 


10.1. Inflation with a non-standard kinetic term 

This class of models includes k-inhation (Atmendariz-Picon 
et ah, 1999; Garriga & Mukhanov, 1999) and Dirac-Born- 
Inheld (DBI) models introduced in the context of brane inha- 
tion (Silverstein & Tong, 2004; Alishahiha et ah, 2004; Chen, 
2005b,a). In these models inhation can take place despite a steep 
potential or may be driven by the kinetic term. 

Moreover, one of the main predictions of inhationary mod¬ 
els with a non-standard kinetic term is that the inhaton pertur¬ 
bations can propagate with a sound speed Cs < 1. We show 
how the Planck combined measurement of the power spectrum 
and the nonlinearity parameter /nl (Planck Collaboration XVII, 
2015) improves constraints on this class of models by breaking 
degeneracies between the parameters determining the observ¬ 
able power spectra. Such degeneracies (see, e.g., Peiris et al. 
(2007), Powell et al. (2009), Lorenz et al. (2008), Agarwal & 
Bean (2009), and Baumann et al. (2015)) are evident from the 
expressions for the power spectra. We adopt the same notation 
as Planck Collaboration XXIV (2014). At leading order in the 
slow-roll parameters the scalar power spectrum depends addi- 
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tionally on the sound speed Cs via (Garriga & Mukhanov, 1999) 




8;r2Mpj Csei ’ 


(84) 


which is evaluated at kcj = aH. Correspondingly, the scalar 
spectral index 

n^-\ - -2e\ - 62 - s (85) 

depends on an additional slow-roll parameter, s - CsKcsH), 
which describes the running of the sound speed. The usual con¬ 
sistency relation holding for the standard single-held models of 
slow-roll inhation (r = - 8 nt) is modihed to r ^ - 8 ntCs, with 
til - -2ei as usual (Garriga & Mukhanov, 1999).'° 

At lowest order in the slow-roll parameters, there are 
strong degeneracies between the parameters (As, Cs, ei, 62 , s)- 
This makes the constraints on these parameters from the power- 
spectrum alone not very stringent, and, for parameters like ei 
and 62 , less stringent compared with the standard case. However, 
combining the constraints on the power spectra observables with 
those on /nl can also result in a stringent test for this class 
of inhationary models. Models where the inhaton held has a 
non-standard kinetic term predict a high level of primordial NG 
of the scalar perturbations for Cs « 1, (see, e.g., Chen et ah, 
2007). Primordial NG is generated by the higher-derivative in¬ 
teraction terms arising from the expansion of the kinetic part of 
the Lagrangian, F(ip, X). There are two main contributions to the 
the amplitude of the NG (i.e., to the nonlinearity parameter /nl) 
coming from the inhaton held interaction terms 5/ (V(5/)^ and 
(5/)^ (Chen et al., 2007; Senatore et al., 2010). The NG from 
the hrst term scales as c~^, while the NG arising from the other 
term is determined by a second parameter C 3 (following the no¬ 
tation of Senatore et al. (2010)). Each of these two interactions 
produces bispectrum shapes similar to the so-called equilateral 
shape (Babich et al. 2004) for which the signal peaks for equi¬ 
lateral triangles with ki - k 2 - k^ (these two shapes are called, 
respectively, “EFTl” and “EET2” in Planck Collaboration XVII 
(2015)). However the difference between the two shapes is such 
that the total signal is a linear combination of the two, leading to 
an “orthogonal” bispectral template (Senatore et al., 2010). 

The equilateral and orthogonal NG amplitudes can be ex¬ 
pressed in terms of the two “microscopic” parameters Cs and C 3 
(for more details see Planck Collaboration XVII (2015)) accord¬ 
ing to 

1 - 

/nl'' = —^ [-0.275-0.0780c?-(2/3) X 0.780C 3 ], (87) 

1 - 

/°f° = —^ [ 0 . 0159 -0.0167c?-(2/3) X 0.0167 C 3 ] . (88) 


Thus the measurements of and obtained in the com¬ 
panion paper (Planck Collaboration XVII, 2015) provide a con¬ 
straint on the sound speed Cj of the inflaton field. Such con¬ 
straints allow us to combine the NG information with the anal¬ 
yses of the power spectra, since the sound speed is the non- 
Gaussianity parameter also affecting the power spectra. 

We use the more accurate relation 



Fig. 35. (ei, 62 ) 68 % and 95 % CL constraints for Planck data 
comparing the canonical Lagrangian case with Cs = 1 to the case 
of varying Cs with a uniform prior 0.024 < Cs < 1 derived from 
the Planck non-Gaussianity measurements. 


In this subsection we consider three cases. In the first case 
we perform a general analysis as described above (focusing on 
the simplest case of a constant sound speed, s - 0 ) improving 
on PCI 13 and Planck Collaboration XXIV (2014) by exploiting 
the full mission temperature and polarization data. The Planck 
constraints on primordial NG in general single-held models of 
inhation provide the most stringent bound on the inhaton sound 
speed (Planck Collaboration XVII, 2015) 

Cs> 0.024 (95% CL). (89) 

We then use this information on Cs as a uniform prior 0.024 < 
Cs < 1 in Eq. ( 86 ) within the HFF formalism, as in PCI13. Fig. 35 
shows the joint constraints on ei and 62 - Planck TT-tlowP yields 
61 < 0.031 at 95 % CL. No improvement in the upper bound on 
61 results when using Planck TT,TE,EEh-1owP. This constraint 
improves that of the previous analysis in PCI 13 and can be com¬ 
pared with the restricted case of Cs = 1, also shown in Fig. 35, 
with 6 \ < 0.0068 at 95 % CL. The limits on the sound speed from 
the constraints on primordial NG are crucial for deriving an up¬ 
per limit on ei, because the relation between the tensor-to-scalar 
ratio and 6 \ also involves the sound speed (see, e.g., Eq. ( 86 )). 
This breaks the degeneracy in the scalar spectral index. 

The other two cases analysed involve DBI models. The de¬ 
generacy between the different slow-roll parameters can be bro¬ 
ken for i = 0 or in the case where s oc 62 - We first consider 
models dehned by an action of the DBI form 

P(/, X) = -/(/)-' Vl - 2/(/)X + /(/)-' - V(/), (90) 

where V)/) is the potential and /(/) describes the warp factor 
determined by the geometry of the extra dimensions. We fol¬ 
low an analogous procedure to exploit the NG limits derived 
in Planck Collaboration XVII (2015) on Cs in the case of DBI 
models: Cj > 0.087 (at 95 % CL). Assuming a uniform prior 


r = , (86) 

accounting for different epochs of freeze-out for the scalar fluctuations 
(at sound horizon crossing kc^ = aH) and tensor perturbations (at hori¬ 
zon crossing k = aH) (Peiris et al., 2007; Powell et al., 2009; Lorenz 
et al., 2008; Agarwal & Bean, 2009; Baumann et al., 2015). 


" This section uses results based on /nl constraints from T + E 
(Planck Collaboration XVII, 2015). In Planck Collaboration XVII 
(2015) it is shown that although conservatively considered as prelim¬ 
inary, the /nl constraints from T + E are robust, since they pass a ex¬ 
tensive battery of validation tests and are in full agreement with T- only 
constraints. 
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0.087 < Cs < 1 and i = 0, Planck TT+lowP gives e\ < 0.024 at 
95 % CL, a 43 % improvement with respect to PCI13. The addi¬ 
tion of high-f” TE and EE does not improve the upper bound on 
ei for this DBI case. 

Next we update the constraints on the particularly interest¬ 
ing case of infrared DBI models (Chen 2005b, a), where /(0) a; 
Al<p‘^. (For details, see Silverstein & Tong (2004), Alishahiha 
et al. (2004), Chen et al. (2007) and references therein). In these 
models the inflaton field moves from the IR to the UV side with 
an inflaton potential 


V(0) = Vo - . (91) 

From a theoretical point of view a wide range of values for (3 
is allowed: 0.1 < P < 10^ (Bean et al., 2008). PCI13 dramati¬ 
cally restricted the allowed parameter space of these models in 
the limit where stringy effects can be neglected and the usual 
held theory computation of the primordial curvature perturba¬ 
tion holds (see Chen (2005a,c) and Bean et al. (2008) for more 
details). In this limit of the IR DBI model, one hnds (Chen, 
2005c; Chen et al., 2007) ^ (/lA,/3)“', - I - -4/A„ 

and dn^ldlnk - -4/A*. (In this model one can verify that 
s ss 1/At ~ ^ 2 / 3 ). Combining the uniform prior on Cs with 
Planck TTh-IowP, we obtain 

/?<0.31 (95% CL), (92) 


10.2. Galileon inflation 


As a further example of the implications of the NG con¬ 
straints on (non-standard) inhationary models we consider the 
case of Galileon inhation (Burrage et al., 2011) (see also 
Kobayashi et al. (2010), Mizuno & Koyama (2010), and Ohashi 
& Tsujikawa (2012)). This represents a well dehned and well 
motivated model of inhation with more general higher deriva¬ 
tives of the inhaton held compared to the non-standard ki¬ 
netic term case analyzed above. The Galileon models of inha- 
tion are based on the so-called “Galilean symmetry” (Nicolis 
et al., 2009), and enjoy some well understood stability prop¬ 
erties (absence of ghost instabilities and protection from large 
quantum corrections). This makes the theory also very predic¬ 
tive, since observable quantities (scalar and tensor power spec¬ 
tra and higher-order correlators) depend on a hnite number of 
parameters. From this point of view this class of models shares 
some of the same properties as the DBI inhationary models 
(Silverstein & Tong, 2004; Alishahiha et al., 2004). The Galileon 
held arises naturally within fundamental physics constructions 
(e.g., de Rham & Gabadadze 2010b, a). These models also offer 
an interesting example of large-scale modihcations to Einstein 
gravity. 

The Galileon model is based on the action (Deffayet et al., 
2009a,b) 


S^fd^x^Fg 


r, 3 3 

jK + JX, 

n=0 


(95) 


and a preference for a high number of e-folds: 78 < A* < 157 at 
95 % CL. 

We now constrain the general case of the IR DBI model, 
including the “stringy” regime, which occurs when the inhaton 
extends back in time towards the IR side (Bean et al., 2008). The 
stringy phase transition is characterized by an interesting phe¬ 
nomenology altering the predictions for cosmological perturba¬ 
tions. A parameterization of the power spectrum of curvature 
perturbations interpolating between the two regimes is (Bean 
et al. (2008); see also Ma et al. (2013)) 





1 

(1 H- x )2 


(93) 


where As = SIAjP I{n^ 0^) is the amplitude of the perturbations 
which depends on various microscopic parameters (hb is the 
number of branes at the B-throat; see Bean et al. (2008) for more 
details), while x = (N™^/Nc)^ sets the stringy phase transition 
taking place at the critical e-fold Ac. (Here is the number 
of e-folds to the end of IR DBI inhation). The spectral index and 
its running are 


^ 4 x^ H- 3x - 2 

~ (x H- l)(x H- 2) ’ 
dn^ 4 x"* H- 6x^ - 55x^ - 96x - 4 

dink ~ (Api^ (x H- l)2(x -H 2)2 


(94a) 

(94b) 


where 


£0 = czA, 

£i = -2(c3/A^)Xa,p, 


£2 = 2(c4/A®)A 
£2 = -2{cslK^)X 


(□<^)2 




•(c4/A®)a2r, 


(U(f,f - 'iUcf, (v^Vv0) + 2 (V/,Vv0)” 


H- 6(c5/A®)a2G/,vV'^V> . 


(96) 


Here A = -V^0V'^0/2, (V^Vv^)^ = and 

(V^Vy0)^ = V^Vv0V''V^(^V’'Vp0. The coupling coefficients c,- 
are dimensionless and A is the cut-off of the theory. The case of 
interest includes a potential term V{(p) - Vo+Acp+fl I2)m^(jp'+... 
to drive inhation. 

The predicted scalar power spectrum at leading order 
is (Ohashi & Tsujikawa (2012); Burrage et al. (201 1); Tsujikawa 
et al. (2013); see also Kobayashi et al. (2011a) and Gao & Steer 
(2011))'^ 










c^k=aH 


82 r 2 A(^o)^c^ ’ 


(98) 


where A = 1 H- c^icfuPpl(2Ep'M^p and = -BjA is the sound 
speed of the Galileon held. Cs is different from the usual slow- 
roll parameter ei and at leading order related according to - 
-2(Z?/1 + 6 c 3 + 18c4 + 30c5)ei. The scalar spectral index 


n.s - 1 = -2ei - ?7s - s (99) 


A prediction for the primordial NG in the stringy regime is not 
available. We assume the standard held-theoretic result for a pri¬ 
mordial bispectrum of the equilateral type with an amplitude 
= -(35/108) [(J3^ (A™')V9) - 1]. By considering the same 
uniform prior on Cj, we obtain p < 0.77, 66 < A™^ < 72, and 
X < 0.41 at 95 % CL, which severely limits the general IR DBI 
model and strongly restricts the allowed parameter space. 


For the following expressions it is convenient to dehne the quanti¬ 
ties 

A = C 2 / 2 - 1 - 6 C 3 + 27 c 4 - 1 - 6 OC 5 , S =-C 2 / 2 - 4 C 3 - 13c4-24^5, (97) 

where c,- = c,Z'“^ for i = 2 to 5, with Z = H^o/A^. In order to have 
a healthy model we require A > 0 (no ghosts) and B < 0 (no gradient 
instabilities). 
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depends on the slow-roll parameters e\, rj^ - esl(Hes), and 
s - Cs/(//Cs). As usual the slow-roll parameter s describes the 
running of the sound speed. In the following we restrict our¬ 
selves to the case of a constant sound speed with s = 0. The 
tensor-to-scalar ratio is 


r — IhejCs = 16eiCs, ( 100 ) 


where we have introduced the parameter Cs = -[ 2 B /(1 + 6 C 3 -H 
18c4 -h 30c5)]Cs, related to the Galileon sound speed. The param¬ 
eter Cs can be either positive or negative. In the negative branch 
a blue spectral tilt for the primordial gravitational waves is al¬ 
lowed, contrary to the situation for standard slow-roll models of 
inflation. We introduce such a quantity so that the consistency 
relation takes the form r « -SutCs, with Ut = - 2 ei analogous 
to Eq. ( 86 ). The measurements of primordial NG constrain Cs, 
which in turn constrain e\ and t/s in Eq. (99). This is analogous 
to the constraints on ei and rf of Eq. (85) in the previous subsec¬ 
tion. 

Galileon models of inflation predict interesting NG sig¬ 
natures (Burrage et al., 2011; Tsujikawa et ah, 2013).'^ We 
have verified (see also Creminelli et al. (2011)) that bispec¬ 
tra can be generated with the same shapes as the “EETl” and 
“EET2” (Senatore et al. (2010); Chen et al. (2007)) constrained 
in the companion paper (Planck Collaboration XVII, 2015), 
which usually arise in models of inflation with non-standard ki¬ 
netic terms, with 


^EFTl 

17 

/nl ” 

VT2 

^EFr2 

1 

/nl ” 



5 30 

Cs Cs 

30/A 




40 

CsCs 

-55 


-H 15 
40 

CsCs 


Cs 30 

320^^- +215 

Cs A 


225cf -H 280- 




( 101 ) 


As explained in the previous subsection, the linear combinations 
of these two bispectra produce both equilateral and orthogonal 
bispectrum templates. Given Eqs. (98)-(101), we can proceed as 
in the previous section to exploit the limits on primordial NG in 
a combined analysis with the power spectra analysis. In Planck 
Collaboration XVII (2015) the constraint Cs > 0.23 (95 % CL) 
is obtained based on the constraints on and One can 
proceed as described in Planck Collaboration XVII (2015) to 
constrain the parameter Cs modifying the consistency relation 
Eq. (100). Adopting a log-uniform prior on A in the interval 
10 < A < 10"^ and a uniform prior 10 < Cs < 1, the Planck 

measurements on and constrain Cs to be 0.038 <Cs< 
100 (95 % CL) (Planck Collaboration XVII, 2015). We also ex¬ 
plore the possibility of the negative branch (corresponding to a 
blue tensor spectral index), finding -100 < q 5 -0.034 (95 % 
CL) (Planck Collaboration XVII, 2015). By allowing a logarith¬ 
mic prior on Cs based on the /nl measurements. Pig. 36 shows 
the joint constraints on ei and for the Ut < 0 branch and for the 
Ut > 0 branch. Planck TTh-IowPh-BAO and the NG bounds on c^ 
constrain ei < 0.036 at 95 % CL for Ut < 0 (and |ei| < 0.041 for 

Ut > 0). 


10.3. Axiom monodromy inflation 

The mechanism of monodromy inflation (Silverstein & 
Westphal, 2008; McAllister et al., 2010; Kaloper et al., 2011; 

See also Mizuno & Koyama (2010), Gao & Steer (2011), 
Kobayashi et al. (2011b), De Felice & Tsujikawa (2013), and Regan 
etal. (2014). 


Planck TT + lowP + BAO 




Fig. 36. Marginalized joint 68 % and 95 % CL for the Galileon 
parameters (ei , rjs) for n^ < 0 (left panel) and Ut > 0 (right 
panel). 


Plauger et al., 2014b) in string theory motivates a broad class of 
inflationary potentials of the form 


V(4>) ^ + A^e 




COS 


^ /_0 V 

/o \ 0 o) 


Pf+l 


. ( 102 ) 


Here ju, Aq, /o, and are constants with the dimension of mass 
and Co, p, p\, Pf, and 70 are dimensionless. 

In simpler parameterizations used in prior analyses of os¬ 
cillations from axion monodromy inflation (Peiris et al., 2013; 
Planck Collaboration XXII, 2014; Easther & Plauger, 2014; 
Jackson et al., 2014; Meerburg et al., 2014b; Meerburg & 
Spergel, 2014; Meerburg et al., 2014a; Meerburg, 2014), one 
assumes p\ - pj - 0, corresponding to a sinusoidal term 
with constant amplitude throughout inflation, taken to be a pe¬ 
riodic function of the canonically-normalized inflaton (p. Taking 
Pa ^ 0 and pj + Q allows the magnitude dxid frequency, re¬ 
spectively, of the modulation to depend on f. Lor example, the 
frequency is always a periodic function of an underlying angular 
axion held, but its relation to the canonically normalized inflaton 
field is model dependent. 

The microphysical motivation for pA + 0 and pf + 0 is 
that in string theory additional scalar fields, known as “moduli,” 
evolve during inflation. The inflationary potential depends on 
a subset of these fields. Because the magnitude and frequency 
of modulations are determined by the vacuum expectation val¬ 
ues of moduli, both quantities are then naturally functions of f. 
The case pA - Pf = 0 corresponds to when these fields are 
approximately fixed, stabilized strongly by additional terms in 
the scalar potential. But in other cases, the axion potential that 
drives inflation also provides a leading term stabilizing the mod¬ 
uli. The exponential dependence of the magnitude in the poten¬ 
tial of Eq. (102) arises because the modulations are generated 
non-perturbatively, e.g., by instantons. Por this reason, the mod¬ 
ulations can be undetectably small in this framework, although 
there are interesting regimes where they could be visible. 

Specific examples studied thus far yield exponents p, pA, 
and Pf which are rational numbers of modest size. Por example, 
models with p - 3, 2, 4/3, 1, and 2/3 have been constructed 
(Silverstein & Westphal, 2008; McAllister et al., 2010, 2014) or 
in another case p = 4/3, pa = -1/3, and pf - -1/3). Pollowing 
Plauger et al. (2014b), we investigate the effect of a drift in fre¬ 
quency arising from pf, neglecting a possible drift in the modu¬ 
lation amplitude by setting pa — Co - 0. Even in this restricted 
model, a parameter exploration using a fully numerical compu¬ 
tation of the primordial power spectrum following the method¬ 
ology of Peiris et al. (2013) is prohibitive, so we follow Plauger 
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et al. (2014b) studying two templates capturing the features of 
the primordial spectra generated by this potential. 

The first template, which we term the “semi-analytic” tem¬ 
plate, is given by 


/ k 


1 + 6ns cos 


/ 



(103) 

The parameter / is higher than /o, the underlying axion de¬ 
cay constant in the potential, by a few percent, but this dif¬ 
ference will be neglected in this analysis. The quantity (pQ is 
some fiducial value for the scalar field, and (pk is the value of 
the scalar field at the time when the mode with comoving mo¬ 
mentum k exits the horizon. At leading order in the slow-roll 
expansion, in units where the reduced Planck mass Mpi = 1, 


cpk = ^j2p {No - MkjK)), where No^N^+ and 0end 

is the value of the scalar field at the end of inflation. 

The second “analytic” template was derived by Flauger et al. 
(2014b) by expanding the argument of the trigonometric func¬ 
tion in Eq. (103) in ln(A:/kt), leading to 


Vnik) = VniK) 



X 


|l H- 6ns cos 


A0 H- a 


In 



N'l 


Ih 


n+\ 



(104) 


The relation between the empirical parameters in the tem¬ 
plates and the potential parameters are approximated by 6ns - 
3b yl{2n)la, where 
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and b is the monotonicity parameter defined in Flauger et al. 
(2014b), providing relations converting bounds on c„ into 
bounds on the microphysical parameters of the potential. 
However, the analytic template can describe more general shapes 
of primordial spectra than just axion monodromy. 

As discussed by Flauger et al. (2014b), there is a degener¬ 
acy between p (or alternatively ns) and /, and for both templates 
we fix p = 4/3 as well as fixing the tensor power spectrum to 
its form in the absence of oscillations. This is an excellent ap¬ 
proximation because tensor oscillations are suppressed relative 
to the scalar oscillations by a factor a(//Mpi)^ « 1. A uniform 
prior —n < A/ < tt is adopted for the phase parameter of both 
templates as well as a prior 0 < 6ns < 0.7 for the modulation 
amplitude parameter. 

In order to specify the semi-analytic template we assume 
instantaneous reheating, which for p - A13 corresponds to 
M* a! 57.5 for k* = 0.05Mpc^'. We set /o = 12.38 Mpi with 
0end = 0.59 Mpi. We adopt uniform priors -4 < log[o(//Mpi) < 
-1 and -0.75 < pj < \ for the remaining parameters. The priors 
0 < ln(a) < 6.9 and -2 < 0^2 < 2 specify the analytic template. 
The single-field effective field theory becomes strongly coupled 
for a > 200. However, in principle the string construction re¬ 
mains valid in this regime. 


10.3.1. Power spectrum constraints on monodromy 
inflation 

We carry out a Bayesian analysis of axion monodromy infla¬ 
tion using a high-resolution version of CAMB coupled to the 


PolyChord sampler (see Sect. 8.2). For our baseline analy¬ 
sis we conservatively adopt the PLIK T-only “binl” likelihood, 
using only low-f polarization data. In addition to the primordial 
template priors specified above, we marginalize over the stan¬ 
dard priors for the cosmological parameters, the primordial am¬ 
plitude, and foreground parameters. 

The marginalized joint posterior constraints on pairs of pri¬ 
mordial parameters for the semi-analytic and analytic templates 
are shown in Figs. 38 and 37, respectively. 

The complex structures seen in these plots arise due to de¬ 
generacies in the likelihood frequency “beating” between un¬ 
derlying modulations in the data and the model (Easther et al., 
2005). Parameter combinations where “beating” occurs over the 
largest k ranges lead to discrete local maxima in the likelihood. 
Fortuitous correlations in the observed realization of the Ce can 
give the same effect. 

The four frequencies picked out by these structures, 
In(Q') {3.5,5.4,6.0,6.81, show improvements of ~ 

(-9.7,-7.1,-12.2,-12.5) relative to ACDM, respectively. 
These frequencies are marked by dotted lines in Fig. 37, and by 
solid lines in Fig. 38 using Eq. (105). The semi-analytic and an¬ 
alytic templates lead to self-consistent results as expected, with 
analogous structures being picked out by the likelihood in each 
template. There is no evidence for a drifting frequency, pf ^ 0 
or c„ 0. Thus, these parameters serve to smooth out structures 
in the marginalized posterior. 

The improvement in is not compelling enough to sug¬ 
gest a primordial origin. Fitting a modulated model to simu¬ 
lations with a smooth spectrum can give rise to ^xlg ~ -10 
improvements (Flauger et al., 2014b). Furthermore, as the mon¬ 
odromy model contains only a single frequency, at least three 
of these structures must correspond to spurious fits to the noise. 
Considering the two models defined by the two templates and the 
parameter priors specified above, the Bayes factors calculated 
using PolyChord favours base ACDM over both templates by 
odds of roughly 8:1. 

Compared to previous analyses of the linear {p - 1) axion 
monodromy model for WMAP9 (Peiris et al., 2013) and the 
2013 Planck data (Planck Collaboration XXII, 2014; Easther 
& Flauger, 2014) the common frequencies are shifted slightly 
higher. The lower frequency in common appears shifted by a fac¬ 
tor of order yfp, from a =» 28.9 to 31.8, and the higher frequency 
in common from a ^ 210 to 223. Flauger et al. (2014b) suggest 
that the lower frequency (which had A^^^ = -9 in PCI13) was 
associated with the 4-K cooler line systematic effects in the 2013 
Planck likelihood. However its presence at similar significance 
in the 2015 likelihood with improved handling of the cooler line 
systematics suggests that this explanation is not correct. The sec¬ 
ond frequency, which appeared with Axlg ~ -20 in WMAP9 
(Peiris et al., 2013) is still present but with much reduced sig¬ 
nificance, suggesting that the high multipoles do not give evi¬ 
dence for this frequency. Additionally, two higher frequencies 
are present, which, if interpreted as being of primordial origin, 
correspond to a regime well beyond the validity of the single- 
field effective field theory. If one of these frequencies were to be 
confirmed as primordial, a significantly improved understanding 
of the underlying string construction would need to be under¬ 
taken. 

In order to check whether the improvement in fit at these 
four modulation frequencies is responding to residual fore¬ 
grounds or other systematics, we examine the frequency residu¬ 
als. Figure 39 shows the residuals of the data minus the model 
(including the best-fit foreground model) for the four PL IK fre¬ 
quency combinations binned at Af = 30 for the lowest mod- 
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Fig. 37. Constraints on the parameters of the analytic template, showing joint 68 % and 95 % CL. The dotted lines correspond to the 
frequencies showing the highest likelihood improvements (see text). 



Fig. 38. Constraints on the parameters of the semi-analytic template, showing joint 68 % and 95 % CL. The solid lines on the 
left-hand panel mark the frequencies showing the highest likelihood improvements (see text). 


ulation frequency, ln(a) 3.5. This plot shows no significant 
frequency dependence, and thus there is no indication that the fit 
is responding to frequency-dependent systematics. Furthermore, 
the plot does not show evidence that the improvement for this 
modulation frequency comes from the feature at ^ « 800, as 
suggested by Easther & Flauger (2014). This feature and another 
21 { X 1500 are apparent at all frequency combinations. Similar 
plots for the three other modulation frequencies also do not show 
indications of frequency dependence. 

In order to confirm whether any of the frequencies picked 
out here have a primordial origin, one can make use of inde¬ 
pendent information in the polarization data to perform a cross¬ 
check of the temperature prediction which minimizes the “look- 
elsewhere” effect (Mortonson et al., 2009). Leaving a complete 
analysis of the independent information in the polarization for 
future work, we now check whether the temperature-only result 
remains stable when high-^ polarization is added in the likeli¬ 
hood. In Fig. 40 we show a preliminary analysis using the PLIK 
temperature and polarization “binl” likelihood plus low-^ po¬ 
larization data. Comparing with the left-hand panels of Figs. 37 
and 38, there are slight differences from the T-only analysis. 
However, all the four frequencies identified in the temperature 
are present when high-f polarization is added. 

10.3.2. Predictions for resonant non-Gaussianity 

The left-hand panel of Fig. 41 presents derived constraints on 
the parameters of the potential, Eq. (102), calculated using the 
analytic template. Another cross-check of primordial origin is 
available since the monodromy model predicts resonant non- 


Gaussianity, generating a bispectrum whose properties would be 
strongly cortelated with that of the power spectrum (Chen et al., 
2008; Flauger & Pajer, 2011). Using the mapping 

/nT = ( 106 ) 

we use the analytic template to derive the posterior probability 
for the resonant non-Gaussianity signal predicted by constraints 
from the power spectram, presented in the middle and right pan¬ 
els of Fig. 41. 

Planck Collaboration XVII (2015) use an improved modal 
estimator to scan for resonant non-Gaussianity. The resolution 
of this scan is curi'ently limited to In(Q') < 3.9, which is poten¬ 
tially able to probe the lowest frequency picked out in the power 
spectrum search. However, the modal estimator’s sensitivity (im¬ 
posed by cosmic variance) of a; 80 is significantly greater 
than the predicted value for this frequency from fits to the power 
spectrum, ~ 10. Efforts to increase the resolution of the 
modal estimator are ongoing and may allow consistency tests of 
the significantly higher levels of resonant non-Gaussianity pre¬ 
dicted by the higher frequencies in the future. 

10.3.3. Power spectrum and bispectrum constraints on 
axion inflation with a gauge fieid coupiing 

We now consider the case where the axion field is coupled to a 
gauge field. Such a scenario is physically well motivated. From 
an effective field theory point of view the derivative coupling is 
natural and must be included since it respects the same shift sym¬ 
metry that leads to axion models of inflation (Anber & Sorbo, 
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Fig. 39. Frequency residuals for the In(Q') ^3.5 likelihood peak, 
binned at A£ - 30. The +1 cr errors are given by the square root 
of the diagonal elements of the covariance matrix. 


2010; Barnaby & Peloso, 2011; Pajer & Peloso, 2013). This 
type of coupling is also ubiquitous in string theory (see, e.g., 
Barnaby et al. (2012) and Linde et al. (2013)). The coupling term 
in the action is (Anber & Sorbo, 2010; Barnaby & Peloso, 2011; 
Barnaby et al., 201 1) 

S D (107) 


where = d/jAy - dyA^, its dual is = e^''°‘^Fc,pl2, and a 
is a dimensionless constant which, from an effective field theory 
perspective, is expected to be of order one. For the potential of 
the axion field, we will not investigate further the consequences 
of the oscillatory part of the potential, focusing on the coupling 
of the axion field to the U(l) gauge field (effectively setting Aq = 
0). 

The coupling of a pseudo-scalar axion with the gauge field 
has interesting phenomenological consequences, both for den¬ 
sity perturbations and primordial gravitational waves (Barnaby 
& Peloso, 2011; Sorbo, 2011; Barnaby et al., 2011, 2012; 
Meerburg & Pajer, 2013; Ferreira & Sloth, 2014). Gauge field 
quanta source the axion field via an inverse decay process 
6A + 6A —> dtp, modifying the usual predictions already at the 
power-spectrum level. Additionally, the inverse decay can gen¬ 
erate a high level of primordial NG. 

The parameter 
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Fig. 40. Constraints on the parameters of the analytic (top) and 
semi-analytic (bottom) templates with the addition of high- 
{ polarization data in the likelihood, showing joint 68 % and 
95 % CL. The lines mark the frequencies showing the highest- 
likelihood improvements identified in the baseline temperature- 
only analysis. 



characterizes the strength of the inverse decay effects. If ^ < 1 
the coupling is too small to produce any modifications to the 
usual predictions of the uncoupled model. For previous con¬ 
straints on ^ see Barnaby et al. (2011, 2012) and Meerburg & 
Pajer (2013). Using the slow-roll approximation and neglecting 
the small oscillatory part of the potential. One can express 


where N is the number of e-folds to the end of inflation, defined 
by eiend = 1 ■ The scalar power spectrum of curvature perturba¬ 
tion is given by 








(108) 


( 110 ) 
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Fig. 41. Derived constraints on the parameters of the potential, Eq. (102), as well as the predicted resonant non-Gaussianity, /j^“, 
using the analytic template, showing joint 68 % and 95 % CL. The dotted lines mark the frequencies showing the highest likelihood 
improvements (see text). 


where (Meerburg & Pajer, 2013) 


m = 



(111) 


Here an asterisk indicates evaluation at the pivot scale k, - 
0.05 Mpc“^ and = H^/(4;r^(^^) and - 1 = -2e« - 77 , 
are the amplitude and spectral index of the standard slow-roll 
power spectrum of vacuum-mode curvature perturbations (the 
usual power spectrum in the absence of the gauge-coupling). 
Numerically evaluating the function / 2 (f) (defined in equa¬ 
tion (3.27) of Barnaby et al. (2011)), we created an analytical 
fit to this function, which is accurate to better than 2 % in the 
range 0.1 < ^ < 7.'"^ In the following, unless stated otherwise, 
we fix p = 4/3 as in the previous subsection and assume instan¬ 
taneous reheating so that N* a; 57.5 and the slow-roll parameters 
e» and 77 * are fixed. For the tensor power spectrum we adopt the 
approximation (Barnaby et al., 201 1) 


The function corresponds to the quantity 1,1) de¬ 
fined in equation (3.29) of Barnaby et al. (201 1). We have com¬ 
puted / 3 (^t) numerically and employed a fit with an accuracy of 
better than 2 We use the observational constraint = 

22.7 + 25.5 ( 68 % CL) obtained in Planck Collaboration XVII 
(2015) from an analysis where only the inverse decay type non- 
Gaussianity is assumed present, so we omit the explicit expres¬ 
sion for the resonant bispectrum Bjgs. since it will not be used 
here. 

We carried out an MCMC analysis of constraints on the 
(scalar and tensor) power spectra predicted by this model with 
the standard PL IK temperature-only likelihood plus low-f po¬ 
larization, marginalizing over standard priors for the cosmolog¬ 
ical parameters and foreground parameters, with uniform priors 
2.5 < ^[10'°!^*] < 3.7 and 0.1 < ^. < 7.0. 

The power spectrum constraint gives 

0.1 <^*<2.3 (95% CL). (116) 


r,(k) = n 



1 + y^t/t.L(f(^))e""^- - 




, ( 112 ) 


where 

Ai(«k)) = 2.6-10-V-’(k). (113) 

Here Vi - and rii - -le, are the “usual” expres¬ 

sions for the tensor amplitude and tensor tilt in standard slow-roll 
inflation. 

The total bispectrum is (Barnaby et al., 2012) 

B{ki) = Bi„v.dec.(^;) + Bres(^i) (114a) 

= /nL^ inv.dec.(^,) + Bres(^,), (114b) 


where the explicit expression for Finv.dec.)^!) (Barnaby et al. 
(201 1); see also Meerburg &. Pajer (2013)) is reported in Planck 
Collaboration XVII (2015). This shows that the inverse decay 
effects and the resonant effects (which arise from the oscillatory 
part of the potential) simply “add up” in the bispectrum. The 
nonlinearity parameter is 


^inv.dec. 

/nl 




(115) 


Given that is exponentially sensitive to this translates 

into the prediction (using Eq. (115)) < 1.2, which is 

significantly tighter than the current bispectrum constraint from 
Planck Collaboration XVII (2015). Indeed, importance sampling 
with the likelihood for , taken to be a Gaussian centred 

on the NG estimate = 22.7 + 25.5 ( 68 % CL) (Planck 

Collaboration XVII, 2015), changes the limit on only at the 
second decimal place. 

We now derive constraints on model parameters using only 
the observational constraint on The constraints thus de¬ 

rived are applicable for generic p and also to the axion mon- 
odromy model discussed in Sect. 10.3, even in the case Aq 3 ^ 0. 
We follow the procedure described in section 11 of Planck 
Collaboration XVII (2015). The likelihood foris taken to 
be a Gaussian centred on the NG estimate = 22.7 ± 25.5 

( 68 % CL) (Planck Collaboration XVII, 2015). We use the ex¬ 
pression of Eq. (115), where is numerically evaluated. To 
find the posterior distribution for the parameter we choose 
uniform priors in the intervals 1.5 x 10 ® < 3.0 x 10“® and 

0.1 < < 7.0. This yields 95 % CL constraints for (for any 

value of p) of 

<2.5 (95% CL). (117) 


The fitting function used is exp{-fl-Zi ln(^)-c [ln(f)]^ -t-d [ln(f)]^ -l- 
e[ln(^)]"'), where the coefficients are a = 10.8, b - 4.58, c = 0.51, 
d = 0 .01, and e = 0 .02. 


The fit has the same expression as the one for / 2 (f) with coefficients 
a = 17.0048, b = 6.6578, c = 0.96479, d = 0.0506098, and e = 
0.039139. 
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If we choose a log-constant prior on we find 

^.<2.2 (95% CL). (118) 

For both cases the results are insensitive to the upper limit cho¬ 
sen for the prior on since the likelihood quickly goes to zero 
for > 3. As the likelihood for is fairly flat, the tighter con¬ 
straint seen for the log-constant case is mildly prior-driven. The 
constraints from the bispectrum are consistent with and slightly 
worse than the result from the power spectrum alone. 

Using a similar procedure and Eq. (109) one can also obtain 
a constraint on a//. Adopting a log-constant prior 2 < ajf < 
100'® and uniform priors 50 < N, < 70 and 1.5 x 10“^ < < 

3.0 X 10“® we obtain the 95 % CL constraints 


a// < 48 Mpj' for p = 1, 

ajf < 35 Mp,' 

for p -2, 



(119) 

and 



cr//<42Mpi' 

for p = 4/3. 

(120) 


For example, for a linear potential, p = 1, if a = (9(1) as sug¬ 
gested by effective held theory, then the axion decay constant / 
is constrained to be 


/ > 0.020 Mpi (95 % CL), (121) 

while for a potential with p - 413 

/> 0.023 Mpi (95% CL). (122) 

These limits are complementary to those derived in Sect. 10.3 
where a gauge coupling of the axion held is taken into account. 

11. Constraints on isocurvature modes 

In PCI13, we presented constraints on a number of simple mod¬ 
els featuring a mixture of the adiabatic (ADI) mode and one 
type of isocurvature mode. We covered the cases of CDM den¬ 
sity isocurvature (CDI), neutrino density isocurvature (NDI), 
and neutrino velocity isocurvature (NVI) modes (Bucher et ah, 
2000), with different assumptions concerning the correlation 
(Langlois, 1999; Amendola et ah, 2002) between the primordial 
adiabatic and isocurvature perturbations. Isocurvature modes, 
possibly correlated among themselves and with the adiabatic 
mode, can be generated in multi-held models of inhation; how¬ 
ever, at present a mechanism for exciting the neutrino velocity 
isocurvature mode is lacking. Section 11.2 shows how these con¬ 
straints have evolved with the new Planck TT-tlowP likelihoods, 
how much including the Planck lensing likelihood changes the 
results, and what extra information the Planck high-^ polariza¬ 
tion contributes. A pure isocurvature mode as a sole source of 
perturbations has been ruled out (Enqvist et ah, 2002), since, as 
can be seen from Fig. 42, any of the isocurvature modes leads 
to an acoustic peak structure for the temperature angular power 
very different from the adiabatic mode, which hts the data very 
well. The different phases and tilts of the various modes also 
occur in the polarization spectra, as shown in Fig. 42 for the E 
mode.*^ 

We give only the results for a log-constant prior on a/f, which is 
well-motivated since it corresponds to a log-constant prior on the axion 
decay constant for some fixed a. 

The transfer function mapping the primordial CDI mode to is 
suppressed by a factor (kjk^y^ ~ (f/Aq)“^ relative to the ADI mode, 
where keq is the wavenumber of matter-radiation equality. As seen in 
Fig. 42, there is a similar damping for the E mode in the CDI versus 



Fig. 42. Angular power spectra for the scale invariant (i.e., n^K = 
1) pure adiabatic mode (ADI, green dashed curves) and for the 
scale invariant (njj =1) pure isocurvature (CDI, NDI, or NVI) 
modes, with equal primordial perturbation amplitudes. The thick 
lines represent the temperature auto-correlation (TT), and the 
thin lines the E mode polarization auto-correlation (EE). 

In Sect. 11.4 we add one extra degree of freedom to the 
generally-correlated ADIh-CDI model by allowing primordial 
tensor perturbations (assuming the inflationary consistency re¬ 
lation for the tilt of the tensor power spectrum and its running). 
Our main goal is to explore a possible degeneracy between ten¬ 
sor modes and negatively correlated CDI modes, tending to tilt 
the large-scale temperature spectrum in opposite directions. In 
Sect. 11.5, we update the constraints on three special cases mo¬ 
tivated by axion or curvaton scenarios. 

The goal of this analysis is to test the hypothesis of adi- 
abaticity, and establish the robustness of the base ACDM 
model against different assumptions concerning initial condi¬ 
tions (Sect. 11.3). Adiabaticity is also an important probe of the 
inflationary paradigm, since any significant detection of isocur¬ 
vature modes would exclude the possibility that all perturbations 
in the Universe emerged from quantum fluctuations of a single 
inflaton field, which can excite only one degree of freedom, the 
curvature (i.e., adiabatic) perturbation.** 

the ADI case. Therefore, to be observable at high f, a CDI mode should 
be (highly) blue tilted. So, if the data favoured as small as possible a 
disturbance by CDI over all scales, then the CDI should have a spectral 
index, ujj, of roughly three. In practice, the lowest-f part of the data 
has very little weight due to cosmic variance, and thus we expect that 
the data should favour njj less than three, but significantly larger than 
one. This should be kept in mind when interpreting the results in the 
CDI case, i.e., one cannot expect strong constraints on the primordial 
CDI fraction at small scales, even if the data are purely adiabatic. The 
imprint of the baryon density isocurvature (BDI) mode, at least at linear 
order, in the CMB is indistinguishable from the CDI case, and hence 
we do not consider it separately as it can be described by = 

JcDi + (Db/D(.)JBDi- The trispectrum, however, can in principle be used 
to distinguish the BDI and CDI modes (Grin et ah, 2014). 

However, conversely, if no isocurvature was detected, the fluctua¬ 
tions could have been seeded either by single- or multi-field inflation, 
since later processes easily wash out inflationary isocurvature pertur- 
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In this section, theoretical predictions were obtained with 
a modified version of the camb code (version Jull4) while 
parameter exploration was performed with the MultiNest 
nested sampling algorithm. 


11.1. Parameterization and notation 


A general mixture of the adiabatic mode and one isocurva¬ 
ture mode is described by the three functions 'PuKik), PiiQr), 
and 'Pniik), describing the curvature, isocurvature, and cross¬ 
correlation power spectra, respectively. Our sign conventions are 
such that positive values for Pki correspond to a positive contri¬ 
bution of the cross-correlation term to the Sachs-Wolfe compo¬ 
nent of the total temperature spectrum. 

As in PCI 13, we specify the amplitudes at two scales k\ < k 2 
and assume power-law behaviour, so that 


Poiik) = exp 


/ ln(^)-ln(^2) \ / an 

iln(^i)-ln(^2)/ 
/ ln(^)-ln(^,) \ . pn 

\ln(k 2 )-ln(ki)i V 


(123) 


where a,b = I and I - Icm, -^ndu or -^nvi- We set 
k\ — 0.002Mpc^' and k 2 — O.lOOMpc^', so that [ki,k 2 \ spans 
most of the range constrained by the Planck data. The positive 
definiteness of the initial condition matrix imposes a constraint 
on its elements at any value of k: 

VPAk)f <r^{k)n-,{k). ( 124 ) 


We take uniform priors on the positive amplitudes, 

e(10-^ 10-^), (125) 

^ (0, 10^®). (126) 

The correlation spectrum can be positive or negative. For a 4^ b 
we apply a uniform prior at large scales (at ki). 


6 (-10^^ 10^^), (127) 


but reject all parameter combinations violating the constraint in 
Eq. (124). To ensure that Eq. (124) holds for all k, we restrict 
ourselves to a scale-independent correlation fraction 

cosAab^ e(-l,l). (128) 

(naPbb)'^" 

Thus is a derived parameter^® given by 


/p(2) 

^ab 


,(1)V 


(129) 


bations (Mollerach, 1990; Weinberg, 2004; Beltran et al., 2005). An 
example is the curvaton model, in which perturbations can be purely 
isocurvature at Hubble exit during inflation, but are later converted to 
ADI if the curvaton or curvaton particles (Linde & Mukhanov, 2006) 
dominate the energy density at the curvaton’s decay. For a summary of 
various curvaton scenarios, see, e.g., Gordon & Lewis (2003). 

Given our ansatz of power-law primordial spectra, if we treated P® 
as an independent parameter as we do with Eq. (124) would always 
be violated somewhere outside [k \, k^}. In PCI13, we dealt with this by 
assuming that when maximal (anti-)correlation is reached at some scale, 
the correlation remains at (-)IOO % beyond this scale. This introduced a 
kink in the cross-correlation spectrum, located at a different wavenum¬ 
ber for each model. Even though the range [ki,k 2 \ was chosen to span 
most of the observable scales, this kink tended to impact the smallest 
(or largest) multipole values used in the analysis. In particular, the kink 
helped fit the dip in the temperature angular power in the multipole 
range £ x 10^0. 


which, in terms of spectral indices, is equivalent to 

«ab = ^(«aa + Wbb) ■ (130) 

The conservative baseline likelihood is Planck TTh-IowP. 
The results obtained with Planck TT,TE,EEh-1owP should be in¬ 
terpreted with caution because the data used in the 2015 release 
are known to contain some low level systematics, in particular 
arising from ?’—>£’ leakage, and it is possible that such system¬ 
atics may be fit by the isocurvature auto-correlation and cross¬ 
correlation templates. (See Planck Collaboration XIII (2015) for 
a detailed discussion.) 

In what follows, we quote our results in terms of derived 
parameters identical to those in PCI13. We define the primordial 
isocurvature fraction as 


Aso(^) = 


Piiik) 

Pmik) + Piiik) ■ 


(131) 


Unlike the primordial correlation fraction cos A defined in 
Eq. (128), ySiso is scale-dependent in the general case. We present 
bounds on this quantity at kio* = ki, kmid = 0.050 Mpc^*, and 

^high ~ ^2- 

We report constraints on the relative adiabatic (ab = 'R'R), 
isocurvature (ab = II), and correlation (ab = 'RI) according to 
their contribution to the observed CMB temperature variance in 
various multipole ranges: 


^ab(^min» ^max) “ 


(A7’)^j^(f’niin, fmax) 

(A7’)2^((f’ndn,Cax)’ 


(132) 


where 

P 

t-max 

(Ar)2,(Ci„, Cax) = X 

The ranges considered are (^min,^max) = (2,20), (21,200), 
(201,2500), and (2,2500), where the last range describes the 
total contribution to the observed CMB temperature variance. 
Here okr measures the adiabaticity of the temperature angular 
power spectrum, a value of unity meaning “fully adiabatic initial 
conditions.” Values less than unity mean that some of the ob¬ 
served power comes from the isocurvature or correlation spec¬ 
trum, while values larger than unity mean that some of the power 
is “cancelled” by a negatively correlated isocurvature contribu¬ 
tion. The relative non-adiabatic contribution can be expressed as 

Q^non-adi = 1 “ Orr - ajl + Uri . 


11.2. Results for generaiiy-correiated adiabatic and one 
isocurvature mode (CDI, NDI, or N VI) 

Results are reported as 2D and ID marginalized posterior prob¬ 
ability distributions. Numerical 95% CL intervals or upper 
bounds are tabulated in Table 15. 

Eigure 43 shows the Planck 68 % and 95 % CL contours for 
various 2D combinations of the primordial adiabatic and isocur¬ 
vature amplitude parameters at large scales (k^ - 0.002 Mpc^') 
and small scales (^2 = O.lOOMpc^*), for (a) the generally- 
correlated ADIh-CDI, (b) ADIh-NDI, and (c) ADIh-NVI mod¬ 
els. Overall, the results using Planck TT-tlowP are consistent 
with the nominal mission results in PCI13, but slightly tighter. 
In the first panels of Eigs. 43 (a), (b), and (c), we also show 
the constraints on the curvature perturbation power in the pure 
adiabatic case. Comparing the generally-correlated isocurvature 
case to the pure adiabatic case with the same data combination 
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Fig. 43. 68 % and 95 % CL constraints on the primordial perturbation power in general mixed ADI+CDl (a), ADl+NDI (b), and 
ADI+NVl (c) models at two scales, ki = 0.002 Mpc“^ (l)and^2 = O.lOOMpc^* (2), for flancA: TT+lowP (grey regions highlighted 
by dotted contours), Planck TT+lowP+lensing (blue), and Planck TT,TE,EE+lowP (red). In the first panels, we also show contours 
for the pure adiabatic base ACDM model, with the corresponding colours of solid lines. 
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Fig. 44. Constraints on the primordial isocurvature fraction, ySjso, at = 0.002 Mpc^' and khigh = O.lOOMpc^', the primordial 
correlation fraction, cos A, the adiabatic spectral index, n<Rji, the isocurvature spectral index, tijj, and the correlation spectral index, 
nKi - (nm + njj)/2, with Planck TT+lowP data (dashed curves) and TT,TE,EE+lowP data (solid curves), for the generally- 
correlated mixed ADI+CDI (black), ADI+NDI (red), and ADI+NVI (blue) models. All these parameters are derived, and the distri¬ 
butions shown here result from a uniform prior on the primary parameters, as detailed in Eqs. (125)-(127). However, the effect of 
the non-flat derived-parameter priors is negligible for all parameters except for njj (and nnj) where the prior biases the distribution 
toward one. With TT+lowP, the flatness of jSiso(^high) in the CDI case up to a “threshold” value of about 0.5 is a consequence of the 
(^/^eq) ^ damping of its transfer function as explained in Eootnote 17. 


summarizes neatly what the data tell us about the initial con¬ 
ditions. If the contours in the -plane were shifted sig¬ 

nificantly relative to the pure adiabatic case, the missing power 
could come either from the isocurvature and postive correla¬ 
tion contributions, or the extra adiabatic power could be can¬ 
celled by a negative cotTelation contribution. We can see that 
these shifts are small. The low-f temperature data continues to 
mildly favour a negative correlation (see in particular the bot¬ 
tom middle panel for each of the three models), since compared 
to the prediction of the best-fit adiabatic base ACDM model, 
the TT angular power at multipoles ^ < 40 is somewhat low. 
But the dotted grey shaded contours in the three middle top pan¬ 
els show that for Planck TT+lowP, the posterior peaks at val¬ 
ues entirely consistent with (0, 0), i.e., the pure adi¬ 

abatic case is preferred. The best-fit values of (P^^, !P®) are 
(1.4 X 10-",4.7x 10-'3)forCDI, (1.2x 10-'^4.6x lO^'®) for 
NDI, and (1.6 x 10-'^2.3 x 10^'°) for NVI, while 
ss (2.4 X 10“®, 2.1 X 10“®). It may appear from the bottom-centre 
panels of Pig. 43 that there is non-zero posterior probability for 
+ 0 when = 0, which would violate the positivity 
constraint, Eq. (124). However, the left-most pixels of the plots 
are actually evaluated at values of large enough so that the 
constraint is satisfied. 

Including the Planck lensing likelihood does not signifi¬ 
cantly affect the non-adiabatic primordial powers, except for 
tightening the constraints on the adiabatic power (see the blue 
versus black contours in the first panels of Pigs. 43 (a), (b), and 
(c)). Including the lensing (C^'^) likelihood constrains the optical 
depth T more tightly than the high-^ temperature and low-^ po¬ 
larization alone (Planck Collaboration XIII, 2015). As there is a 
strong degeneracy between r and the primordial (adiabatic) per¬ 
turbation power Vkk (denoted in the other sections of this paper 
by As), it is natural that adding the lensing data leads to stronger 
constraints on 'Pkr. Moreover, replacing the low-t’ likelihood 


Planck lowP by Planck lowP+WP constrains r better (Planck 
Collaboration XIII, 2015). In the ADI+CDI case the effect of 
this replacement was very similar to adding the Planck lensing 
data (see also Table 15). Although the Planck lensing data do not 
directly constrain the isocurvature contribution,"® they can shift 
and tighten the constraints on some derived isocurvature param¬ 
eters by affecting the favoured values of the standard parameters 
(present even in the pure adiabatic model). However, this effect 
is small, as confirmed by Table 15. Therefore, in the figures we 
do not show ID posteriors of the derived isocurvature parame¬ 
ters for Planck TT+lowP+lensing, since they would be (almost) 
indistinguishable from Planck TT+lowP, as we see in Pig. 43 for 
the primary non-adiabatic parameters. 

In contrast, the high-f polarization data significantly tighten 
the bounds on isocurvature and cross-correlation parameters, as 
seen by comparing the dotted grey and red contours in Pig. 43. 
The significant negative correlation previously allowed by the 
temperature data in the ADI+CDI and ADI+NDI models is now 
disfavoured. This is also clearly visible in the ID posteriors 
of primordial and observable isocurvature and cross-correlation 
fractions, shown, respectively, in Pig. 44 and 45; note how the 
cos A and a-Ri parameters are driven towards zero by the inclu¬ 
sion of the high-f TE,EE data (from the dashed to the solid lines) 
in the ADI+CDI and ADI+NDI cases. We have also checked 
that when the lowP data are replaced by a simple Gaussian prior 
on the reionization optical depth (r = 0.078 + 0.019), the trend 
is similar: the high-f {{ > 30) Planck TT data allow a large 
negative correlation, while the high-^ Planck TE,EE data pre¬ 
fer positive cotTelation. This is clearly seen in Pig. 46 for the 
ADI+CDI case. The best-fit values show an even more dramatic 
effect; we find cos A = -0.55 with TT+lowP, and +0.15 with 
TT,TE,EE+lowP. 


This is expected, since already with Planck TT+lowP, the allowed 
isocurvature fraction is so small that it hardly affects the lensing poten¬ 
tial 
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Fig. 45. Constraints on the fractional contribution of the adiabatic ('R'R), isocurvature (II), and correlation ('RI) components to 
the CMB temperature variance in various multipole ranges, Eq. (132), with Planck TT+lowP data (dashed curves) and with Planck 
TT,TE,EE+lowP data (solid curves). These are shown for the generally-correlated mixed ADI+CDI (black), ADI+NDI (red), or 
ADI+NVI (blue) models. 


Hence, there is a competition between the temperature and 
polarization data that balances out and yields almost symmetric 
results about zero correlation (except in the ADI+NVI case). The 
isocurvature auto-correlation amplitude is also strongly reduced, 
especially in the ADI+CDI case. The best-fit values are slightly 
offset from (!^jj, !^jj) = (0,0), but the pure adiabatic model is 
still well inside the 68 % CL (for ADI+CDI and ADI+NDI) or 
95 % CL (for ADI+NVI) regions. In summary, the high-^ polar¬ 
ization data exhibit a strong preference for adiabaticity, although 
one should keep in mind the possibility of unaccounted-for sys¬ 
tematic effects in the polarization data, possibly leading to artifi¬ 
cially strong constraints. Eor example, the tendency for polariza¬ 
tion to shift the constraints towards positive correlation may be 
due to particular systematic effects that mimic modified acoustic 
peak structure, as we discussed in Sect. 11.1. 

We also performed a parameter extraction with the Planck 
TT,TE,EE+lowP+lensing data, but this combination did not 
provide interesting new constraints—only a tightening of 
bounds on the standard adiabatic parameters, as in the Planck 
TT+lowP+lensing case. 

We provide 95 % CL upper limits or ranges for /3isoj cos A, 
and a-RK in Table 15. With Planck TT+lowP, the constraints 
on the non-adiabatic contribution to the temperature variance. 


1 - aRRi2,250Q), are (-1.5 %, 1.9 %), (-4.0%, 1.4%), and 
(-2.3%, 2.4%) in the ADI+CDI, ADI+NDI, and ADI+NVI 
cases, respectively.^' With Planck TT,TE,EE+lowP these tighten 
to (0.1%, 1.5%), (-0.1%,2.2%), and (-2.0%,0.8%). In the 
ADI+CDI case, zero is not in the 95 % CL interval, but this 
should not be considered a detection of non-adiabaticity. Eor 
example, as mentioned above, (^P^j, Pjj) = (0,0) is in the 
68 % CL region, and the best-fit values are (!P^j.,!P®) = (1.0 x 
10“'^, 3.5 X 10*^). Moreover, the improvement in with re¬ 
spect to the adiabatic model is only 5.3 with 3 extra parame¬ 
ters, so this is not a significant improvement of fit. Indeed, for 
all generally-correlated mixed models the improvement in is 
very small. In particular, with Planck TT+lowP it does not even 
exceed the number of extra degrees of freedom, which is three 
(see Table 15). 

Einally, we checked whether there is any Bayesian evidence 
for the presence of generally-correlated adiabatic and isocurva- 


It should be noted that these numbers can be positive even if the cor¬ 
relation contribution is negative. This happens whenever ajj > la^jl. 
Thus, in the observational non-adiabaticity estimator, 1 - aR«(2,2500), 
the negative numbers are not as pronounced as in the primordial corre¬ 
lation fraction, cos A. 
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Fig. 46. Constraints on the primordial correlation fraction cos A 
in the mixed ADI+CDI model with Planck TT+lowP data 
(dashed black curve) compared to the case where Planck lowP 
data are not used, but replaced by a Gaussian prior t = 0.078 + 
0.019 (dashed red curve). The same exercise is repeated with 
Planck TT,TE,EE data (solid curves), demonstrating that, to a 
great extent, the preferred value of cos A is driven by the high-^ 
data. 


ture modes. In all cases, and with all data combinations stud¬ 
ied, the Bayesian model comparison supports the null hypothe¬ 
sis, i.e., adiabaticity. Indeed, the logarithm of the evidence ra¬ 
tio is InB = ln(Piso/7^ADi) < “5 (i.e., “decisive” evidence, 
or odds of much greater than 150:1, in favour of pure adia¬ 
baticity within Planck’s accuracy and given the parameteriza¬ 
tion and prior ranges used in our analysis), except for ADl-rNDI 
with Planck TT-rlowP-rlensing, for which the evidence ratio is 
slightly larger, -4.6, corresponding to odds of 1:100 for the 
ADIh-NDI model compared to the pure adiabatic model. 

11.3. Robustness of the determination of standard 
cosmoiogicai parameters 

Another spectacular outcome of our analysis is the robustness 
of the determination of the standard cosmological parameters 
against assumptions on initial conditions. Eigure 47 shows the 
ID marginalized posteriors for several cosmological parameters 
(not all independent of each other) with the Planck TT-rlowP 
data alone. Eor the first time, we observe that in the presence of 
one generally-correlated isocurvature mode (CDl, NDI, or NVI), 
predictions for these parameters remain very stable with respect 
to the pure adiabatic case. Except for the ADIh-NDI case, the 
posteriors neither broaden nor shift significantly. A small broad¬ 
ening is only observed in the case of the sound horizon angle, 
0MC. which is naturally the most sensitive parameter to tiny dis¬ 
turbances of the acoustic peak structure. In the ADIh-NDI case, 
the peak of the posterior distribution for some parameters shifts 
slightly, but the largest shift (for Qck^) is less than 1 cr. 

It is particularly striking that the case of a scale-invariant 
adiabatic spectrum {n-R-R - 1) is excluded at many cr even 
when isocurvature modes are present: at 4.7 cr (ADIh-CDI), 5.0 cr 
(ADIh-NDI), 5.4 cr (ADIh-NVI). This illustrates how much the 
constraining power of the CMB has improved. With WMAP 
data, there was still a strong degeneracy between, e.g., the pri¬ 
mordial isocurvature fraction and the adiabatic spectral index 
(Valiviita & Giannantonio, 2009; Savelainen et ak, 2013). This 
degeneracy nearly disappears with Planck TTh-IowP, and even 
more so with Planck TT,TE,EEh-1owP, as shown in the upper 
panel of Pig. 48. Contours in the {n^-j^, cos A) space also shrink 
considerably, with some correlation remaining between these pa¬ 


rameters in the ADIh-CDI and ADIh-NVI cases (Pig. 48, lower 
panel). 

11.4. CDl and primordial tensor perturbations 

A primordial tensor contribution adds extra temperature angu¬ 
lar power at low multipoles, where the adiabatic base ACDM 
model predicts slightly more power than seen in the data. Hence, 
allowing for a non-zero tensor-to-scalar ratio, r, might tighten 
the constraints on positively correlated isocurvature, but degrade 
them in negatively correlated models. We test how treating r as 
a free parameter affects the constraints on isocurvature, and, on 
the other hand, how allowing for the generally-correlated CDl 
mode affects the constraints on r. These cases are denoted as 
“CDiH-r.” Eor comparison, we examine the pure adiabatic case 
in the same parameterization, and call it “ADIn-r.” We also con¬ 
sider another approach where we fix r = 0.1. These cases are 
named as “CDlH-r=0.1” and “ADlH-r=0.1.” 

In the pure adiabatic case (where the curvature and tensor 
perturbations stay constant on super-Hubble scales) the primor¬ 
dial r is the same as the tensor-to-scalar ratio at the Hubble 
radius exit of perturbations during inflation, which we call r. 
However, in the presence of an isocurvature component, 
is not constant in time even on super-Hubble scales (Garcfa- 
Bellido & Wands, 1996). Instead, the isocurvature component 
may source P-r-r, for example, if the background trajectory in 
the field space is curved between Hubble exit and the end of 
inflation (Langlois, 1999; Langlois & Riazuelo, 2000; Gordon 
et al., 2001; Amendola et al., 2002). As a result, we will have at 
the primordial time Prr - PrrHI - cos^ A), where Prr is the 
curvature power at Hubble exit. That is, by the primordial time 
the curvature perturbation power is larger than at the Hubble ra¬ 
dius exit time (Bartolo et al., 2001; Wands et al., 2002; Byrnes &. 
Wands, 2006). Thus we find a relation (Savelainen et al., 2013; 
Valiviita et al., 2012; Kawasaki & Sekiguchi, 2008): 

r = (l-cos^A)f, (134) 

i.e., the tensor-to-scalar ratio at the primordial time (r) is smaller 
than the ratio at the Hubble radius exit time (r). 

The derivation of Eq. (134) assumes that the adiabatic and 
isocurvature perturbations are uncorrelated at Hubble radius exit 
(cos A = 0), and that all the possible primordial correlation 
(cos A it 0) appears from the evolution of super-Hubble per¬ 
turbations between Hubble exit and the primordial time. This is 
true to leading order in the slow-roll parameters, but inflation¬ 
ary models that break slow roll might produce perturbations that 
are strongly correlated already at the Hubble radius exit time. 
In these cases the correlation would depend on the details of 
the particular model, such as the detailed shape of the potential 
and the interactions of the fields. However, a generic prediction 
of slow-roll inflation is that, at Hubble radius exit, the cross- 
correlation Pri is very weak, and indeed is of the order of the 
slow-roll parameters compared to the auto-correlations Prr and 
Pii (see, e.g., Byrnes & Wands, 2006). Thus, for slow-roll mod¬ 
els, I cos A| = (9(slow-roll parameters) « 1. 

In our analysis, we fix the tensor spectral index by the 
leading-order inflationary consistency relation, which now reads 
(Wands et al., 2002) 


f r 

8 8(1 - cos^ A) 


(135) 


Assuming a uniform prior for r would lead to huge negative n^ 
whenever cos^ A was close to one. Therefore, when studying the 
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Fig. 47. Constraints on selected “standard” cosmological parameters with Planck TT+lowP data for the generally-correlated 
ADI+CDl (black), ADl+NDI (red), and ADI+NVl (blue) models compared to the pure adiabatic case (ADI, green dashed curves). 
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Fig. 48. Dependence of the determination of the adiabatic spec¬ 
tral index hk-r (called in the other sections of this paper) on 
the primordial isocurvature fraction ySiso and correlation fraction 
cos A, with Planck TT+lowP data (dashed contours) and with 
Planck TT,TE,EE+lowP data (shaded regions). 


CDI+r case we assume a uniform prior on f at k = 0.05 Mpc ' 
(for details, see Savelainen et al., 2013). 

Surprisingly, allowing for a generally-correlated CDI mode 
(i.e., three extra parameters) hardly changes the constraints on r 
from those obtained in the pure adiabatic model. In Eig. 49 we 
demonstrate this in a “standard” plot of ro.oo2 versus adiabatic 
spectral index. 

Erom Table 15 we notice that, with Planck TT+lowP and 
TT,TE,EE+lowP, fixing r to 0.1 tightens constraints on the pri¬ 
mordial isocurvature fraction at large scales. This is as we ex¬ 
pected, since both tensor and isocurvature perturbations add 
power at low {, and the data do not prefer this. However, the 
shape of the tensor spectrum and correlation spectrum are such 
that negative correlation cannot efficiently cancel the unwanted 
extra power over all scales produced by tensor perturbations (at 
{ < 70). Therefore, the correlation fraction cos A is almost unaf¬ 
fected. However, when we allow r to vary, the cancelation mech¬ 
anism works to some degree when using Planck TT+lowP data, 
leading to more negative cos A than without r: with varying r 
we have cos A in the range (-0.43,0.20), while without r it is 
in (-0.30,0.20), at 95 % CL. As there is now some cancella¬ 
tion of power at large scales, the constraint on /3iso(kiow) weakens 
slightly from 0.041 without r to 0.043 with r. On the other hand, 
the high-^ polarization data constrain the correlation to be so 
close to zero that with Planck TT,TE,EE+lowP the results for 
cos A with and without r are almost identical. 

The mean value of cos A in the CDI+r cases is -0.071 
(TT+lowP) and -0.076 (TT,TE,EE+lowP). Therefore, 1 - 
cos^ A 0.99, and so we do not expect a large difference be¬ 
tween the primordial r and the Hubble radius exit value, f. The 
smallness of the difference is evident in Table 14. To summarize, 
CDI hardly affects the determination of r from the Planck data, 
and allowing for tensor perturbations hardly affects the determi¬ 
nation of the non-adiabaticity parameters. 

11.5. Special CDI cases 

Next we study three one-parameter CDI extensions to the adia¬ 
batic model. In all these extensions the isocurvature mode mod¬ 
ifies only the largest angular scales, since we either fix njj 
to unity (“axion”) or to the adiabatic spectral index (“curvaton 
I/II”). As can be seen from Pig. 42, the polarization E mode at 
multipoles £ > 200 will not be significantly affected by this type 
of CDI mode. Therefore, these models are much less sensitive 
to residual systematic effects in the high-f” polarization data than 
the generally correlated models. 
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n-mz 

Fig. 49. 68 % and 95 % CL constraints on the primordial adia¬ 
batic spectral index and the primordial tensor-to-scalar ratio 
r (more accurately, in the CDI+r model, the primordial tensor- 
to-curvature power ratio) ?A k - 0.002 Mpc '. Filled contours 
are for generally-correlated ADI+CDI and solid contours for the 
pure adiabatic model. 


Model (and data) 

ro.05 

ro.05 

^tens /^scal 
'^10 /'^lO 

CDI+r (TT+lowP) 

0.086 

0.089 

0.041 

ADI+r (TT+lowP) 

0.101 

0.101 

0.048 

CDI+r (TT,TE,EE+lowP) 

0.092 

0.092 

0.043 

ADI+r (TT,TE,EE+lowP) 

0.094 

0.094 

0.044 


Table 14. 95 % CL upper bounds on the tensor-to-scalar ratio 
(actually the tensor-to-curvature power ratio) at the primordial 
time, r, and earlier, at the Hubble radius exit time during in¬ 
flation, r, ?Ak - 0.05 Mpc '. (In the pure adiabatic case r and 
r are equal.) In the last column indicates the tensor 

contribution to the temperature angular power at ^ = 10 rela¬ 
tive to the temperature power from scalar perturbations (Cj™* = 

qRI qIR qII-. 

*^10 ^ ^ *^10 ^ *^10 


11.5.1. Uncorrelated ADI+CDI (“Axion”) 

We start with an uncorrelated mixture of adiabatic and CDI 
modes CPri - 0), and make the additional assumption that 
^Ji “ assume unit isocurvature spectral index, 

nil - L Constraints in the (n^K^Aso) plane are presented in 
Fig. 50. This model is the only case for which our new results do 
not improve over bounds from PCI13. At kmid = 0.050Mpc *, 
we find /?iso < 0.038 (95 % CL, TT,TE,EE+lowP; see Table 15), 
compared with ySjso < 0.039 using Planck 2013 and low-^ 
WMAP data. This is not surprising, since fixing «jj to unity im¬ 
plies that bounds are dominated by measurements on very large 
angular scales, { < 30, as can easily be understood from Pig. 42. 
Hence the results are insensitive to the addition of better high-f 
temperature data, or new high-f polarization data. 

We summarized in PCI 13 why an uncorrelated CDI mode 
with Mjj ss 1 can be produced in axion models, under a number 
of restrictive assumptions: the Peccei-Quinn symmetry should 
be broken before inflation; it should not be restored by quantum 



6 0.95 0.96 0.97 0.98 
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Fig. 50. Uncorrelated ADI+CDI with njj = 1 (“axion”). 

fluctuations of the inflaton or by thermal fluctuations when the 
Universe reheats; and axions produced through the misaligne- 
ment angle should contribute to a sizable fraction (or all) of the 
dark matter. Under all of these assumptions, limits on ySiso can 
be used to infer a bound on the energy scale of inflation, using 
equation (73) of PCI13. This bound is strongest when all the 
dark matter is assumed to be in the form of axions. In that case, 
the limit on j0iso(kmid) for Planck TT,TE,EE+lowP gives 

/ r \ 0.408 

//inf < 0.86 X 10" Q^^ LoiiGev ) 

where //inf is the expansion rate at Hubble radius exit of the scale 
corresponding to k^a - 0.050 Mpc * and fa is the Peccei-Quinn 
symmetry breaking energy scale. 

11.5.2. Fully correlated ADI+CDI (“Curvaton I”) 

Another interesting special case of mixed adiabatic and CDI (or 
BDI) perturbations is a model where these perturbations are pri- 
mordially fully correlated and their power spectra have the same 
shape. These cases are obtained by setting 
which, by condition (129), implies that the corresponding state¬ 
ment holds at scale k^ and indeed at any scale. In addition, we 
set From this it follows 

that ySiso is scale independent. Therefore, this model has only one 
primary non-adiabaticity parameter, 

A physically motivated example of this type of model is the 
curvaton model (Mollerach, 1990; Linde & Mukhanov, 1997; 
Enqvist & Sloth, 2002; Moroi & Takahashi, 2001; Lyth & 
Wands, 2002; Lyth et al., 2003) with the following assumptions. 
(1) The average curvaton field value, is sufficiently below 
the Planck mass when cosmologically interesting scales exit the 
Hubble radius during inflation. (2) At Hubble radius exit, the 
curvature perturbation from the inflaton is negligible compared 
to the perturbation caused by the curvaton. (3) The same is true 
for any inflaton decay products after reheating. This means that, 
after reheating, the Universe is homogeneous, except for the spa¬ 
tially varying entropy (i.e., isocurvature perturbation) due to the 
curvaton field perturbations. (4) Later, CDM is created from the 
curvaton decay and baryon number after curvaton decay. This 
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Fig. 51. Fully correlated ADI+CDI with tijj - (“curvaton 
!”)■ Since the spectral indices are equal, the primordial isocurva¬ 
ture fraction is scale independent. 


corresponds to case 4 presented in Gordon &. Lewis (2003). (5) 
The curvaton contributes a significant amount to the energy den¬ 
sity of the Universe at the time of the curvaton’s decay to CDM, 
i.e., the curvaton decays late enough. (6) The energy density of 
curvaton particles possibly produced during reheating should be 
sufficiently low (Bartolo & Liddle, 2002; Linde & Mukhanov, 
2006). (7) The small-scale variance of curvaton perturbations, 
- {^X^)slx ^7 is negligible, so that it does not significantly 
contribute to the average energy density on CMB scales; see 
equation (102) in Sasaki et al. (2006). The last two conditions are 
necessary in order to have an almost-Gaussian curvature pertur¬ 
bation, as required by the Planck observations. Namely, if they 
are not valid, a large follows, as discussed below. Indeed, 
the conditions (6) and (7) are related, since curvaton particles 
would add a homogeneous component to the average energy 
density on large scales, and hence we can describe their effect 
by = P;(-, particles/P;P, field, where P;f, field is the average energy 
density of the classical curvaton held on large scales; see equa¬ 
tion (98) in Sasaki et al. (2006). Then the total energy density 
earned by the curvaton will be - p^, field + P;t', particles- 

The amount of isocurvature and non-Gaussianity present af¬ 
ter curvaton decay depends on the “curvaton decay fraction”. 


ro 


^Px 


4" 4pradiation 


(137) 


evaluated at curvaton decay time. If conditions (6) and (7) do not 
hold, then the isocurvature perturbation disappears. 

The curvaton scenario presented here is one of the sim¬ 
plest to test against observations. It should be noted that at least 
the conditions (l)-(5) listed at the beginning of this subsection 
should be satished simultaneously. Indeed, if we relax some of 
these conditions, almost any type of correlation can be produced. 

Indeed, if curvaton particles are produced during reheating, they 
can be expected to survive and outweigh other particles at the moment 
of curvaton decay, but by how much depends on the details of the model. 
As the curvaton field (during its oscillations) and the curvaton parti¬ 
cles have the same equation of state and they decay simultaneously, no 
isocurvature perturbations are produced. 


Fig. 52. Fully anticorrelated ADIh-CDI with njj - (“curva¬ 
ton 11”). 


For example, the relative correlation fraction can be written as 
cos A = -\//l/(l -H A), where A - (8/9)r^et(Mpi/^,)^. Therefore, 
the model is fully correlated only if T » 1. If the slow-roll pa¬ 
rameter e» is very close to zero or the curvaton held value x* is 
large compared to the Planck mass, this model leads to almost 
uncorrelated perturbations. 

As seen in Fig. 51 and Table 15, the upper bound on the pri¬ 
mordial isocurvature fraction in the fully correlated ADIh-CDI 
model weakens slightly when we add the Planck lensing data 
to Planck TTh-IowP, whereas adding high-^ TE,EE tightens the 
upper bound moderately. With all of these three data combina¬ 
tions, the pure adiabatic model gives an equally good best-ht 
as the fully-correlated ADIh-CDI model. Bayesian model com¬ 
parison strengthens the conclusion that the data disfavour this 
model with respect to the pure adiabatic model. 

The isocurvature fraction is connected to the curvaton decay 
fraction, Eq. (137), by 


9(1 - rp) 

D + 9(1 - rof 


(138) 


(see case 4 in Gordon & Lewis, 2003). Reading the constraints 
on ySiso from Table 15, we can convert them into constraints on 
ro and further into the non-Gaussianity parameter assuming a 
quadratic potential for the curvaton and instantaneous decay“^ 
(Sasaki et al., 2006): 


^local 

/nl 


= (1 + A?) 
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(139) 


If conditions (6) and (7) hold, i.e., A^ = 0, as implicitly assumed, 
e.g., in Bartolo et al. (2004a,b), then the smallest possible value 


It should be noted that, in particular, in the older curvaton literature 
a formula = 4 ^ is often quoted or utilized. This result, which fol¬ 
lows from considering only squares of first order perturbations, is valid 
when td is close to zero, i.e., when is very large. However, when 
To is close to unity or < 10, which is the case with the Planck 
measurements, the second and third terms in Eq. (139) are vitally im¬ 
portant. These follow from genuine second order perturbation theory 
calculations. Coincidentally, if one erroneously uses the expression ^ 
in the limit —> 1, one obtains a result h- 5/4, whereas the correct 
formula (139), with Aj = 0, leads to -5/4, when ro —> 1. 
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of is -5/4, which is obtained when ro - 1, and Eqs. (138) 
and (139) yield for the various Planck data sets (at 95 % CL):^'^ 

TT+lowP: ySiso < 0.0018 ^ 0.9860 < rz, < 1 

^ -1.250 < /i°“' < -1.220, (140) 

TT+lowP+lensing: < 0.0022 ^ 0.9845 < < 1 

^-1.250 <-1.217, (141) 

TT,TE,EE+lowP: < 0.0013 ^ 0.9882 < rz, < 1 

^ -1.250 < <-1.225. (142) 

Thus the results for the simplest curvaton model remain un¬ 
changed from those presented in PCI13, i.e., in order to pro¬ 
duce almost purely adiabatic perturbations, the curvaton should 
decay when it dominates the energy density of the Universe 
(rz, > 0.98), and the non-Gaussianity parameter is constrained 
to close to its smallest possible value (-5/4 < < -1.21), 

which is consistent with the result = 2.5 + 5.7 (68 % CL, 
from T only) found in Planck Collaboration XVII (2015). 


11.5.3. Fully anticorrelated ADI+CDI (“Curvaton M”) 


The curvaton scenario or some other mechanism could also 
produce 100% anticorrelated perturbations, with njj = 

The constraints in the («■«•;« >Aso) plane are presented in Pig. 52. 
Examples of this kind of model are provided by cases 2, 3, and 6 
in Gordon & Lewis (2003). These lead to a hxed, large amount 
of isocurvature, e.g., in case 2 to /liso = 9/10, and are hence 
excluded by the data at very high signihcance. However, case 
9 in Gordon & Lewis (2003), with a suitable ro (i.e., rz> > Rc, 
where - pdipz + Ph)), leads to fully anticorrelated perturba¬ 
tions and might provide a good ht to the data. In this case CDM 
is produced by curvaton decay while baryons are created ear¬ 
lier from inflaton decay products and do not carry a curvature 
perturbation. We obtain a very similar expression to Eq. (138), 
namely 




9(1 - rp/R^f 
rl + 9(l-rD/Rc)^ ' 


(143) 


We convert this approximately to a constraint on ro by hxing Rc 
to its best-fit value, Rc - 0.8437 (Planck TT+lowP), within this 
model. The results for the various Planck data sets are: 


TT+lowP: ySiso < 0.0064 0.8437 < rz, < 0.8632 

^ -0.9379 < < -0.8882, (144) 

TT+lowP+lensing: < 0.0052 0.8437 < rz, < 0.8612 

-0.9329 < < -0.8882, (145) 

TT,TE,EE+lowP: Aso < 0.0008 0.8437 < rz, < 0.8505 

=> -0.9056 < < -0.8882. (146) 

After all the tests conducted in this section, both for the gen¬ 
erally correlated GDI, NDI, and NVI cases as well as for the 
special GDI cases, we conclude that, within the spatially flat base 
AGDM model, the initial conditions of perturbations are consis¬ 
tent with the hypothesis of pure adiabaticity, a conclusion that is 
also overwhelmingly supported by the Bayesian model compar¬ 
ison. Moreover, Planck Gollaboration XVII (2015) report a null 

However, if was non-negligible, then all the constraints on 
would shift upward. For example, with A^ = 1, our constraints 
on ySiso would translate to 0 < < 0.03. On the other hand, 

the Planck constraint of ///“’ can be converted to an upper bound 
A? = P^,panicie.s/P^,fieid < 8.5 (95 % CL from T only) as shown in Planck 
Collaboration XVII (2015). 


detection of isocurvature non-Gaussianity, with polarization im¬ 
proving constraints significantly. 


12. Statistical anisotropy and inflation 

A key prediction of standard inflation, which in the present 
context includes all single held models of inflation as well as 
many multi-field models, is that the stochastic process gener¬ 
ating the primordial cosmological perturbations is completely 
characterized by its power spectrum, constrained by statistical 
isotropy to depend only on the multipole number {. This state¬ 
ment applies at least to the accuracy that can be probed using the 
GMB, given the limitations imposed by cosmic variance, since 
all models exhibit some level of non-Gaussianity. Nevertheless, 
more general Gaussian stochastic processes can be envisaged for 
which one or more special directions on the sky are singled out, 
so that the expectation values for the temperature multipoles take 
the form 

= (147) 

rather than the very special form 

(aL(4m')*) = Crdz.rd,„,„z, (148) 

which is the only possibility consistent with statistical isotropy. 

The most general form for a Gaussian stochastic process 
on the sphere violating the hypothesis of statistical isotropy, 
Eq. (147), is too broad to be useful, especially given the fact that 
we have only one sky to analyse. Eor ( < fmax, there are 0{{^ca) 
multipole expansion coefficients, compared with model 

parameters. Therefore, in order to make some progress on test¬ 

ing the hypothesis of statistical isotropy, we must restrict our¬ 
selves to examining only the simplest models violating statis¬ 
tical isotropy, for which the available data can establish mean¬ 
ingful constraints and for which one can hope to And a simple 
theoretical motivation. 

12.1. Asymmetry: observations versus model building 

In one simple class of statistically anisotropic models, we start 
with a map produced by a process respecting statistical isotropy, 
which becomes modulated by another held in the following man¬ 
ner to produce the observed sky map: 

ST^^yiCl) = (l + M(Q)) dT^dCl), (149) 

where Q denotes a position on the celestial sphere and A7’s.i(Q) is 
the outcome of the underlying statistically isotropic process be¬ 
fore modulation. Roughly speaking, where the modulating held 
M(Cl) is positive, power on scales smaller than the scale of vari¬ 
ation of M(Cl) is enhanced, whereas where M(Cl) is negative, 
power is suppressed. We refer to this as a “power asymmetry.” 
If M(Q) = Ad ■ n, we have a model of dipolar modulation with 
amplitude A and direction d, but higher-order or mixed modula¬ 
tion may also be considered, such as a quadrupole modulation or 
modulation by a scale-invariant held M(£l), to name just a few 
special cases. 

In Planck Gollaboration XXIII (2014) and Planck 
Gollaboration XVI (2015), the details of constructing effi¬ 
cient estimators for statistical anisotropy, in particular in the 
presence of realistic data involving sky cuts and possibly in¬ 
completely removed foreground contamination, are considered 
in depth. In addition, the question of the statistical significance 
of any detected “anomalies” from the expectations of base 
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Model (and data) 

100/ii,„(k,„„) 

100/3i,Afem.d) 

100/3iso(khigh) 

100 cos A 

100a:«K(2,2500) 

A« 

At" 

InS 

General models: 

GDI (TT+lowP) 

4.1 

35.4 

56.9 

[-30:20] 

[98.1:101.5] 

3 

-2.1 

-8.8 

GDI (TT+lowP+WP) 

4.2 

35.5 

57.2 

[-31:23] 

[97.9:101.4] 

3 

-1.8 

-9.1 

GDI (TT,TE,EE+lowP) 

2.0 

[3.4:28.1] 

[3.1:51.8] 

[ -6:20] 

[98.5: 99.9] 

3 

-5.3 

-8.8 

GDI (TT,TE,EE+lowP+WP) 

2.1 

[2.3:28.4] 

[2.6:52.1] 

[ -7:21] 

[98.5: 99.9] 

3 

-5.5 

-8.2 

GDI (TT+lowP+lensing) 

4.5 

37.9 

59.4 

[-28:17] 

[98.1:101.1] 

3 

-1.2 

-8.8 

NDI (TT+lowP) 

14.3 

22.4 

27.4 

[-33: 1] 

[98.6:104.0] 

3 

-2.0 

-5.3 

NDI (TT,TE,EE+lowP) 

7.3 

[3.4:19.3] 

[3.5:26.7] 

[ -9:10] 

[97.8:100.1] 

3 

-5.5 

-5.5 

NDI (TT+lowP+lensing) 

15.8 

[1.4:24.1] 

[0.3:28.4] 

[-32: 0] 

[98.6:104.0] 

3 

-2.8 

-4.6 

NVI (TT+lowP) 

8.3 

[0.1:10.2] 

11.9 

[-26: 6] 

[97.6:102.3] 

3 

-2.8 

-6.3 

NVI (TT,TE,EE+lowP) 

7.4 

[0.9: 7.4] 

[0.4: 8.8] 

[-22:-4] 

[99.2:102.0] 

3 

-6.2 

-6.5 

NVI (TT+lowP+lensing) 

9.7 

[0.4:11.6] 

13.1 

[-23: 7] 

[97.1:102.0] 

3 

-2.5 

-6.5 

General models + r: 

GDI+r=0.1 (TT+lowP) 

3.4 

38.7 

63.9 

[-33:24] 

[98.1:101.4] 

3 

-5.4 

-8.9 

GDI+r=0.1 (TT,TE,EE+lowP) 

1.6 

[4.4:31.7] 

[6.9:59.2] 

[ -6:22] 

[98.6: 99.9] 

3 

-6.3 

-8.1 

GDI+r (TT+lowP) 

4.3 

34.9 

56.2 

[-43:20] 

[97.9:102.4] 

3 

-3.3 

-7.7 

GDI+r (TT,TE,EE+lowP) 

1.7 

[3.9:29.0] 

[5.8:53.8] 

[ -5:21] 

[98.6: 99.9] 

3 

-5.1 

-7.2 

Special GDI cases: 

Uncorrelated, njj = 1 
“axion” (TT+lowP) 

3.3 

3.7 

3.8 

0 

[98.5:100] 

1 

0.0 

-5.2 

“axion” (TT,TE,EE+lowP) 

3.5 

3.8 

3.9 

0 

[98.4:100] 

1 

-0.2 

-4.9 

“axion” (TT+lowP+lensing) 

3.9 

4.3 

4.4 

0 

[98.3:100] 

1 

0.0 

-5.0 

Eully correlated, njj = 

“curvaton 1” (TT+lowP) 

0.18 

0.18 

0.18 

100 

[97.5:100.0] 

1 

-0.1 

-8.1 

“curvaton I” (TT,TE,EE+lowP) 

0.13 

0.13 

0.13 

100 

[97.8: 99.9] 

1 

0.0 

-7.8 

“curvaton I” (TT+lowP+lensing) 

0.22 

0.22 

0.22 

100 

[97.3: 99.7] 

1 

0.0 

-8.5 

Eully anti-correlated, njj = 

“curvaton 11” (TT+lowP) 

0.64 

0.64 

0.64 

-100 

[100.5:105.1] 

1 

-1.1 

-5.4 

“curvaton IT’ (TT,TE,EE+lowP) 

0.08 

0.08 

0.08 

-100 

[100.1:101.8] 

1 

0.0 

-8.9 

“curvaton IT’ (TT+lowP+lensing) 

0.52 

0.52 

0.52 

-100 

[100.4:104.4] 

1 

-0.6 

-6.3 


Table 15. Constraints on mixed adiabatic and isocurvature models. For each mixed model, we report 95% CL bounds on the frac¬ 
tional primordial contribution of isocurvature modes at three comoving wavenumbers (k\ov, = 0.002 Mpc ' , - 0.050 Mpc ^ , 

and A:high = 0.100 Mpc '), as well as the scale-independent primordial correlation fraction, cos A. The fractional adiabatic contribu¬ 
tion to the observed temperature variance is denoted by aKK(2,2500), and from this the nonadiabatic contribution can be calculated 
as anon-adi = 1 - Q'rk( 2, 2500). The number of extra parameters compared with the corresponding pure adiabatic model is denoted 
by An, and is the difference between the x^ of the best-htting mixed and pure adiabatic models. (A negative Ax^ means that 
the mixed model is a better ht to the data.) In the last column we give the difference between the logarithm of Bayesian evi¬ 
dences. (A negative InZ? = ln(Piso/FADi) means that Bayesian model comparison disfavours the mixed model. With our settings of 
MultiNest the uncertainty in these numbers is about ±0.5.) 


ACDM is examined in detail. Importantly, in the absence of 
a particular inflationary model for such an observed anomaly, 
the signihcance should be corrected for the “multiplicity of 
tests” that could have resulted in similarly signihcant detections 
(i.e., for the “look elsewhere effect”), although applying such 
corrections can be ambiguous. In this paper, however, we 
consider only forms of statistical anisotropy that are predicted 
by specihc inflationary models, and, hence, such corrections 
will not be necessary. 

Several important questions can be posed regarding the link 
between statistical isotropy and inflation. In particular, we can 
ask the following questions. (1) Does a statistically significant 
finding of a violation of statistical isotropy falsify inflation? (2) 
If not, what sort of non-standard inflation could produce the re¬ 
quired departure from statistical isotropy? (3) What other per¬ 
haps non-inflationary models could also account for the violation 
of statistical isotropy? In this section, we begin to address these 
questions by assessing the viability of an inflationary model for 
dipolar asymmetry, as well as by placing new limits on the pres¬ 
ence of quadrupolar power asymmetry. 


For the case of the dipolar asymmetry reported in Planck 
Collaboration XVI (2015), there are two aspects that make in¬ 
flationary model building difficult. First is the problem of ob¬ 
taining a significant amplitude of dipole modulation. In Planck 
Collaboration XVI (2015) the asymmetry was found to have am¬ 
plitude A 6-7 % on scales 2 <{ < 64. This compares with 
the expected value of A = 2.9 % on these scales due to cosmic 
variance in statistically isotropic skies. One basic strategy for in¬ 
corporating the violation of statistical isotropy into inflation is 
to consider some form of multi-field inflation and use one of the 
directions orthogonal to the direction of slow roll as the held 
responsible for the modulation. Obtaining the required modula¬ 
tion is problematic because most extra fields in multi-field in¬ 
flation become disordered in a nearly scale-invariant way, just 
like the fluctuations in the held parallel to the direction of slow 
roll. What is needed resembles a pure gradient with no fluctu¬ 
ations of shorter wavelength. In Liddle & Cortes (2013) it was 
proposed that such a held could be produced using the super¬ 
curvature mode of open inflation. (See, however, the discussion 
in Kanno et al. (2013).) Also, in order to respect the /nl con- 
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straints, one must avoid that the modulating field leave a direct 
imprint on the temperature anisotropy. 

The second aspect which makes model building difficult for 
dipolar asymmetry is that the measured amplitude is strongly 
scale dependent, and on scales i > 100 no significant detection 
of a dipolar modulation amplitude is made (Planck Collaboration 
XVI, 2015), once our proper motion has been taken into account 
(Planck Collaboration XXVII, 2014). On the other hand, the 
simplest models are scale-free and produce statistical anisotropy 
of the type described by the ansatz in Eq. (149), for which the 
bulk of the statistical weight should be detected at the reso¬ 
lution of the survey. To resolve this difficulty, Erickcek et al. 
(2009) proposed modulating CDI fluctuations generated within 
the framework of a curvaton scenario, because, unlike adiabatic 
perturbations, CDI perturbations entering the horizon after last 
scattering contribute negligibly to the CMB fluctuations (recall 
Fig. 42). 

The situation for the quadrupolar power asymmetry is dif¬ 
ferent from the dipolar case in that no detection is currently 
claimed. Model building is easier than the dipolar case since 
no pure gradient modes are required, but also more difficult in 
that anisotropy during inflation is needed. While the isotropy 
of the recent expansion of the Universe (i.e., since the CMB 
fluctuations were first imprinted) is tightly constrained, bounds 
on a possible anisotropic expansion at early times are much 
weaker. Ackerman et al. (2007) proposed using constraints on 
the quadrupolar statistical anisotropy of the CMB to probe the 
isotropy of the expansion during inflation—that is, during the 
epoch when the perturbations now seen in the CMB first exited 
the horizon. Assuming an anisotropic expansion during infla¬ 
tion, Ackerman et al. (2007) computed its impact on the three- 
dimensional power spectrum on superhorizon scales by integrat¬ 
ing the mode functions for the perturbations during inflation and 
beyond. Several sources of such anisotropy have been proposed, 
such as vector fields during inflation (Dimastrogiovanni et al., 
2010; Soda, 2012; Maleknejad et al., 2013; Schmidt & Hui, 
2013; Bartolo et al., 2013; Naruko et al., 2014), or an inflating 
solid or elastic medium (Bartolo et al., 2013). 


12.2. Scale-dependent modulation and idealized 
estimators 


The ansatz in Eq. (149) expressed in angular space may be 
rewritten in terms of the multipole expansion and generalized 
to include scale-dependent modulation by means of Wigner 3 j 
symbols: 


+L 


(4.4,,) = Z Z '^4 


L=0 M=-L 


{ L 
1 m m! M 


(150) 


Because of the symmetry of the left-hand side, the coefficients 
M acquire a phase under interchange of { and {'. 

This is the most general form consistent with the hypothesis of 
Gaussianity. The usual isotropic power spectrum, which is the 
generic prediction of simple models of inflation, includes only 
the L - 0 term, where g = and the Wigner 3 j symbol 

provides the factor. The coefficients with L > 0 

introduce statistical anisotropy. 

If we assume that there is a common vector (corresponding 
to L = 1 on the celestial sphere) that defines the direction of the 
anisotropy of the power spectrum for all the terms of L = 1, we 
may adopt a more restricted ansatz for the bipolar modulation, 
so that 


where we assume that Xm is normalized (i.e., - !)■ 

In such a theory, supposing that C] ^ is theoretically determined, 
but the orientation of the unit vector Xm is random and isotrop¬ 
ically distributed on the celestial sphere, we may construct the 
following quadratic estimator for the direction: 


y(L) ^ _ 

" (2i+l)(Q)F2(Q0F2 

le {' L\ J T 

^\m m' Mj ’ 


(152) 


where the weights for the unbiased minimum variance estimator 
are given by 


we,e'-,L 


ce/'-.L 


Yj^e.e'-L 


e.c y 


(153) 


This construction, which for the L - I case may be found in 
Moss et al. (2011) and Planck Collaboration XVI (2015), may 
be readily generalized to L > 1 in the above way. 


12.3. Constraining inflationary models for dipolar 
asymmetry 

In this section, we confront with Planck data the modulated cur¬ 
vaton model of Erickcek et al. (2009), which attempts to explain 
the observed large-scale power asymmetry via a gradient in the 
background curvaton field. In this model, the curvaton decays 
after CDM freeze-out, which results in a nearly scale-invariant 
isocurvature component between CDM and radiation. In the vi¬ 
able version of this scenario, the curvaton contributes negligibly 
to the CDM density. A long-wavelength fluctuation in the curva¬ 
ton field initial value cr, is assumed, with amplitude Act* across 
our observable volume. This modulates the curvaton isocurva¬ 
ture fluctuations according to Sa-y ^ 26cr,lcr,. The curvaton pro¬ 
duces all of the final CDI fluctuations, which are nearly scale in¬ 
variant, as well as a component of the final adiabatic fluctuations. 
Hence both of these components will be modulated, and the pa¬ 
rameter space of the model will be constrained by observations 
of the power asymmetry on large and small scales, as well as the 
full-sky CDI fraction. In practice, the very tight constraints on 
small-scale power asymmetry obtained in Planck Collaboration 
XVI (2015) imply a small curvaton adiabatic component, which 
implies that the CDI and adiabatic fluctuations are only weakly 
correlated. This model easily satisfies constraints due to the 
CMB dipole, quadrupole, and non-Gaussianity (Erickcek et al., 
2009). 

There are two main parameters which we constrain for this 
model. Eirst, the fraction of adiabatic fluctuations due to the cur¬ 
vaton, is defined as 





(154) 


Here, ’P-Rm and Pa- are the inflaton and curvaton primordial 
power spectra, respectively, and Eo- is the coupling from cur¬ 
vaton isocurvature to adiabatic fluctuations. (Up to a sign, ^ 
is equal to the correlation parameter.) Next, the coupling of 
curvaton to CDI, Mqdio-, is determined by the constant k = 
McDia-IR ^ “1^ where 


R = 


3Q^ 


tt _ ^1 ytP 


(151) 


4Qy + "I" 3 Qcdm 
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and all density parameters are evaluated just prior to curvaton 
decay. The isocurvature fraction can be written in terms of these 
two parameters by 

- rr^- 

These parameters determine the modulation of the CMB temper¬ 
ature fluctuations via AQ/Q = IKeAa-flo-t, where (Erickcek 
et ah, 2009) 


Kt = 




Cf + 9 k^O^° + 3/rq” 
Cf +^{9K^Op + 3 kC]°^)' 


(157) 


Here C’f, C™, and are the adiabatic, CDI, and coiTelated 
power spectra calculated for unity primordial spectra. 

Note that this modulated curvaton model contains some sim¬ 
ple special cases. For /c = 0, we have a purely adiabatic, i.e., 
scale-invariant, modulation. This is equivalent to a modulation 
of the scalar amplitude As. On the other hand, if we take the 
limit K ^ oo, with fixed (i.e., with fixed isocurvature frac¬ 
tion ySiso), we obtain a pure CDI modulation. For k - ^ - 0 we 
have no modulation, i.e., we recover base ACDM. Therefore this 
model is particularly useful for examining a range of possible 
modulations within the context of a concrete framework. 

In order to constrain this model, we use a formalism which 
was developed to determine the signatures of potential gradi¬ 
ents in physical parameters in the CMB (Moss et ah, 2011), and 
which is used to examine dipolar modulation and described in 
detail in Planck Collaboration XVI (2015). This approach is well 
suited to testing the modulated curvaton model since it can ac¬ 
commodate scale-dependent modulations. Briefly, we write the 
temperature anisotropy covariance given a gradient AXm in a pa¬ 
rameter X as 


C, 


£mf't 


= Ctda-S,„,„ + (-1)'”^ [{2(+l){2(' + 1)]'^2 


{ {' 1 
0 0 0 




{ {’ 1 
-m m' Ml’ 


(158) 


where 6Cw = dCddX+dCp/dX. Note that this covariance takes 
the form of Eqs. (150) and (151), with 

Cit' = ^ [(2^ + 1)(2^' + 1)]'^' (o 0 o) • 


We then construct a maximum likelihood estimator for 
the gradient components. We use C ^ Altered data (Planck 
Collaboration XV, 2015) and perform a mean-field subtraction, 
giving 


AXm = 


flA 




■CmC'm' 


( (' 1 
■m m' M 


tt’ 


(160) 


Here f\M is a normalization correction due to the applied mask, 
M(Q.), and is given by 


flM 




£/Q Tr„(Q)M(Q). 


(161) 


The T(m are the Altered data and Ft = In practice, ex¬ 

ploring the parameter space of the model is sped up dramatically 
by binning the estimator, Eq. (160), into bins of width A£ - 1, 



Fig. 53. 68 % and 95 % CF regions in the modulated isocurva¬ 
ture model parameter space using the Planck temperature data 
up to ^max = 1000 (contours). The region above the dashed curve 
is ruled out by the Planck constraint on an uncorrelated, scale- 
invariant isocurvature component. 


which means that the estimators only need to be calculated once 
(Planck Collaboration XVI, 2015). Finally, for the modulated 
curvaton model we identify 


dCe 


= 2KeCt. 


(162) 


Note that, for our constraints, we Ax the curvaton gradient to 
its maximum value, Aa-^jcr^ - 1. Therefore, our constraints are 
conservative, since smaller Aa-,lcr, would only reduce the mod¬ 
ulation that this model could produce. 

The temperature anisotropies measured by Planck constrain 
the modulated curvaton parameters k and ^ via Eqs. (157) and 
(160). Figure 53 shows the constraints in this parameter space 
evaluating the estimator to fmax = 1000. The maximum like¬ 
lihood region corresponds to a band at /c > 3. For parame¬ 
ters in this region, the model produces a large-scale asymme¬ 
try via a mainly-CDI modulation. However, the amplitude of 
this large-scale asymmetry is lower than the 6-7 % actually ob¬ 
served (Planck Collaboration XVI, 2015). The reason is that, had 
a CDI modulation produced all of the large-scale asymmetry, 
the consequent small-scale asymmetry (due to the shape of the 
scale-invariant CDI spectrum) would be larger than the Planck 
observations allow. The allowed CDI modulation is further re¬ 
duced by the Planck 95 % upper limit on an uncorrelated, scale- 
invariant (“axion”-type) isocurvature component, ySiso < 0.0331, 
from Sect. 11. Imposing this constraint reduces the available pa¬ 
rameter space in the plane via Eq. (156), as illustrated in 
Fig. 53. 

The best At in Fig. 53 corresponds to Axlg - -6.8 rela¬ 
tive to base ACDM, for two extra paramters. In order to as¬ 
sess how likely such an improvement would be in statistically 
isotropic skies, we note that the best-fit CDI modulation ampli¬ 
tude is very close to the mean amplitude expected due to cos¬ 
mic variance, as calculated directly from Eq. (160). More pre¬ 
cisely, since the amplitude is distributed with three degrees 
of freedom, i.e., Maxwell-Boltzmann distributed, we conclude 
that about 44 % of statistically isotropic skies will exhibit a mea¬ 
sured [via Eq. (160)] isocurvature modulation larger than that of 
the actual sky. 

To summarize, the modulated curvaton model can only pro¬ 
duce a small part of the observed large-scale asymmetry, and 
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what it can produce is entirely consistent with cosmic variance in 
a statistically isotropic sky. Hence we must favour base ACDM 
over this model. Finally, note that further generalizing the model 
(e.g., to allow non-scale-invariant CDI spectra) may allow more 
large-scale asymmetry to be produced and hence result in an im¬ 
proved at the expense of more parameters. On the other 

hand, the neutrino isocurvature modes are not expected to fit the 
observed asymmetry well, due to their approximate scale invari¬ 
ance (see Fig. 42). 


12.4. Constraints on quadrupolar asymmetry generated 
during inflation (L-2) 

In this section, we assume a quadrupolar direction dependence 
in the primordial power spectrum, which we may expand using 
spherical harmonics according to 

g,(k) (ic ■ ^ Y2MCk). (163) 

M 

With this expansion, we then write the power spectrum as 


Pik) = Paik) 


1 + ^ glMfk) Y2M(k) 
M 


(164) 


where we have absorbed gt{k)ll> into the normalization of the 
isotropic part, Po{k) = f’o(^)(l + g*(k)l3), and defined 


giuik) = 


87t g,{k) 

15 1 + g.(k)/3 


nuiEu) 


8 ;r 

15 


g^Yl^iKi), 


(165) 


with g 2 M{k) satisfying g 2 ,-M{k) = (-1)^^ g^j^jik). The approx¬ 
imation on the rightmost side is justified by the fact that the 
amplitude g* is observationally constrained to quite a small 
value, gt < 10“^ (Kim & Komatsu, 2013). In this analysis, 
we shall model the scale dependence of gt{k) with a power 
law, gt {klk^y, and consider three values of the spectral index, 
namely ^ = 0, 1, and 2. Here, ko is a pivot wavenumber, which 
we set to 0.05 Mpc *. Given the anisotropic power spectrum of 
Eq. (164), the CMB temperature and polarization fluctuations 
have the following expectation (Ma et al., 201 1): 
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where = {2In) f k^dk A^fk) A'’^',(k) Po(k)(klko)‘>, and A[(k) 
and A^(k) denote the temperature and E-mode radiation transfer 
functions, respectively. 

The analysis is made using the foreground-cleaned CMB 
temperature maps SMICA, NILC, SEVEM, and Commander, 
where we apply the common mask from the previous Planck re¬ 
lease (Planck Collaboration XII 2014), which has more conser¬ 
vative foreground masking than the newly available mask. The 
implementation details of the optimal estimator can be found in 
the appendix of Planck Collaboration XVI (2015). We estimate 
gLM from the data at the multipoles 2 < ( < 1200. The range of 
multipoles is chosen such that the effect of residual foregrounds 
on the analysis is insignificant. For the assessment, we use re¬ 
alistic simulations containing residual foregrounds. Apart from 
that, we estimate the statistical uncertainty of g^M with various 
f’max values using simulations. We find that the temperature data 


Table 16. Constraints on g*, the quadrupolar asymmetry 
amplitude, determined from the SMICA, NILC, SEVEM, and 
Commander foreground-cleaned maps. Error bars are 68 % CL. 


g, X 10^ ^ = 0 q = i n - 2 


Commander . . . -0.09)|g3 -0.27)| [q 

NTT.r r)(sn+>''3 (^]^+i.5i _f^f^4+o.73 

. ^•'^'■'- 1.62 - 0.99 - 0.66 

SEVEM. 0.13+];^^ 0.031];^ -0.01)“;^, 

SMICA. 0.23!];™ 0.16!];^ 0-13!^;“] 


at f > 1500 make little contribution to the reduction of statistical 
uncertainty, due to the point-source masking of numerous holes. 

Table 16 shows the constraints on g* obtained after marginal¬ 
izing over the direction E^i. Our limits provide a stringent test of 
rotational symmetry during inflation, with g* = 0.23j^[ 24 x 10 “^ 
(68 % CL) for the scale-independent case (i.e. q - 0),^^ using the 
SMICA temperature map data, where the known effects of sys- 
tematics including the asymmetry of the Planck beams are taken 
into account. We note that the central values for g* in Table 16 
are closer to zero than expected given the errors. The reason for 
this is not clear, but it might simply be due to chance. 

Unlike the results from the temperature data, the analysis of 
foreground-cleaned polarization data shows a large, apparently 
highly statistically significant departure from a null value; how¬ 
ever, this result is likely due to systematic errors, most notably 
a temperature-to-polarization leakage associated with the mis¬ 
match of elliptical beams, for which at present a reliable estimate 
is lacking. 


13. Combination with BiCEP2/Keck Array-Planck 
cross-correiation 

We discuss the implications of the recent constraints on the pri¬ 
mordial B-mode polarization from the cross-correlation of the 
BICEP2 and Keck Array data at 150 GHz with the Planck maps 
at higher frequencies to characterize and remove the contri¬ 
bution from polarized thermal dust emission from our Galaxy 
(BICEP2/Keck Array and Planck Collaborations (2015), here¬ 
after BKP). On its own, the BKP likelihood leads to a 95 % CL 
upper limit r < 0 . 12 , compatible with and independent of the 
constraints obtained using the 2015 Planck temperature and 
large angular scale polarization in Sect. 5 (note, however, that 
the BKP likelihood uses the Hamimeche-Lewis approximation 
(Hamimeche & Lewis, 2008) which requires the assumption of 
a fiducial model). The BKP results are also compatible with the 
Planck 2013 Results (Planck Collaboration XVI, 2014; Planck 
Collaboration XXII, 2014). The posterior probability distribu¬ 
tion for r obtained by BKP peaks away from zero at r 0.05, 
but the region of large posterior probability includes r -0. 

Here we combine the baseline two-parameter BKP likeli¬ 
hood using the lowest five B-mode bandpowers with the Planck 
2015 likelihoods. The two BKP nuisance parameters are the B- 
mode amplitude and frequency spectral index of the polarized 
thermal dust emission. The combined analysis yields the follow¬ 
ing constraint on the tensor-to-scalar ratio: 

ro .002 < 0.08 (95 % CL, Planck TTh-IowPh-BKP) , (167) 

The constraints from the Planck 2013 data by Kim & Komatsu 
(2013) should be multiplied by a factor of V2 in our normalization. 

















Planck Collaboration: Constraints on inflation 


55 



i^H Planck TT+lowP 

Planck TT+lowP+BKP 
Planck TT+lowP+BKP+BAO 
Natural inflation 
Hilltop quartic model 
A a attractors 

— - Power-law inflation 

- Low scale SB SUSY 

- inflation 

V 

- V (xcf)^ 

— y ot 
- V QC(f) 

- y OC 02/3 

• At* =50 
m A^*=60 


0.94 


0.96 

Primordial tilt (n„) 


0.98 


1.00 


Fig. 54. Marginalized joint 68% and 95% CL regions for and ro.oo 2 from Planck alone and in combination with its cross¬ 
correlation with BICEP2/Keck Array and/or BAO data compared with the theoretical predictions of selected inflationary models. 


further improving on the upper limits obtained from the different 
data combinations presented in Sect. 5. 

By directly constraining the tensor mode, the BKP likeli¬ 
hood removes degeneracies between the tensor-to-scalar ratio 
and other parameters. Adding tensors and running, we obtain 

'• 0.002 < 0.10 (95 % CL, Planck TTh-IowPh-BKP) , (168) 

which constitutes almost a 50 % improvement over the Planck 
TTh-IowP constraint quoted in Eq. (28). These limits on tensor 
modes are more robust than the limits using the shape of the 
spectrum alone owing to the fact that scalar perturbations 
cannot generate B modes irrespective of the shape of the scalar 
spectrum. 

13.1. Implications of BKP on selected inflationary models 

Using the BKP likelihood further strengthens the constraints 
on the inflationary parameters and models discussed in Sect. 6 , 
as seen in Pig. 54. If we set 63 = 0, the first slow-roll pa¬ 
rameter is constrained to ei < 0.0055 at 95 % CL by Planck 
TTh-IowPh-BKP. With the same data combination, concave po¬ 
tentials are preferred over convex potentials with logZ? = 3.8, 
which improves on logZ? = 2 obtained from the Planck data 
alone. 

Combining with the BKP likelihood strengthens the con¬ 
straints on the selected inflationary models studied in Sect. 6. 
Using the same methodology as in Sect. 6 and adding the BKP 
likelihood gives a Bayes factor preferring over chaotic in¬ 
flation with monomial quadratic potential and natural inflation 
by odds of 403:1 and 270:1, respectively, under the assumption 
of a dust equation of state during the entropy generation stage. 
The combination with the BKP likelihood further penalizes the 
double-well model compared to inflation. However, adding 


Table 17. Results of inflationary model comparison using the 
cross-correlation between BICEP2/Keck Array and Planck. This 
table is the analogue to Table 6, which did not use the BKP like¬ 
lihood. 


Inflationary Model 

In Box 

Wint = 0 Wint t 0 

R + R-/6M^ 


+0.3 

n = 2 

-6.0 

-5.6 

Natural 

-5.6 

-5.0 

Hilltop {p - 2) 

-0.7 

-0.4 

Hilltop {p = 4) 

-0.6 

-0.9 

Double well 

-4.3 

-4.2 

Brane inflation (p = 2) 

+0.2 

0.0 

Brane inflation (p = 4) 

+0.1 

-0.1 

Exponential inflation 

-0.1 

0.0 

SB SUSY 

-1.8 

-1.5 

Supersymmetric or-model 

-1.1 

+0.1 

Superconformal (m = 1) 

-1.9 

-1.4 


BKP reduces the Bayes factor of the hilltop models compared 
to R^, because these models can predict a value of the tensor-to- 
scalar ratio that better fits the statistically insignificant peak at 
r 0.05. See Table 17 for the Bayes factors of other inflationary 
models with the same two cases of post-inflationary evolution 
studied in Sect. 6. 

13.2. Implications of BKP on scalar power spectrum 

The presence of tensors would, at least to some degree, require 
an enhanced suppression of the scalar power spectrum on large 
scales to account for the low-/” deficit in the spectrum. We 
therefore repeat the analysis of an exponential cut-off studied 
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Fig. 55. Impact of BKP likelihood on scalar primordial power 
spectrum reconstruction. We show how including the BKP like¬ 
lihood affects the reconstruction in Sect. 8.3. The top panel is to 
be compared with the reconstructions in Fig. 27, and we observe 
that including BKP has a minimal impact given the uncertainty 
in the reconstruction. The middle panel is to be compared with 
Fig. 31, and here we notice that including BKP excludes the tra¬ 
jectories with large values of e. The bottom panel shows how 
the inflationary potential reconstructions are modified by BKP 
(to be compared with Fig. 32). 


in Sect. 4.4 with tensor perturbations included and the standard 
tensor tilt (i.e., - -r/8). Allowing tensors does not signif¬ 

icantly degrade the improvement found in Sect. 4.4 for 
Planck TTh-IowP with a best fit at r « 0. When the BKP like¬ 
lihood is combined, we obtain Ax^ - -4 with respect to the 
base ACDM model with a best fit at r 0.04. However, since 
this model contains 3 additional parameters, it is not preferred 
over the base ACDM. 

In Fig. 55 we show how the scalar primordial power spec¬ 
trum reconstruction discussed in Sect. 8.3 is modified when 
the BKP likelihood is also included. While the power spec¬ 
trum reconstruction hardly varies given the uncertainties in the 
method, the trajectories of the slow-roll parameters are signif¬ 
icantly closer to slow roll. When the 12-knot reconstruction 
is carried out, the upper bound on the tensor-to-scalar ratio is 
r < 0.11 at 95 % CL. The per degree of freedom for the 5 
low-k and 6 high-k knots are 1.14 and 0.22, respectively, corre¬ 
sponding to p-values of 0.33 and 0.97. 

13.3. Relaxing the standard single-field consistency 
condition 

We now relax the consistency condition (i.e., Wt = -^/8) and 
allow the tensor tilt to be independent of the tensor-to-scalar ra¬ 
tio. This fully phenomenological analysis with the BKP likeli¬ 
hood is complementary to the study of inflationary models with 
generalized Lagrangians in Sect. 10, which also predict mod¬ 
ifications to the consistency condition n^ = -r/8 for a nearly 
scale-invariant spectrum of tensor modes. In this subsection we 
adopt a phenomenological approach, thereby including radical 
departures from < 0, including values predicted in alterna¬ 
tive models to inflation (Gasperini & Veneziano, 1993; Boyle 
et ah, 2004; Brandenberger et ah, 2007). In Sect. 10 we fold in 
the Planck /nl constraints (Planck Collaboration XVII, 2015), 
whereas here we consider Planck and BKP likelihoods only. 
Complementary probes such as pulsar timing, direct detection of 
gravitational waves, and nucleosynthesis bounds could be used 
to constrain blue values for the tensor spectral index (Stewart 
& Brandenberger, 2008), but here we are primarily interested in 
what CMB data can tell us. 

We caution the reader that in the absence of a clear detec¬ 
tion of a tensor component, joint constraints on r and Wt depend 
strongly on priors, or equivalently on the choice of parameteri¬ 
zation. Nevertheless, the BKP likelihood has some constraining 
power over a range of scales more than a decade wide around 
k 0.01 Mpc \ so the results are not entirely prior driven. 

The commonly used (r, nd parameterization suffers from 
pathological behavior around r = 0, which could be prob¬ 
lematic for statistical sampling. We therefore use a parame¬ 
terization specifying r at two different scales, (analo¬ 

gous to the treatment of primordial isocurvature in Sect. 11) as 
well as the more familiar (r, nd parameterization. We present 
results based on kj = 0.002 Mpc^’ and ^2 = 0.02Mpc“\ 
also quoting the amplitude at k = 0.01 Mpc ' for both pa- 
rameterizations. This scale is close to the decorrelation scale 
for (r, nd for the Planck h-BKP joint constraints. We obtain 
^ 0.002 < 0.07 (0.06) and ro .02 < 0.29 (0.31) at 95 % CL from the 
two-scale parameterization with Planck TTh-IowPh-BAOh-BKP 
(TT,TE,EEh-1owPh-BKP). Eigure 56 illustrates the impact of 
the BKP likelihood on the one-dimensional posterior prob¬ 
abilities for these two parameters. The derived constraint at 
k = 0.01 Mpc-' is rooi < 0.12 (0.12) at 95% CL with 
Planck TTh-IowPh-BAOh-BKP (TT,TE,EEh-1owPh-BKP). The up¬ 
per panel of Fig. 57 shows the relevant 2D contours for the 
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Fig. 56. Posterior probability density of the tensor-to-scalar ra¬ 
tio at two different scales. The inflationary consistency rela¬ 
tion is relaxed and ro .002 ro .020 are used as sampling pa¬ 

rameters, assuming a power-law spectrum for primordial ten¬ 
sor perturbations. When the BKP likelihood is included in the 
analysis, the results with Planck TTh-IowPh-BAO and Planck 
TT,TE,EEh-1owP coincide (red dashed and red solid curves, re¬ 
spectively). 
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Fig. 57. 68 % and 95 % CL constraints on tensors when 
the inflationary consistency relation is relaxed, with Planck 
TT-pIowP-pBAO (blue dashed contours) and TT,TE,EE-p1owP 
(blue shaded regions). The red colours are for the same data plus 
the BKP joint likelihood. The upper panel shows our indepen¬ 
dent primary parameters ro .002 and ro .020 ■ The lower panel shows 
the derived parameters n, and ro.oi. The scale k - 0.01 Mpc^' is 
near the decorrelation scale of (n,, r) for the Planck -pBKP data. 
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Fig. 58. The same as Fig. 57 lower panel, but using n, and ro .002 
as primary parameters. 

tensor-to-scalar ratios at the two scales and the improvement due 
to the combination with the BKP likelihood. The lower panel 
shows the 2d contours in (ro.oi, nt) obtained by sampling with the 
two scale parameterization. Figure 58 shows the 2D contours in 
(ro.oi, nt) obtained by the (ro. 002 , nt) parameterization. 

We conclude that positive values of the tensor tilt nt are 
not statistically significantly preferred by the BKP joint mea¬ 
surement of B-mode polarization in combination with Planck 
data, a conclusion at variance with results reported using the 
BICEP2 data (Lewis, 2014; Gerbino et al., 2014). However, the 
now firmly established contamination by polarized dust emis¬ 
sion could explain the discrepancy. Values of tensor tilt consis¬ 
tent with the standard single-field inflationary consistency rela¬ 
tion are compatible with the Planck -i-BKP constraints. 

14. Conclusions 

The Planck full mission temperature and polarization data are 
consistent with the spatially flat base ACDM model, whose per¬ 
turbations are Gaussian and adiabatic with a spectrum described 
by a simple power law, as predicted by the simplest inflationary 
models. For this release, the basic Planck results do not rely on 
external data. The first Planck polarization release at large an¬ 
gular scales from the LEI 70 GHz channel determines an optical 
depth of T = 0.067 + 0.022 (68 % CL, Planck low multipole 
likelihood), a value smaller than the previous Planck 2013 re¬ 
sult based on the WMAP9 polarization likelihood as delivered 
by the WMAP team. This Planck value of r is consistent with an 
analysis of WMAP9 polarization data cleaned for polarized dust 
emission using the Planck 353 GHz data (Planck Collaboration 
XV, 2014; Planck Collaboration XI, 2015). The estimates of cos¬ 
mological parameters from the full mission temperature data and 
polarization on large angular scales are consistent with those of 
the Planck 2013 release. The TE and EE spectra at ^ > 30 to¬ 
gether with the lensing power spectra lead to cosmological con¬ 
straints in agreement with those obtained from temperature. 

The Planck full mission temperature and large angular scale 
polarization data rule out an exactly scale-invariant spectrum of 
curvature perturbations at 5.6 cr. For the base ACDM model, the 
spectral index is measured to be Us = 0.965 + 0.006 (68 % CL, 
Planck TTh-IowP). No evidence for a running of the spectral in¬ 
dex is found, with dnttldlnk - -0.008 + 0.008 (68 % CL, Planck 
TT+lowP). 














58 


Planck Collaboration: Constraints on inflation 


The Planck full mission data improve the upper bound on 
the tensor-to-scalar ratio to ro .002 < 0.10 (95% CL, Planck 
TT+lowP), a bound that changes only slightly when including 
the Planck lensing likelihood, the high-f polarization likelihood, 
or the likelihood from the WMAP large angular scale polar¬ 
ization map (dust-cleaned with the Planck 353 GHz map). We 
showed how the low-f deficit in temperature contributes to the 
Planck upper bound on ro. 002 , but this deficit is not a statistically 
significant anomaly within the base ACDM cosmology. Using 
the full mission Planck data, we find the upper bound on ro .002 
stable, even when extended cosmological models or models with 
CDM isocurvature are considered. The Planck bound on ro .002 is 
consistent with the recent result ro .002 < 0.12 at 95 % CL ob¬ 
tained by the BKP cross-correlation analysis which accounts for 
contamination by polarized dust emission (Planck Collaboration 
XXX, 2014). By combining Planck TT-blowP with the BKP 
cross-correlation likelihood, we obtain ro .002 < 0.08 at 95 % CL. 

The increased precision of the Planck full mission data re¬ 
duces the area enclosed by the 95 % confidence contour in the 
(Hs, r)-plane by 29 %. We performed a Bayesian model compari¬ 
son with the same methodology as in PCI 13, taking into account 
reheating uncertainties by marginalizing over two extra param¬ 
eters: the energy scale at thermalization, pth, and the parame¬ 
ter Wint characterizing the average equation of state between the 
end of inflation and thermalization. Among the models consid¬ 
ered using this approach, the inflationary model proposed by 
Starobinsky (1980) is the most preferred. Due to its high tensor- 
to-scalar ratio, the quadratic model is now strongly disfavoured 
with respect to inflation for Planck TT-nlowP in combination 
with BAG data. By combining wth the BKP likelihood, this trend 
is confirmed, and natural inflation is also disfavoured. 

We reconstructed the inflaton potential and the Hubble pa¬ 
rameter evolution during the observable part of inflation using a 
Taylor expansion of the inflaton potential or H{4>). This analysis 
did not rely on the slow-roll approximation, nor on any assump¬ 
tion about the end of inflation. When higher-order terms were 
allowed, both reconstructions led to a change in the slope of the 
potential at the beginning of the observable range, thus better fit¬ 
ting the low-^ temperature deficit by turning on a non-zero run¬ 
ning of running and accommodating ro .002 ~ 0.2. These models, 
however, are not significantly favoured compared to lower order 
parameterizations that lead to slow-roll evolution at all times. 

Three distinct methods were used to reconstruct the primor¬ 
dial power spectrum. All three methods strongly constrain devia¬ 
tions from a featureless power spectrum over the range of scales 
0.008Mpc^' < k < 0.1 Mpc^'. More interestingly, they also in¬ 
dependently find common patterns in the primordial power spec¬ 
trum of curvature perturbations V-Rik) at k < 0.008 Mpc '. These 
patterns are related to the dip at f 20-40 in the temperature 
power spectrum. This deviation from a simple power-law spec¬ 
trum has weak statistical significance due to the large cosmic 
variance at low (. 

This direct reconstruction of the power spectrum is com¬ 
plemented by a search for parameterized features in physically 
motivated models. The models considered range from the mini¬ 
mal case of a kinetic energy dominated phase preceding a short 
inflationary stage (with just one extra parameter), to a model 
with a step-like feature in the potential and in the sound speed 
(with five extra parameters). As with the Planck 2013 nomi¬ 
nal mission data, these templates lead to an improved fit, up to 
lS.x^ ~ 12. However, neither the Bayesian evidence nor a fre- 
quentist simulation-based analysis shows any statistically signif¬ 
icant preference over a simple power law. 


We have updated the analysis that combines power spectrum 
constraints with those derived from the /nl parameters (Planck 
Collaboration XVII, 2015). New limits on the sound speed in¬ 
ferred from the full mission temperature and polarization data 
further constrain the slow-roll parameters for generalized mod¬ 
els, including DBI inflation. For the first time, we derived com¬ 
bined constraints on Galileon inflation, including the region of 
parameter space in which the predicted spectrum of gravitational 
waves has a blue spectral index. 

Several models motivated by the axion monodromy mecha¬ 
nism in string theory predict oscillatory modulations and corre¬ 
sponding non-Gaussianities, potentially detectable by Planck. A 
TT -only analysis picks up four possible modulation frequencies, 
which remain present when the high-^ polarization likelihood is 
included. An inspection of frequency residuals in the high-^ TT 
likelihood does not reveal evidence of foreground-related sys- 
tematics at similar frequencies. However, a Bayesian evidence 
comparison prefers the smooth base ACDM model over modu¬ 
lated models, suggesting that the latter could simply be fitting the 
noise in the data. The monodromy model predicts resonant non- 
Gaussian features correlated to power spectrum features. A par¬ 
tial analysis beyond the power spectrum was presented. We also 
constrained a possible pseudo-coupling of the axion to gauge 
fields by requiring that non-Gaussianities induced by inverse de¬ 
cay satisfy the Planck bounds on /nl- 

Section 11 reports on a search for possible deviations from 
purely adiabatic initial conditions by studying a range of mod¬ 
els including isocurvature modes as well as possible correlations 
with the adiabatic mode. The Planck full mission temperature 
data are consistent with adiabaticity. The Planck TT data place 
tight constraints on three-parameter extensions to the flat adi¬ 
abatic base ACDM model, allowing arbitrarily correlated mix¬ 
tures of the adiabatic mode with one isocurvature mode (of ei¬ 
ther the CDM, baryon, neutrino density, or neutrino velocity 
type). Adding the high-f TE and EE polarization data further 
squeezes the constraints, since polarization spectra contain ad¬ 
ditional shape and phase information on acoustic oscillations. 
The likelihood with polarization included is in agreement with 
adiabatic initial conditions. However, the tightening of the con¬ 
straints after including polarization must be interpreted with cau¬ 
tion because of possible systematic effects. For this reason we 
emphasize the more conservative Planck TT-rlowP bounds in 
Table 15. The constraints on the six base-ACDM cosmologi¬ 
cal model parameters remain stable when correlated isocurva¬ 
ture modes are allowed. The largest shifts occur for the neutrino 
density mode, but these shifts are not significant (i.e., < 1 cr). 
The constraints on the tensor-to-scalar ratio also remain stable 
when isocurvature modes are allowed. 

Finally we examined the connection between inflation and 
statistical isotropy, a key prediction of the simplest inflation¬ 
ary models. We tested separately the two lowest moments of an 
anisotropic modulation of the primordial curvature power spec¬ 
trum. We found that a modulated curvaton model proposed to 
explain the observed large-scale dipolar power asymmetry can¬ 
not account for all of the asymmetry, and hence is not preferred 
over statistically isotropic base ACDM. The full mission temper¬ 
ature data place the tightest constraints to date on a quadrupolar 
modulation of curvature perturbations. 
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